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Introduction 


Nuclear  explosion  seismology,  as  a  result  of  the  international 
consequences  of  nuclear  testing,  has  a  strong  political  compo- 
nenL  The  time  period  1987-19%  was  marked  by  increased  po¬ 
litical  interest  in  this  subject  as  well  as  rapid  scientific  develop¬ 
ments.  Joint  Verification  Experiments  (JVE’s)  between  the 
Soviet  Union  and  the  United  States  provided  new  opportunities 
for  joint  testing  of  seismic  instrumentation  and  capabilities.  A 
new  set  of  protocols  were  exchanged  between  the  U.S.  and  the 
Soviet  Union  concerning  the  enforcement  of  the  1974  Thresh¬ 
old  Test  Ban  Treaty  (TTBT)  and  1976  Peaceful  Nuclear  Ex¬ 
plosion  Treaty  (PNEri-  These  and  other  political/scientific  de¬ 
velopments  are  disctissed  in  the  first  main  seaion  of  this  review. 

Research  has  continued  on  verification  issues  including  detec¬ 
tion,  location,  discrimination  and  yield  determination.  These  is¬ 
sues  are  addressed  in  the  second  seaion.  Much  emphasis  was 
placed  on  the  utilization  of  arrays  in  deteaion  and  location 
problems.  Yield  determination  research  has  focused  upon  the 
Lf  phase  as  a  stable  indicator  of  yield  as  well  as  ways  to  com¬ 
bine  multiple  yield  estimates  from  a  single  explosioa 

Explosion  source  theory  is  the  third  area  of  research  that  is 
reported.  Work  has  been  divided  primarily  between  teleseisraic 
and  regional  studies  although  a  small  amount  of  research  has 
continued  with  close-in  data  extending  into  the  hydrodynamic 
regime.  A  number  of  researcher  have  begun  to  quantify  the 
secondary  source  process,  spall  (.ensile  failure  of  near  surface 
layers  above  the  explosion)  as  to  its  contribution  to  regional 
seismograms.  As  identified  in  the  Office  of  Technology  As¬ 
sessment  (OTA)  Report  on  Seismic  Verification  (1988),  dis¬ 
crimination  difficulties  may  develop  for  low-yield  lest  ban 
treaties  with  large  chemical  explosions  detonated  for  mining 
purposes.  Following  the  work  in  the  early  60's  in  this  area,  a 
number  of  researchers  have  returned  to  study  the  source  signa¬ 
tures  from  chemical  explosions  at  regional  distances. 

Wave  propagation  problems  continue  to  play  a  dominant  role 
in  many  of  the  explosion  studies.  This  subjea  comprises  the 
fourth  seaion.  Regional  waveforms  have  been  studied  with  re¬ 
newed  intensity  in  response  to  the  increased  political  interest  in 
treaty  verificatioa  The  development  of  a  physical  understand¬ 
ing  of  the  propagation  of  the  various  regional  phases,  induding 
Pp,  Pg,  mam'e  P,  Sn,  Lj,  and  Rj^  has  been  the  focus  of  many 
papers. 
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Issues  of  attenuation  at  teleseismic  distances  and  thus  bios  in 
yield  estimation  remain  of  interest  although  the  JVEs  provided 
an  independent  assessment  of  these  effects  at  Shagan  River. 
The  physical  description  of  transverse  motions  from  explosions 
continues  to  be  explored  in  terms  of  teaonic  strain  release  al¬ 
though  new  results  charaaerizing  wave  propagation  in 
anisotropic  materials  offer  an  alternate  explanation  for  these 
observations. 

Political/Scientific  Developments 

The  political  aspects  of  nuclear  seismology  are  no  better  por¬ 
trayed  than  in  the  Reagan  Administrations  claim  that  the  Soviet 
Union  had  ’likely  violations*  of  the  Threshold  Test  Ban  Treaty 
[Arms  G}ntrol  Today,  July/August  1990].  In  an  attempt  to  ad¬ 
dress  questions  associated  with  sehmic  measurements  from 
underground  nuclear  tests,  the  Office  of  Technology  Assess¬ 
ment  (OTA)  in  1988  published  SEISMIC  VERIHCATION  of 
Nuclear  Testing  Treaties.  This  document  culminated  a  na¬ 
tional  review  of  seismic  capabilities  by  government,  private,  and 
academic  geophysicist!  seismologist!  geologist!  statistician! 
and  physidsu.  The  report  identified  the  tasks  in  monitoring 
underground  nudear  explosions  induding  deteaion,  identifica¬ 
tion,  evasion,  and  yield  determinatioa 

It  discussed  Soviet  compliance  and  the  ability  to  verify  the 
TTBT  and  PNET.  The  report  concluded  that  all  yield  esti¬ 
mates  of  U.S.  and  Soviet  lesu  are  within  the  90  percent  confi¬ 
dence  level  for  yields  of  ISO  ki  or  less.  A  number  of  authors  in 
the  past  four  years  have  looked  dosely  at  the  yield  estimation 
question  and  come  to  similar  conclusions  [Evemdun  and 
Marsh,  1987;  Sykes  and  Davi!  1987],  The  report  recom¬ 
mended  a  phased  approach  to  further  beaties  restriaing  or 
eliminating  the  testing  of  nuclear  weapons.  It  was  concluded 
that  a  threshold  between  S  and  10  kt  could  be  reliably  moni¬ 
tored  seismically  with  lower  levels  possible.  Restricted  testing 
praaices  that  limit  detonations  to  below  the  water  table,  avoid 
decoupling  scenario!  and  provide  spedai  handling  for  large 
chemical  explosions  would  make  even  lower  threshold  limits 
possible  with  existing  seismic  Instrumemation.  As  noted  by 
Richards  (1990)  not  all  government  officials  agreed  with  the 
condusions  of  the  OTA  Report  Robert  Barker  a  spokesman 
for  the  Department  of  Defense  stated  ’(the  OTA  Report)  has 
homogenized  fact  with  fiaion.*  Richards  and  Lindh  (1987) 
have  argued  for  the  low-yield  threshold  test  ban  treaty  while 
Hamm  (1987)  has  presented  arguments  against  such  a  testing 
change.  A  further  debate  of  testing  issues  can  be  found  in  the 
discussions  of  Miller  et  at  ( 1987)  and  Feiveson  et  al.  (1987). 

As  noted  in  the  OTA  Report,  the  Federal  Government  has 
continued  to  request  independent  estimaia  of  yields  to  assess 
uncertainties  in  seismic  yield  estimates.  CORRTEX 
(Continuous  Reflcaometty  for  Radius  venus  Time)  is  a  non- 
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seismic  method  suggested  to  fit  this  verification  requirement. 
This  hydrodynamic  methodology  requires  on-site  measure¬ 
ments  and  deduces  the  yield  of  the  explosion  from  the  velocity 
of  the  shock  wave  and  the  properties  of  the  surrounding  media. 
In  a  recent  intetvievi»  [Arms  Control  Today,  July/August  1990] 
Ambassador  C  Paul  Robinson  claims  that  the  U.  S.  experience 
with  CORRTEX  gives  an  error  in  yield  of  12%  which  in  a  new 
geological  environment  may  become  as  large  as  30%.  Lamb 
(1988)  suggests  that  under  certain  drcumsiances  that  these  er¬ 
rors  may  be  larger.  Robinson  attributes  the  seismic  yield  error 
faaors  to  be  in  the  range  of  1.7  to  1.75.  He  goes  onto  remark 
that  prior  to  the  1  June  1990  testing  protocols  exchanged  be¬ 
tween  the  Soviet  Union  and  the  United  States  that  the  U.  S. 
could  not  be  certain  if  the  Soviets  had  tested  over  400  to  SOO 
kilotons.  These  statements  about  seismic  and  CORRTEX 
yield  estimates  were  criticized  by  Richards  and  Lamb  in  Arms 
Control  Today,  September  1990. 

A  set  of  Joint  Verification  Experiments  with  U.  S.  and  Soviet 
cooperation  were  conduaed  in  1988  [Robinson,  1989].  Hy¬ 
drodynamic  measurements  were  made  by  both  countries  at  the 
Nevada  Test  Site  (NTS)  (Anderson  and  Manning,  1988]  and 
the  Soviet  test  site  (Manning  and  Anderson,  1988]  at  Shagan 
River  in  Eastern  Kazakhstan.  Seismic  measurements  of  these 
explosions  in  both  the  U.S.  and  the  Soviet  Union  [Priestley  et 
aL,  1990]  Lndicate  that  hydrodynamic  yields  determined  from 
these  explosions  are  in  good  agreement  with  seismic  yield  esti¬ 
mates  [Sykes  and  Ekstrom,  1989]. 

Seisffiic  experiments  involving  the  Soviet  Union  and  the 
United  States  have  included  installation  of  seismic  stations 
within  the  borders  of  both  nations  [Science,  25  August  1989]. 
This  instrumentation  program  began  in  19^  as  a  cooperative 
enterprise  between  the  Natural  Resources  Defense  Council 
(NRDC)  and  the  Soviet  Union  [von  Hippel,  1989].  It  has  since 
expanded  to  include  the  Incorporated  Research  Institutions  for 
Seismology  (IRIS)  and  the  Federal  government  [Science,  25 
August  1989].  The  IRIS  sponsored  instrumentation  has  been 
entitled  the  Eurasian  Seismic  Studies  Program  (ESSP).  The 
program  is  designed  to  provide  data  for  constraining  the  seis¬ 
mic  properties  of  the  crust  and  continental  lithosphere,  seismic 
propagation  characteristics,  noise  properties,  and  differences  in 
seismic  sources.  Six  broadband  Global  Seismic  Network 
(GSN)  stations  have  been  installed  with  a  target  of  twenty  sta¬ 
tions.  Two  regional  networks  are  planned  for  Caucasus  and 
Kirghizia  with  an  additional  small  aperture  (2  km)  array  at 
Kirghizia  [IRIS  Proposal,  1991-1995].  The  coming  years  will 
see  much  activity  in  the  study  of  crust  and  upper  mantle  struc¬ 
ture  [Priestley  et  al.,  1988],  seismic  noise  (Berger  et  at,  1988], 
seismicity  [Thurber  et  aL,  1989]  and  regional  vtave  propagation 
in  the  Sc^et  Union. 

A  set  of  protocols  was  signed  at  the  summit  meeting  between 
the  Soviet  Union  and  the  United  States  on  June  1, 1990.  These 
led  to  consent  by  the  Senate  of  the  TTDT  (1974)  and  PNET 
(1976)  in  Oaober  of  1990.  The  Treaties  came  into  force  on 
December  1 1, 1990.  The  protocols  establish  on-site  inspeaions 
of  tests  as  well  as  the  use  of  in-country  seismic  stations  for 
monitoring  certain  underground  tests  [Chemical  and  Engi¬ 
neering  News,  October  8,  1990].  These  protocols  were  dis¬ 
cussed  with  Ambassador  C  Paul  Robinson  in  an  interview  in 
Arms  Control  Today  (Juty/August  1990).  The  responsibility 
fbr  seismic  in-country  monitoring  went  to  the  On-Site  Inspec¬ 
tion  Agency.  For  explosions  with  anticipated  yields  above  35 
kilotons  the  U.  S.  has  the  right  for  on-site  inspection  of  the  test 


location.  For  explosions  that  are  planned  to  exceed  50  kilotons, 
CORRTEX  measurements  are  possible.  In  addition,  in-coun¬ 
try  seismic  measurements  are  allowed  for  these  larger  explo¬ 
sions. 

Under  the  United  Nations  sponsored  Conference  on  Disar¬ 
mament  (CD),  the  Group  of  Scientific  Experts  (GSE)  has  up- 
dertaken  a  series  of  experiments  to  investigate  the  exchange 
and  interpretation  of  seismograms  on  an  international  scale 
(Lamb,  1988  and  Tirman,  1988].  These  experiments  are  de¬ 
signed  to  investigate,  develop  and  deploy  tools  necessary  for 
electronic  data  transfer  and  management  and  provide  experi¬ 
ence  that  might  be  necessary  in  an  international  effort  to  moni¬ 
tor  underground  nuclear  testing.  DARPA  has  taken  the  US 
lead  in  this  effort  with  the  establishment  of  an  Experimental  In¬ 
ternational  Data  Center  in  Washington,  DC.  Seismic  stations 
that  will  be  participating  in  the  program  are  Blacksburg,  VA 
Lajitas,  TX;  North  Pole,  AK;  Pinon  Flat  Observatory,  CA; 
Pinedale,  WY;  and  Berkeley,  CA 

Verification 

The  initial  task  of  signal  deteaion  in  any  verification  process 
was  well  reviewed  in  the  OTA  report  Quantification  of  back¬ 
ground  noise  in  Eastern  Kazakhstan  [Berger  et  at,  1988]  and 
the  United  States  [Rodgers  et  al.,  1987]  continue  in  an  effort  to 
improve  the  deteaion  capabilities  of  new  systems.  Given 
(1990)  reports  on  noise  levels  at  the  four  IRIS  installations  in 
the  Soviet  Union.  Renewed  interest  in  high  frequency  signal 
propagation  has  resulted  in  broadband  noise  chatoaenzation 
and  consideration  of  instruments  buried  e:  great  depths 
[Bantow  et  aL,  1990].  As  noted  by  Sertno  (1990)  detection 
capabilities  are  dependent  on  signal  attenuation  as  well  as  noise 
levels.  Detection  capabilities  of  existing  networks  were  con¬ 
cluded  to  be  adequate  for  well  coupled  explosions  with  yields  of 
a  few  kilotons  or  less  in  USSR[OTA  Report].  Small  aperture 
regional  arrays  have  also  been  investigated  as  to  their  daection 
capabilities  as  exemplified  by  the  study  at  NORE3S  [Sereno 
and  Bratt,  1989].  Problems  with  deteaion  of  decoupled  explo¬ 
sions  [OTA  Report;  Glenn  and  Rial,  1987]  which  wxxild  re¬ 
quire  more  extensive  in-country  networks  for  monitoring  to  the 
lowest  yield  levels  were  identified. 

The  second  task  in  any  verification  operation  is  event  loca- 
tioa  Work  in  this  area  has  extended  to  the  utilization  of  single 
station  three-component  data  for  both  event  detection  and 
location  [Magotra  e:  aL,  1987].  Similar  analysis  was  conduaed 
using  a  sparse  array  of  singie  three-component  stations  in  East¬ 
ern  Kazakhstan  [Thurber  et  at,  1989].  These  studies  indicate 
the  strength  of  polarization  analysis  [Nakanishi  and  Jarpe, 
1988]  along  with  multiple  phase  identification  in  conduaing 
earthquake  and  explosion  locations.  Inaeased  reliance  on 
small  arrays  in  verification  has  resulted  in  the  investigation  of 
arrays  as  location  tools.  Bratt  and  Bache  (1988)  discuss  back- 
azimuth  and  arrival  time  determination  using  the  NORESS 
and  FINESS  arrays.  Same  et  at  (1990)  report  on  azimuth  es¬ 
timation  using  the  NORESS  regional  array.  Harris  (1990)  and 
Suteau-Henson  (1990)  compare  direction  estimation  from  high 
frequency  arrays  and  three-component  stations  and  find  that 
the  single  station  estimates  are  more  susceptible  to  errors  in¬ 
troduced  by  low  signal  to  noise  ratios.  Recent  work  has  been 
reported  that  attempts  to  take  explicit  consideration  of  the  two- 
and  three-dimensional  nature  of  the  eanh  in  producing  loca¬ 
tions  [TralU  and  Johnson,  1987;  Nelson  and  Vidale,  1990]. 
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Discriiomatioa  studies  during  the  past  Tour  yean  have  fo¬ 
cused  primarily  on  events  from  the  Western  U.  S.  and  the  So¬ 
viet  Union.  Utilizing  the  four  broadband  stations  operated 
around  NTS  by  Lawrence  Livermore  National  Laboratory 
(LLNL),  Taylor  et  aL  (1988)  investigated  low  (1-2  Hz)  and 
high  (6-8  Hz)  frequency  contributions  to  Pn,  Pg,  and  Lg  as  pos¬ 
sible  earthquakc/explosion  discriminants.  The  earthquakes 
were  found  to  be  richer  in  high  frequencies  than  the  explosions 
although  overburied  explosions  weie  enriched  in  high  frequen¬ 
cies  and  often  misclassifred.  The  authors  note  that  it  is  difficult 
to  separate  the  effects  of  a  depth  dependent  Q  model  resulting 
in  stronger  attenuation  in  the  shallower  explosions  from  a  true 
difference  between  explosion  and  earthquake  sources.  Chaet 
(1988)  also  reports  enriched  high  frequencies  from  Western 
U.S.  earthquakes  compared  to  explosions  above  10  Hz,  in  con- 
tradiaion  to  the  claims  of  Evemden  et  al.  (1986).  Denny  et  aL 
(1987)  investigated  the  Ms/mb  discriminant  in  the  Western  U. 
S.  and  found  its  success  was  dependent  on  the  equation  used 
for  determining  magnitude.  A  comprehensive  discrimination 
study,  again  with  data  from  the  LLNL  array,  reported  by 
Taylor  et  aL  (1989).  AO  discriminants  suggested  by  Pomeroy  et 
ai  (1982)  were  tested  in  a  systematic  away  against  a  data  set 
consisting  of  233  NTS  explosions  and  130  Western  U8.  eanh- 
quakes.  The  explosions  were  recorded  at  ranges  of  200-100  km 
while  the  earthquakes  were  observed  at  ranges  of  17S-1300  km. 
Multi-station  discrimination  led  to  misclassification  of  up  to  2% 
of  the  explosions  and  4%  of  the  earthquakes,  with  most  prob¬ 
lems  occurring  for  events  below  mb  4.0.  Suteau-Henson  and 
Bache  (1988)  report  on  a  discrimination  study  using  Pq  and  Lg 
recorded  at  NORESS.  They  find  that  spectral  ratios  can  some¬ 
times  separate  explosions  from  earthquakes  but  that  in  general 
the  spectral  ratios  vary  as  much  within  a  class  of  events  as  they 
do  between  classes,  i^umgardt  and  Young  (1990)  report  on  a 
discrimination  study  using  the  NORESS  array  and  small 
sources  ( magnitude  2-3)  in  the  300-S00  km  distance  range.  In 
contrast  to  studies  in  the  Western  United  States  they  find  the 
best  discriminants  to  consist  of  ratios  between  compressional 
and  shear  energy  (Pn/Sn).  At  the  lowest  yields  (1  kt),  the  dis¬ 
crimination  of  nuclear  explosions  from  quarry  explosions  may 
become  a  problem.  A  number  of  researchers  have  begun  to 
study  these  sources  (see  source  theory  discussion).  Baumgardt 
and  Ziegler  (1988),  Smith  (1989)  and  Hedlin  et  aL  (1989, 
1990)  have  suggested  that  spearal  scalloping  may  be  a  useful 
discriminant  for  ripple-fired  quarry  explosions. 

Yield  determination  has  been  one  of  the  most  active  areas  of 
research  in  the  past  four  years,  with  focus  on  application  to  re¬ 
gional  observations  in  both  the  U.  S.  and  the  U.S.S.R.  This 
emphasis  is  a  natural  result  of  the  newly  deployed  stations  and 
data  that  have  been  exchanged  between  the  two  countnes.  In 
the  case  of  regional  observations,  primary  focus  has  been 
placed  upon  the  utilization  of  Lg  amplitudes  as  a  stable  mea¬ 
sure  of  explcjsion  yield.  Nuttli  pioneered  the  Lg  yield  work  and 
reported  on  yield  estimates  for  explosions  in  E.  Kazakhstan 
[Nuttli,  1987]  and  Novaya  Zemtya  [NuttU,  1988].  The  Lg 
method  of  yield  estimation  has  been  scrutinized  by  a  number  of 
other  researchers.  Patton  (1988)  reported  on  a  comprehensive 
study  of  U.  S.  explosions  from  the  NTS  recorded  at  the  LLNL 
broadband  array  around  the  test  site.  He  also  notes  small  scat¬ 
ter  in  these  Lg  yield  estimates.  Har,sen  et  aL  (1990)  report  on 
the  stability  of  log  RMS  Lg  to  0.03  units  for  IRIS,  COSN,  and 
NORSAR  data  from  Shagan  River  explosions.  Refinement  of 
yield  estimates  using  mb  measurements  has  included  the  as¬ 
sessment  of  data  sensoring,  clipping,  and  noise  contamination 
[McLaughlin  et  al.,  1989;  Jih  and  Shumway,  1989].  A  number 


of  authors  have  made  maximum  likelihood  magnitude  esti¬ 
mates  [McLaughlin,  1988].  Scatter  in  mb  measure.mcnts,  re¬ 
sulting  from  either  near-source  or  near-recenfcr  structure, 
have  been  analyzed  and  modeled  [McLaughlin  and  Jih,  1988]. 
Specific  attention  has  been  applied  to  understanding  the  varia¬ 
tion  of  magnitudes  for  explosions  detonated  in  Yucca  Flats  at 
NTS  where  the  systematic  variations  in  magnitude  can  be 
linked  to  die  basin  struaure  [Ferguson,  1988;  McLaughlin  et 
aL  1987].  These  studies  indicate  a  growing  modeling  capability 
for  two-  and  three-dimensional  wave  propagation  and  an  appli¬ 
cation  of  these  tools  to  yield  determination  questions.  Work 
continues  in  other  areas  of  yield  determination  as  exemplified 
by  the  utilization  of  hydrophone  records  to  make  yield  esti¬ 
mates  [McCreety,  1987]. 

Source  Theory 

Although  not  strictly  a  seismic  wave  propagation  problem, 
one  of  the  biggest  source  measurement  developments  in  the 
past  four  years  has  been  the  utilizatioa  of  time  of  arrival  mea¬ 
surements  in  the  hydrodynamic  regioa  around  explosions  to  de¬ 
termine  yield  [OTA  Report].  CORRTEX  measurements 
played  an  important  role  during  the  JVEs  in  1988.  King  et  al., 
1989  report  on  a  detailed  equation  of  state  analysis  that  de¬ 
scribes  the  coupling  of  energy  in  the  hydrodynamic  region.  The 
study  focuses  upon  the  claim  that  enhanced  coupling  can  result 
from  small  cavities,  1-3J  m/kt^.  Their  numerical  results  sug¬ 
gest  that  there  is  no  enhanced  oxipling  prediaed  when  radia¬ 
tion  transport  of  energy  is  included  in  the  calculations.  These 
types  of  studies  can  be  seen  as  an  initial  step  in  accounting  for 
energy  transport  from  explosions.  Such  data,  linked  to  inter¬ 
mediate  stress  observations,  ofrer  opportunity  to  study  devel¬ 
opment  of  the  seismic  source  function  from  nudear  explosions. 
Denny  and  Goodman  ( 1990)  report  a  re-analysis  of  the  Salmon 
and  Sterling  data.  They  argue  that  the  intermediate  stress 
P'gime  data  show  a  two-wave  system  that  yields  a  high  fre¬ 
quency  reduced  velodty  potential  decay.  They  also  find  the 
&imon-Slerling  decoupling  value  to  be  7Z 

In  attempts  to  further  utilize  teleseismic  measurements, 
Murphy  et  aL  (1989)  and  Murphy  (1989)  have  focused  upon 
network  averaged  teleseismic  spectra  for  recovery  of  source  in¬ 
formation,  such  as  yield  and  burial  depth.  They  suggest  that 
the  network  averaging  process,  which  includes  frequency  de¬ 
pendent  station  correaions,  reduces  the  variances  of  speara 
and  emphasizes  source  properties.  The  above  analysis  has  al¬ 
lowed  the  determination  of  scaling  relations  and  source  depth 
efrects.  Der  et  al.  (1987)  .-eport  a  multi-channel  deconvolution 
technique  which  factors  teleseismic  spectra  into  source  and  re¬ 
ceiver  efrects,  again  attempting  to  separate  these  contributions 
to  isolate  the  source.  At  regional  arrays  Der  et  aL  (1990)  sug¬ 
gest  the  application  of  coherent  processing  techniques  to  iden¬ 
tify  interevent  differences  separate  from  intersite  effects.  Wal¬ 
ter  et  aL  (1988)  investigated  regional  P  waves  at  frequencies  as 
high  as  30  Hz  to  identify  explosion  and  earthquake  source  dif¬ 
ferences.  They  report  an  increase  in  spectral  -implitude  with 
moment  from  earthquakes  that  is  consistent  with  an  source 
model  The  explosion  data  they  investigated  showed  greater 
scatter  than  earthquake  data  but  indicated  a  spectral  decay  that 
was  in  excess  of  r  .  As  noted  by  Chael  (1988),  this  effect  could 
be  a  result  of  either  source  or  attenuation  processes. 

Inversion  techniques  designed  to  charaaerize  the  explosion 
source  in  terms  of  moment  tensors  continue  to  be  investigated 
with  emphasis  on  isolating  isotropic  and  deviatoric  components. 
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Faaon  (198Sb)  describes  an  inversion  of  regional  fundamental 
and  higher  mode  surface  waves  to  determine  the  moment  ten¬ 
sor.  He  finds  that  a  spall  correaion  is  necessary.  Johnson 
(1988)  reported  on  second  order  moment  tensor  inversions  of 
near-source  data  (1-11  km)  from  two  Pahute  Mesa  nuclear  ex¬ 
plosions.  The  source  was  dominated  by  the  isotropic  moment 
tensor  although  a  secondary,  long-period  contribution  (OJ-0.4 
Mz)  was  identified  in  the  vertical  dipole.  This  secondary  source 
might  be  attributed  to  near  surface  spallation.  Stump  (1987) 
found  a  similar  secondary  source  contribution  in  the  vertical 
dipole  of  the  moment  tensor  when  invening  data  from  chemi¬ 
cal  explosions.  In  this  case,  acceleration  measurements  within 
the  spall  zone  were  able  to  corroborate  the  spall  interpretation. 
In  an  attempt  to  investigate  the  robustness  of  isotropic  source 
estimates,  Vasco  and  Johnson  (1989)  conduaed  an  investiga¬ 
tion  of  extremal  solutions.  Their  results  indicate  that  an 
isotropic  contribution  of  10%  or  more  is  necessary  for  unique 
identification  in  data  with  10%  random  noise  or  2%  latent  het¬ 
erogeneity.  In  a  search  for  other  basis  funaions  for  the  seismic 
source  representation,  Vasco  has  suggested  decomposing  the 
source  into  orthogonal  source  time  funaions  (1989).  He 
[Vasco,  1990]  also  suggests  an  orthogonal  decomposition  of  the 
moment  density  tensor  which  might  be  applicable  to  sources 
with  finite  dimensions. 

As  noted  in  the  verification  discussion,  one  future  discrimina¬ 
tion  task  under  a  Reduced  Threshold  Test  Ban  Treaty  is  the 
discrimination  of  waveforms  generated  by  large  chemical  explo¬ 
sions  and  small  nuclear  explosions.  This  motivation  and  the 
more  basic  desire  to  use  chemical  explosions  as  small  scale 
analogs  of  nuclear  explosions  has  led  to  a  moderate  amount  of 
work  in  charaaerizing  chemical  explosions  as  sources  of  seismic 
waves.  Slump  and  Reinke  (1987)  discuss  the  design  and  im¬ 
plementation  of  relatively  low-cost  contained  che  nical  explo¬ 
sions  for  source  studies.  They  [Stump  and  Reinke,  1988]  go  on 
to  use  multiple  chemical  explosions  to  experimentally  test  linear 
superposition.  Baumgardt  and  Ziegler  (1988)  report  on  re- 
gionai  observations  from  ripple-fired  quarry  explosions  that 
suggest  that  spearal  scalloping  can  be  used  as  a  discriminant 
[Smith,  1989]  for  such  blasts.  Hedlin  et  at  (1989)  report  that 
seismograms  firom  suspeaed  quarry  explosions  have  energy  or¬ 
ganized  in  discrete,  time-independent  frequency  bands  while  a 
single  explosion  calibration  event  shows  no  such  modulation. 

Teaonic  release  is  another  secondary  source  accompanying 
explosions  which  may  contaminate  observations  [Cohee  and 
Lay,  1988].  The  phase  reversal  of  some  Rayleigh  waves  gener¬ 
ated  by  underground  explosions  from  E.  Kazakhstan  has  been 
used  as  evidence  for  teaonic  release.  Day  et  al.  (1987)  repli¬ 
cate  these  observed  phase  reversals  coupled  with  time  delays  on 
the  order  of  a  few  seconds.  Their  two-dimensional  models  in¬ 
clude  tectonic  prestress  which  is  released  during  the  explosion 
process.  Pre-stress  greater  than  100  bars  was  required  to  gen¬ 
erate  synthetics  that  matched  observations. 

The  secondary  source  process  known  as  spall,  which  results 
from  the  tensile  failure  of  near  surface  materials  above  the  ex¬ 
plosion,  has  continued  to  be  investigated  as  a  source  of  inter¬ 
mediate  frequency  (0.1-5  Hz)  regional  waves.  Taylor  and  Ran¬ 
dall’s  (1989)  modeling  of  waveforms  suggest  that  the  spall 
source  may  have  important  regional  contributions.  Patton 
(1990),  in  an  attempt  to  put  better  constraints  on  the  size  of 
the  secondary  spall  source,  has  begun  the  investigation  of  data 
from  within  the  spall  zone  using  it  to  develop  scaling  relations 
for  spall  mass  and  momentum. 


Wave  Propagation  and  Slruaure 

Waves  from  underground  exolosions  are  observed  at  near¬ 
source,  regional,  and  teleseismic  distances.  Regional  and  tele- 
seismic  observations  have  application  to  verification  issues 
while  near-source  observations  provide  the  opponunity  to  re¬ 
solve  questions  about  explosion  source  charaaerization  and 
possibly  unique  separation  of  propagation  and  source  contribu¬ 
tions.  We  review  wave  propagation  studies  that  have  devel¬ 
oped  for  all  three  observation  ranges  from  explosions.  The  re¬ 
gional  and,  to  a  lesser  extent,  the  teleseismic  observations  have 
received  the  most  attention  oy  the  research  community.  These 
studies  are  motivated  by  a  desire  to  produce  physical  models 
that  describe  the  generation  and  propagation  of  seismic  energy 
from  the  explosion  to  the  receiver.  Such  a  physical  basis  must 
be  developed  for  waves  that  are  uses  in  verification  tasks  of 
discrimination  and  yield  determination. 

The  scarcity  of  near-source  research  reflects  the  faa  that  such 
observations  have  not  been  seen  to  play  a  verification  role  and 
further,  that  most  such  data  is  recorded  and  analyzed  by  per¬ 
sons  responsible  for  underground  testing.  Murphy  and  Shah 
(1988)  studied  near-field  Rayleigh  waves  generated  by  atmo¬ 
spheric  explosions.  They  found  that  these  waves  were  strongly 
afreaed  by  the  shallow  shear  wave  struaure  near  the  explosion. 
Leonard  and  Johnson  (1987),  in  an  attempt  to  refine  the  ve¬ 
locity  struaure  at  Pahute  Mesa,  conduaed  one-dimensional 
travel  time  inversions  for  P  and  S  struaure  utilizing  near¬ 
source  observations.  Their  philosophy  was  to  begin  the  devel¬ 
opment  of  simple  one-dimensional  models  that  could  be  re¬ 
vised  as  two-  and  three-dimensional  struaural  information  be¬ 
comes  available.  The  effea  of  random  struaure  in  the  near¬ 
source  region  and  quantification  of  scattering  efrects  was  dis¬ 
cussed  by  Stump  and  Reinke  (1988).  They  argue  that  to  quan¬ 
tify  source  models  with  near-source  data  it  is  necessary  first  to 
separate  stochastic  and  deterministic  wave  propagation  effects. 

Regional  wave  propagation  studies  have  resulted  from:  (1) 
the  desire  to  monitor  explosions  of  all  sizes  throughout  the 
world;  (2)  the  development  of  regional  arrays  which  are  pro¬ 
viding  abundant  regional  data  [Sereno  et  al.,  1987];  (3)  the  ap¬ 
parent  stability  of  the  Lg  phase  as  a  measure  of  yield  [Patton, 
1988;  Hansen  et  aL,  1990];  and  (4)  the  recent  establishment  of 
new  regional  stations  in  the  Soviet  Union  [Priestley  et  aL,  1990]. 
The  primary  regional  phases  that  have  been  studied  include  Pg, 
Pr,  Sr,  Lg,  and  Rg.  Charaaerization  of  regional  propagation 
effects  in  the  Western  U.S.  is  one  area  where  there  has  been 
much  work  since  NTS  provides  an  abundant  set  of  artificial 
sources  [Hough  et  al.,  1989]. 

Chael  (1988)  reports  on  Pg  observations  at  Nelson,  NV  of 
NTS  explosions  and  nearby  earthquakes.  As  mentioned  in  the 
discrimination  discussion,  he  uses  these  observations  to  illus¬ 
trate  that  the  eanhquakes  have  more  energy  above  10  Hz  than 
explosions.  Walter  et  al.  (1988)  also  report  on  regional  P  wave 
observations  from  NTS.  Sereno  and  Given  (1990)  desaioe 
numerical  models  intended  to  replicate  P„  attenuation  in  the 
crust  and  uppermost  mantle  to  frequencies  as  high  as  15  Hz. 
They  find  that  a  frequency  independent  Pa  spreading  is  a  bad 
assumption  for  data  reduaion.  Sereno  et  al.  (1988)  have  used 
Pr  and  Lg  observations  at  NORESS  to  inven  for  both  attenua¬ 
tion  and  source  moment.  Their  approach  al'ows  for  a  quantifi¬ 
cation  of  regional  wave  propagation  effects.  Vogfjord  and 
Langston  (1990)  use  f-k  anal^is  of  regional  seismograms  (50- 
300  km)  recorded  at  NCRESS  to  identify  coherent  phases. 
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Composite  seismograms  are  then  construaed  by  piecing  to¬ 
gether  stacked  records  at  the  phase  velocities  for  which  coher¬ 
ent  energy  is  identified.  Wave  propagation  in  the  oceanic  litho¬ 
sphere  was  investigated  by  Sereno  and  Orcutt  ( 1987)  with  focus 
00  Pn  and  Sn  phases.  They  find  that  the  coda  associated  with 
these  ocean  phases  can  be  explained  in  terms  of  a  leaky  organ- 
pipe  mode  in  the  sediment  layer/oceanic  water  column. 

One  of  the  most  interesting  but  east  understood  regional 
phases  is  Lj.  The  primary  motivation  for  a  physical  under¬ 
standing  of  the  phase  comes  from  its  utilization  as  an  estimator 
of  yield  (NuttU,  1987  and  1988;  Patton,  1988],  Work  continues 
in  an  attempt  to  quantify  this  phase  [Frankel  et  al.,  1990]  and, 
in  particular,  understand  its  attenuation  characteristics 
[Mitchell  and  Hwang,  1987;  Xie  and  Mitchell,  1990],  A  major 
outstanding  problem  in  understanding  Lj  '  he  insuffioent 
excitation  of  the  phase  in  one  dimensional  velooty  models. 
Baumgardt  (1990)  has  reported  on  observational  data 
supporting  the  blockage  and  scattering  of  the  phase  funher 
indicating  the  importance  of  two  and  three  dimensional 
structures  on  this  phase.  Dainty  and  ToksOz  (1990)  also 
investigate  the  scattered  energy  in  Lf  as  well  as  P  phases. 

Shallow  explosions  and  earthquakes  in  the  Eastern  U.  S.  have 
been  found  to  generate  relatively  high  frequency  (0.4-2J  Hz) 
fundamental  mode  Rayleigh  waves  [Kafka  and  Reiter,  1987], 
Kafka  (1990)  has  suggested  that  these  waves  may  be  a  good 
depth  discriminant  in  areas  where  they  are  observed. 

Analysis  and  modeling  of  teleseismic  observations  from  nu¬ 
clear  explosions  has  explored  the  effea  of  multipathing  and 
scattering  on  these  waveforms.  Lay  (1987  a  and  b),  Lay  and 
Welc  (1987),  and  Lynnes  and  Lay  (1989)  have  investigated 
teleseismic  P  waves  and  early  coda  for  evidence  of  near-source 
scattering.  Gupta  and  Blandford  (1987)  suggest  that  spectral 
ratios  between  P  and  P  oxia  can  be  used  to  charaaerize  near- 
source  material  properties.  Gupta  et  aL,  (1990a)  through  a 
broadband  f-k  analy^  at  the  EKA  array  of  U.S.  and  Soviet  ex¬ 
plosions,  quantify  the  effects  of  near-receiver  scattering.  This 
same  group  of  workers  [Gupta  et  aL  1990b|  suggest  that  both 
near-receiver  and  near-source  information  can  be  recovered 
from  the  scattered  waveforms.  Oirmier  ( 1987),  in  a  modeling 
study,  attempted  to  determine  the  effects  a  three-dimensional 
velocity  model  beneath  NTS  with  scale  lengths  of  20-100  km 
and  a  few  percent  velocity  fluctuation  would  have  on  teleseis¬ 
mic  mb.  He  finds  that  such  variations  lead  to  as  much  as  a  fac- 
tor  of  three  variation  in  amplitudes  and  if  considered  could  re¬ 
duce  the  mb  variances  by  25%.  He  aiso  suggests  that  P  coda 
may  be  a  more  stable  measure  if  scattering  in  the  crust  and  up¬ 
per  mantle  is  important.  McLaughlin  and  Jih  (1988)  and 
McLaughlin  et  al.  (1987)  investigated  the  effea  of  near-source 
scattering  on  teleseismic  mb  through  a  series  of  two-dimen¬ 
sional  linear-elastic,  finite-difference  calculations.  For  explo¬ 
sions  at  the  French  Sahara  Test  Site  the  maximum  effea  is 
±0.15  mb  units.  Al  Yucca  Flats  variations  as  large  as  ±030 
CBb  units  were  modelled.  These  variations  can  be  reduced  by 
averaging  over  a  network  of  obseivaiions.  Determination  of 
teleseismic  receiver  funaions,  using  multiple  component  data, 
continues  in  an  attempt  to  further  quantify  recerver  effecu. 
Broadband  receiver  funaions  for  seismic  stations  in  E.  Kaza¬ 
khstan  were  reported  by  Priestley  et  aL,  (1988). 

One  of  the  newest  developments  in  wave  propagation  has 
been  a  number  of  papers  beginning  to  quantify  effects  of 
anisotropic  media  on  seismic  waves  observed  from  explosions. 


This  work  has  developed  out  of  an  emphasis  on  anisotropic  ef¬ 
fects  in  the  seismic  exploniioa  industiy.  A  series  of  papers  by 
Ben-Menahem  (1990)  and  Ben-Menahem  and  Sena  (1990)  re¬ 
port  Green’s  funaions  for  anisotropic  media  with  spedal  con¬ 
sideration  of  SH  waves.  Other  studies  of  explosions  in 
anisotropic  media  have  been  reported  by  Mandal  and  ToksOz 
(1990).  They  indicate  that,  with  a  moderate  amount  of 
anisotropy  around  the  explosion,  SH  waves  might  be  gener¬ 
ated. 

Future  Issues 

The  progress  in  a  field  over  a  time  period  such  as  the  last  four 
years  is  motivated  by  important  issues  that  are  confronting  a 
panicular  area.  In  an  attempt  to  provide  a  glimpse  into  the  fu¬ 
ture  of  nuclear  seismology,  five  outstanding  areas  of  research 
are  identified.  These  include  the  development  of  a  physical  ba¬ 
sis  for  regional  wave  propagation,  the  investigation  of  the  seis¬ 
mic  source  charaaerization  of  small  nuclear  explos.jns,  quan¬ 
tification  of  explosion  phenomenology,  the  exploration  of 
quarry  and  engineering  explosions  as  seismic  sources,  and  uti¬ 
lization  of  automated  signal  processing  tools. 

Physical  Basis  for  Regional  Wave  Prcpagation 

This  topic  is  particularly  iraponant  as  small  yield  explosions 
will  most  likely  be  monitored  at  regional  distances.  A  physical 
understanding  of  the  processes  is  needed  if  techniques  deve.- 
oped  in  one  geographical  area  are  to  be  transported  to  a  new 
area.  Much  of  the  work  that  has  been  completed  to  date  on  re- 
^onal  waves  bos  been  primarily  experimental  This  work  must 
be  extended  to  the  development  of  physical  models  which  can 
be  tested.  This  testing  may  involve  extensive  instrumentation, 
not  unlike  that  suggested  by  PasscaL  i  i  order  to  resoNe  impor¬ 
tant  aspects  of  the  regional  wave  propagation.  Some  key  ques¬ 
tions  to  be  addressed  indude;  (I)  How  is  energy  distributed 
among  the  phases  Pn,  Pj,  Sn,  I^?;  (2)  How  homogeneous  or 
inhomogeneous  is  the  crust  and  upper  mantle  (1D,2D3D)?; 
(3)  How  do  we  charaaerize  this  homogeneity/inhoraogeneity 
experimentally?;  (4)  What  are  the  source  depth  effects  on  each 
of  these  phases?;  (S)  How  does  the  regional  energy  compare  to 
that  observed  locally  and  leleseismically? 

Small  Source  Phvsics 

The  seismological  community  has  spent  a  good  deal  of  its  ef¬ 
forts  and  time  charaaerizing  large  (>100  kt)  nuclear  explo¬ 
sions.  This  focus  has  resulted  from  the  150  kt  TTBT  and 
PNET  and  the  good  signal  to  noise  ratios  resulting  from  such 
events.  These  sources  are  often  deeply  buried  (122m/kt*'^).  A 
number  of  techniques  exist  for  separating  these  events  from  the 
eanhquake  population  as  well  as  determining  yield  ftom  the 
wave  amplitudes  (mb  -  logfYield)  relations).  LJttie  work  has 
been  completed  in  charaaerizing  small  yield  nuclear  explosions 
including  careful  resolution  of  such  basic  questions  as  the  ap¬ 
propriate  slope  for  a  mb  -  Iog(Yield)  relation  and  the  effea  ma¬ 
terial  propenies  have  on  such  relations.  As  a  result  of  their 
shallow  depth,  these  sources  can  be  detonated  in  highly 
porous,  dry  materials  as  well  as  in  cavities  offering  an  opfxsnu- 
mty  for  decoupling.  The  seismological  community  has  spent  lit¬ 
tle  time  studying  such  events  either  experimentally  or  theoreti¬ 
cally.  No  consideration  has  been  given  to  how  one  might  char- 
aae.-ize  an  explosive  site  before  a  test  so  mat  some  of  these  ef- 
feos  might  be  anticipated  and  mitigated.  Some  of  the  ques¬ 
tions  that  arise  are:  (1)  What  are  earthquake/expiosion  source 


6 


STUMP:  Nuclear  Ew-osion  Seismology 


ftmctions  at  snuiU  magnitude?;  (2)  What  coupling  variations 
can  be  expected  in  shallow,  weathered,  dry  materials?;  (3)  How 
effective  it  decoupling  in  cavities  or  previously  explosively  dis¬ 
turbed  materials?;  (4)  What  are  the  effecu  of  source  region  in- 
homogeneity?;  (S)  How  do  we  complete  geophysical/geological 
site  characterization  to  constrain  explosion  coupling? 


dear  explosions.  Little  work  has  been  completed  In  describing 
these  sources  as  generators  of  regional  seismic  energy.  Some 
key  questions  that  must  be  addressed:  (1)  What  are  the  differ¬ 
ent  types  of  industrial  explosions?;  (2)  How  does  each  type 
generate  seismic  waves  (P  and  S  waves)?;  (3)  Do  the  individual 
explosions  interact  linearly  or  nonlineariy?;  (4)  Can  multiple 
explosions  be  discriminated  from  a  single  explosion? 


Physical  processes  by  which  chemical  and  nuclear  explosions 
couple  energy  into  the  earth  are  not  well  quamifieo  from  first 
prindples.  As  one  might  expect,  these  processes  are  highly  de¬ 
pendent  upon  the  material  properties  of  the  media  of  interest. 
In  many  cases  the  complete  data  sets  for  resolving  these  cou¬ 
pling  processes  are  not  available.  Source  related  research 
might  indude:  (1)  Quantification  of  chemical/nuclear  source 
differences;  (2)  Determination  of  the  imponance  of  tensile 
bilure  of  near-surfhce  layers  (spall);  (3)  Development  of  a 
physical  model  for  spall;  (4)  Determination  of  energy  balance  in 
the  three  spatial  dimensions;  (5)  A  quantification  of  energy 
coupling  from  the  hydrodynamic  to  the  linear  regimes. 


Many  industrial  explosions  are  detonated  each  year  for  pur¬ 
poses  ranging  from  excavation  for  construction  to  the  recovery 
of  subsurface  resources.  Each  of  these  sources  generate  seis¬ 
mic  waves  and  must  be  separated  from  the  population  of  nu- 


Arrav  Processlng/Automatie  Data  Analysis 

Increasingly  rigorous  verification  requirements  may  require 
increased  data  processing  and  analysis  capabilities  [Bache  et  ai., 
1990;  Bratt  et  aU  1990).  Automation  of  these  tasks  is  and  will 
remain  an  important  area  of  researah  [Anderson,  1990;  Dysart 
and  Pulli,  1990J. 

Much  work  has  been  completed  in  the  past  four  years  thanks 
to  continuing  support  for  basic  research  in  nuclear  .seismology. 
Many  questions  still  remain,  particularly  those  associated  with 
the  ph]^cal  processes  of  energy  coupling  and  wave  propaga¬ 
tion.  As  new  requirements  are  imposed  under  friture  treaty 
scenarios,  an  understanding  of  these  processes  will  be  neces¬ 
sary  for  developing  refined  yield  determination,  discrimination 
and  detection  tools. 
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Abstract.  Near-source  waveforms  from  explosions  detonated  at  the 
Nevada  Test  Site  (NTS)  are  reviewed.  Dau  are  separated  into  four  types: 
free-field  strong  (FFS),  representative  of  the  region  where  material 
strength  dominates;  free-field  weak  (FFV/),  the  region  where  weakly 
nonlinear  properties  and  transition  (tom  elastic  to  elastic  response  are 
important;  free  surface  spall  (FSS),  where  material  tensile  strength  is 
important;  and  free  surface  elastic  (FSE),  where  most  seismic 
observations  begin.  Dau  from  Pahute  Mesa  (FSS  &  FSE),  Rainier  Mesa 
(FFW  3t  FSE),  and  Yucca  Flau  (FSS  8t  FSE)  are  specifically  consideied. 
Each  of  the  dau  types  is  explored  as  to  iu  resolution  of  important 
physical  processes  in  the  source  region  and  resultant  seismic  radiauoa 
Specific  auention  is  paid  to  the  variability  of  these  motions.  Single  and 
scaled  multiple  explosion  peak  accelerations  from  Pahute  Mesa  and  Yucca 
Flau  show  as  much  as  a  factor  of  6-8  scatter  with  range.  Large  scatter  in 
single  explosion  dau  suggesu  a  propagation  path  effect  while  Urge 
amplitudes  for  a  scaled  explosion  from  below  the  water  table  supports  a 
coupling  difference  between  explosions.  Dau  scatter  decreases  at  long 
periods  as  exemplified  by  long  period  momenu  which  have  a 
multiplicative  error  of  1.49  for  Pahute  Mesa.  Numerical  models  of  body 
and  surface  wave  propagation  in  realistic  one-dimersional  Pahute  N'esa 
model]  indicate  strong  effecu  of  velocity  strucoire  near  the  shot  point  (or 
body  waves  traveling  to  the  free  surface  at  short  offset  (<2  depths  of 
burial).  Synthetic  waveform  difference  between  a  site  specific  model  and 
an  average  model  decrease  with  increasing  source-receiver  offset  or 
increasing  period  where  near-source  surface  waves  are  emphasiaed. 
Comparison  of  frec-field  and  free  surface  dau  bom  the  same  explosion  at 
Rainier  Mesa  supports  significantly  reduced  scatter  in  free-field  data. 
Removal  of  the  weathered  layer  as  a  dominant  effect  in  the  free-field  dau 
can  explain  the  reduced  scatter.  Analysts  of  dau  spanning  the  transition 
from  FSS  to  FSE  regions  indicates  that  explosion  geometry  plays  a 
strong  role  In  the  decay  of  free  surface  data.  These  dau  and  extended 
rellecuviiy  calculations  appropnaie  for  Pahute  Mesa  predict  that  strong 
spall  zone  motions  com:  from  a  region  out  to  a  free  surface  range  just 
beyond  one  depth  of  burial  for  explosions  with  standard  scaled  depths  of 
bunal. 

Introduction 

A  physically  based  understanding  of  the  nuclear  explosion  source 
function  is  necessary  for  improving  discrimtnauon  between  earuhquakes 
and  explosions  as  well  as  refining  yield  estimates  for  explosions.  The 
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physical  model  allows  prcdicuon  and  exuapolation  u>  new  environments 
where  expenmenul  dau  may  not  exist  [Murphy,  1977],  It  provides 
correlauon  between  dau  observed  at  different  ranges  such  as  near-souice, 
regional,  and  teleseismic  [Taylor  and  Rambo.  1990).  An  tmpiricai 
approach  to  problems  of  discriminatirM  and  yield  determination  may  be 
adequate  for  limned  conditions,  but  is  not  acceptable  for  the  wide  range  of 
environments  and  distances  that  exist  on  the  earth.  The  focus  of  this 
study  is  near-source  dau  seu  and  their  utility  in  physical  source  model 
resolution  [Stump  and  Johnson.  1984;  Johnson.  1988).  Emphasis  will 
be  placed  on  the  successes  and  faiiiirts  in  separating  source  and 
propagauon  path  effects.  Near-source  is  loosely  defmed  for  this  suidy  as 
source-receiver  separations  of  hundreds  of  meters  to  tens  of  kilometers. 

Explosions,  unlike  natural  events,  can  be  controlled.  Time  of 
occurrence  and  location  are  determined  by  the  experimenter,  which  allows 
for  design  of  g-ound  motion  expcrimeia  whose  expressed  purpose  is  the 
characterization  of  seismic  source  properties  [Suimp  and  Reinke,  1987], 
Not  only  can  the  elastic  field  be  deimined.  but  development  of  this  field 
from  the  near-source  region  where  ihs  maieral  is  nonlinear  can  be 
documenied  Quantification  of  the  notion  field  in  this  environment 
allows  the  linkage  of  :.eismic  observadons  u>  measuremenis  of  pressure, 
stress,  and  shock  front  propagauoa  in  the  hydrodynamic,  strongly 
nonlinear,  and  weakly  nonlinear  regiona.  Each  mouon  environment  is 
related  m  the  yield  of  the  explosion  and  maienal  in  which  the  explosion  is 
detonated.  Comparison  of  yield  estimiies  made  from  dau  within  each 
region  is  dependent  upon  a  physical  understanding  of  appropriate 
cumulauve  physical  processes. 

Seismic  discnminauon  and  yield  deienninaiion  studies  rely  upon 
regionaJ  and  teleseismic  dau  suice  moss  monitoring  scenarios  are  limited 
to  such  dau  (Bache.  1982;  Pomeroy  et  aL,  1982;  U.S.  Congress.  19881. 
Physical  source  constraints  provided  by  near-source  dau  can  be  used  as 
consistency  checks  against  models  developed  from  regional  and 
teleseismic  dau  where  propagation  path  effects  may  be  more  severe 
(Taylor  and  Rambo.  1990). 

Neor-source  dau  sets  discussed  in  this  paper  are  from  the  Nevada  Test 
Site  (NTS)  where  materials  such  as  loosely  consolidated  alluvium, 
volcanic  tuffs  and  rhyolites,  granites,  and  other  sedimentary  rocks  con  be 
found.  Particular  aueniion  will  be  paid  to  the  documenuuon  of  free 
surface  and  subsurface  ground  motions  recorded  from  a  variety  of 
explosions  ai  NTS.  These  motions  will  be  used  io  constrain  the 
equivalent  elastic  source  funcuon  and  size  of  the  explosion.  Special 
aticntion  will  be  paid  u>  the  spatial  vaiiabiliiy  of  such  motions  and 
resulting  errors  in  source  strength  esomaies. 

This  paper  begins  with  a  division  of  ihe  near-source  motion  fields 
around  the  explosion  followed  by  a  review  of  previous  work. 
Observational  dau  seu  bom  Pahute  Mesa,  Yucca  Fbis.  and  Rainier  Mesa 
are  introduced  and  explored.  Wavefield  modeling  and  source  inversion  are 
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piesoued  as  the  observational  data  are  used  to  constrain  physical  processes 
around  the  explosion.  Implications  and  conclusions  are  given  as 
developed  from  data  analysis. 

Near-Souree  Waves 

Ground  motion  around  nuclear  explosions  can  be  observed  at  the  free 
surface  where  instrumentation  is  relatively  inexpensive  or  at  depth  (free- 
field)  where  the  effects  of  the  free  surface  are  a  minimum.  Both  free-field 
(FF)  and  free  surface  (F5)  data  will  be  discussed  in  this  paper. 

For  the  purposes  of  our  discussion  the  motion  field  around  the 
explosion  is  divided  inur  five  regions  [Rodean.  1971.  1981].  The 
strongest  motions  from  contained  explosions  are  found  in  the  free-field 
hydrodynamic  (FFH)  region  where  shock  waves  propagate.  The  velocity 
of  this  :hock  front  is  used  to  make  yield  estimates  utilizing  the 
CORRTEX  methodology  (King  et  al.,  1989].  Moving  farther  from  the 
explosion  one  encounters  a  region  of  strongly  nonlinear  motions  where 
material  strength  dominates  the  response.  This  region  is  given  the 
nomenclature  FFS  (free-field  strong).  Pressure  and  stress  measurements 
are  typically  used  to  characterize  this  region.  A  weakly  nonlinear  region 
is  encountered  as  the  motion  field  decays  further.  .More  subtle  material 
properties  such  as  the  transition  from  plastic  u>  elastic  response  become 
importanL  For  free-field  dau  we  designate  this  region  as  FFW  with  the 
W  for  weak. 

As  the  motion  field  interacts  with  the  free  surface,  tensile  failure  of 
near-surface  layers  can  result  in  ballistic  motions  or  spall.  This  motion 
environment  is  designated  as  FSS  (free  surface  spall).  The  spall  process 
extends  to  depth  as  the  reflected  tensile  wave  propagates  back  inu)  the 
material.  The  cost  of  fielding  subsurface  gauges  is  often  prohibitive,  so 
most  spall  data  has  been  taken  al  the  free  surface  (Patton,  1990],  As 
waves  propagate  further,  they  transition  to  purely  elastic/anelastic 
motions  which  are  designated  as  free-surface  elastic  (FSE). 

This  paper  will  document  motions  within  the  FFW,  FSS,  and  FSE 
regions  from  NTS  explosions.  Each  of  the  motion  fields  will  be 
illustrated  with  specific  examples  from  Pahute  Mesa  (FSS  ^  FSE), 
Rainier  Mesa  (FFW  <t  FSE).  and  Yucca  Flats  (FSS  A  FSE). 

Previous  Work 

Free-field  and  free  surface  ground  motions  and  stresses  have  been 
measured  since  the  inception  of  underground  testing.  Free-field  motions 
(FFS  A  FFW)  have  been  used  to  assess  the  importance  of  material 
properties  on  strong  ground  motions.  Such  assessments  are  particularly 
important  when  large  engineering  structures  are  placed  close  to  the 
explosion.  The  tunnel  shots  in  Rainier  Mesa  are  an  example  of  this  type 
of  test  .Measurements  are  often  associated  with  these  structures  which 
might  bias  free-field  motions.  Some  of  the  earliest  estimates  of  free-field 
decay  rates  in  alluvium,  tuff,  granite,  and  salt  were  summarized  by 
Wheeler  and  Preston  (1968],  Ferret  and  Bass  (1975]  extended  the  datasets 
to  include  nuclear  detonations  in  dry  alluvium,  dry  uiff,  wet  tuff,  and  hard 
rock  (granite,  salt,  dolomite,  and  other  sedimentary  rocks).  These  studies 
focused  upon  attenuauon  of  waveforms  with  slant  range  scaled  by  the 
cube  root  of  explosive  yield.  Germain  (1986]  has  reworked  the  Pcrrei  and 
Bass  analysis  including  more  recent  data. 

Applicaonn  of  free-field  dau  to  calculation  of  seismic  source  functions 
in  the  form  of  reduced  displacement  potentials  for  a  number  of  materials 
was  made  by  Wenh  and  Herbst  (1963]  followed  by  the  analytic  modeling 
of  Haskell  (1967).  Murphy  and  Bennett  (1979]  review  free-field  seismic 
dau  for  alluvium,  tuff,  dolomite,  sandstone-shale,  and  interbedded  lave 
flows. 

Near-source,  free  surface  ground  motion  gauges  have  typically  been 
placed  at  ranges  within  2  depths  of  bunal  (OOB)  where  recording 
equipment  is  located.  The  pnmary  purpose  of  these  gages  is  in  assessing 
equipment  survivability.  These  gages  are  usually  within  the  spall  zone 


(Patton,  1990]  and  ore  useful  in  characterizing  motions  designated  as  FSS. 
Free-field  decay  rates  are  not  strictly  applicable  since  these  motions 
involve  interaction  with  the  free  surface.  Bemreuter  et  al.  (1970] 
developed  a  set  of  decay  curves  for  free  surface  data.  Patton  (19901  has 
discussed  the  uidizaiion  of  these  spall  waveforms  and  their  decay  rates  in 
constraining  the  secondary  seismic  source  from  spall. 

A  second  set  of  common  free  surface  dau  involves  measurements 
designed  to  quantify  effects  on  populated  areas  surrounding  NTS  (Hays. 
1974;  O'Brien  and  Lahoud.  1982).  These  observations  and  models 
involve  propagation  distances  of  tens  to  hundreds  of  kilometers.  A  linal, 
small  dau  set  exists  for  free  surface  ground  motion  from  just  outside  the 
spall  zone  to  tens  of  kilometers  (Stump  and  Johnson.  1984;  Johnson, 
1988].  These  dau  sets  have  been  used  to  constrain  the  equivalent  seismic 
source  (Helmberger  and  Hadley.  1981].  Such  measurements  are  less 
expensive  and  easier  to  field  than  free-field  measurements. 

This  work  is  designed  to  compare  and  contrast  the  FFW.  FSS,  and 
FFE  dau  sets  for  purposes  of  constraining  different  parts  of  the  explosive 
source  function.  Dau  sets,  which  have  been  separately  studied  in  the 
past,  are  investigated  as  an  integrated  set.  Each  contains  different 
information;  for  example  FSS  is  dominated  by  interaction  of  the  upgoing 
stress  wave  with  the  free  surface.  An  important  part  of  this  discussion  is 
the  identification  of  the  variability  of  each  dau  set  The  goal  is  to 
quantify  how  each  different  type  of  rear-source  dau  censinuns  the  seismic 
souree  function. 

Observational  Dau 

Dau  from  three  areas  of  NTS  are  reviewed  including  Pahute  .Mesa, 
Yucca  Flats,  and  Rainier  Mesa.  In  ui  attempt  to  emphasize  the  different 
physical  processes  leading  to  these  motions  the  dau  is  divided  into  that 
represented  by  weakly  nonlinear  motions  in  the  free-deid  (FFW),  spall 
motions  at  the  free  surface  (FSS)  and  eiasuc  motions  from  the  free  surface 
(FSE). 

PaMuu  Mtsa  (rSS  <i  FZZ) 

Los  Alamos  National  Laboratory  (LANL)  routinely  makes  three- 
component,  free  surface  digital  acceleration  measurements  within  the  spall 
zone  of  contained  nuclear  explosions.  A  typical  array  of  accelerometers 
for  the  AMARILLO  explosion  with  accompanying  acceleradon  records  is 
given  in  Figure  1.  Spall  zone  accelerograms  are  characterized  by  an  iraiial 
compressive  wave  (positive  xceleraiion).  -I  %  dwell  during  free-fall 
(ballistic  motion)  following  tensile  failure,  and  a  large  spall  clusure  signal 
(positive  acceleration).  Dwell  times  for  muJiiple  kiloton  explosions  are 
near  I  s  as  iUusirared  by  the  AMARILLO  ground  zero  (GZ)  record  with  a 
dwell  lime  of  0.65  s.  The  spatial  effects  of  spall  are  illustrated  in  Figure 
I  where  acceleration  records  at  free  surface  ranges  of  15  m,  238  m,  and 
640  m  are  reproduced.  Spall  initiation  time,  idenufied  by  the  sun  of  -I  j 
dwell,  increases  with  range  while  spall  rejoin,  oesignaied  by  the  impulsive 
rejoin  signal,  decreases  with  range.  Dwell  time  decreases  from  0,65  to 
0.23  s  over  the  observational  range.  As  a  result  of  the  temporal  and 
spatial  finiieness  of  the  spall  secondary  source,  its  equivalent  source 
function  is  peaked  in  the  frequency  domain  (Day  et  al.,  1983;  Slump. 
1983], 

Free  surfxe,  peak  velocity  measurements  from  FSS  are  summarized  in 
Figure  2  where  they  are  plotted  against  scaled  free  surface  range  (r  » 
R/kt'^).  The  dau  are  derived  through  numerical  integration  from  the 
otigtnal  accelerograms.  Six  nuclear  explosions  (Table  I)  from  Pahute 
Mesa  are  represented  in  the  plot  Dau  such  as  that  displayed  in  Figure  2 
can  be  plotted  against  free  surface  ran^e  as  done  in  the  figure  or  slant 
range  as  suggested  by  Ferret  and  Bass  (1973)  or  Paiuan  (1990],  If 
propagation  path  effects  dominaie  the  mouon  field,  then  slasts  range  is  the 
more  appropriate  variable.  If  free  surface  inieiaclion  is  more  important 
then /re*  surface  range  may  be  more  appropriaie.  As  source-receiver  range 
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Fig.  1.  Three,  free-suriace  venical  acceleration  records  from  the  nuclear 
explosion  AMAREXO.  The  slant  ranges  (GM4-912ni.  GM2-696m. 
GMl-640m)  and  free  surfxe  ranges  &om  ground  ato  (CM4.13m,  CM2- 
23Sm,  GM6-6i0ni)  are  given.  The  horizontal  time  scale  is  2.S7  s.  The 
initial  compressive  wave  reaches  each  successive  station  at  greater  time 
while  spall  rejoin,  indicated  by  the  impulsive  secondary  signal,  occurs  a< 
earlier  times  for  the  farther  stations.  The  surface  geometry  of  the 
accelerometer  array  is  given  at  the  booom  of  the  figure. 

increases  relative  us  source  depth,  the  two  distance  measures  converge. 
The  data  are  displayed  as  a  function  of  frtt  lurfact  nuigt  in  this  case  smce 
this  form  emphasizes  the  spall  zone  which  commonly  eaiends  us  2-} 
DOB  range  {244-3dd  mAt*^)  on  the  sealed  ptai.  This  represenution 


Tablet.  Event Chancterisiics 


Event  nanie 

Date 

KEIS  mh 

BACXfiEACH 

780*11 

3.3 

SHEEPSHEAD 

790926 

3.6 

NEBBKXO  ■ 

82062* 

3.6 

CHANCELLOR 

830901 

3.* 

CYBAR 

860717 

5.7 

AMARILLO 

890627 

4.9 

Free  Surface  Range  (nO 

Fig.  2.  Peak  vertical  velocity  from  the  Pahuie  Mesa  spall  zone  dau 
(Table  1)  is  plotted  against  scaled  free  surface  range.  The  free  surface 
range  is  scaled  by  the  cube  root  of  die  ezplosion  yield  while  the  peak 
velocity  values  are  unsealed.  Power  law  models  were  fit  to  the  dau  for 
ranges  less  tlun  and  greater  than  lOO  nvta^^.  This  free  surface  range 
corresponds  ns  appioaimaiely  one  scaled  depth  of  burial,  122  for 

NTS  explosions. 

emphasizes  the  strong  irapea  geometry  has  on  spall  zone  data.  From  CZ 
to  nfrtt  surfact  range  of  apposimaielv  1  DOB,  veioaiy  (acceleration  and 
displacement  also)  decays  little  far  veloeity  dau)  followed  by  a 
faster  decay  it  farther  ranges  This  central  region  with  little 

spatial  decay  indicaies  that  out  to  tfrn  oufaet  range  equal  to  one  DOB 
there  is  liule  change  in  slant  range.  The  implicauon  of  such  decay 
patterns  is  that  the  central  portion  of  the  mall  zone  will  have  nearly 
constant  escape  velocities  or  momemuffl  contnuujon  U]  a  free  surface 
range  equal  U3  I  DOB.  Beyond  this  range  ui  the  edge  of  the  spall  zone, 
mouons  decay  rapidly,  thus  reducing  the  spall  eoniribuuon  from  this  area. 
A  second  implication  of  this  decay  pattern  is  that  sating  variations  or 
uymmeuies  in  the  edge  of  the  ipall  zone  found  from  overhead 
photography  of  die  peaces  may  nm  be  u  imponam  in  the  central  region 
where  mouons  decay  liule  P^alker,  1982).  This  velocity  dau  shows 
scatter  of  between  a  factor  of  3  snd  5  for  the  expiosions  represented. 

Moving  ouuide  ihe  spall  zone,  a  typical  FSE  array  for  characterizing  a 
nuclear  explosion  is  given  In  Figure  3x  The  dimension  of  the  array  is  on 
the  order  of  10  km.  Careful  attention  is  paid  to  coverage  in  both  range 
and  azimuth.  Most  giges  are  accelerometers  although  velwiiy  transducers 
are  sometims  used  at  more  distant  ranges.  Peak  vertical  accelerations 
from  a  tingle  explosion  observed  at  a  nunfaer  of  azimuths  and  ranges  are 
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Fig.  3.  (a)  A  typical  free  surface  instrainent  atny  for  characteroaiion  of 
the  explosion  souice  fiinction.  The  xceleromeiers  lie  all  outside  the  spall 
zone  (FSE  region)  and  designed  to  give  good  azimuthal  coverage  of  the 
source,  (b)  Peak  vertical  accelerations  are  plotted  against  free  surfxe 
range  (unsealed)  for  a  single  nuclear  explosion.  The  large  scatter  (6-8)  in 
these  observations  with  range  is  evident. 

given  in  Figure  3b.  Accelerations  span  over  an  order  of  magnitude  in 
amplitude  and  slightly  less  than  an  order  of  magnitude  in  range.  A 
distina  amplitude  decay  with  range  is  observed  althmgh  it  is  obscured  by 
a  factor  of  6-8  scatter  in  peak  accelerations. 

The  single  explosion  data  can  be  compared  to  scaled  f-nit  xeelerations 
(FSE)  from  a  number  of  Pahute  Mesa  explosions,  one  of  which  was 
detonated  below  the  water  table  (Figure  4).  The  same  characiertsucs 
observed  in  the  single  explosion  dau  are  replicated  in  the  multiple 
explosion  data  set.  There  is  an  approximate  order  of  magnitude  scatter  in 
peak  acceleration  data  with  the  shot  from  below  the  water  table  exhibiting 


SCALED  RANGE  (m) 

Fig.  4.  Peak  accelerations  fiom  a  number  of  Pahute  Mesa  explosions  are 
sealed  by  the  cube  root  of  yield  and  plotted  against  scaled  free  surface 
range.  All  data  is  from  outside  toe  spall  zone  (FSE  region).  A  large 
amount  of  scaiier  (6-8)  is  observed  in  this  data.  There  is  some  indication 
of  a  coupling  effen  as  the  data  from  an  explosion  detonated  below  the 
water  table  (solid  squares)  picu  at  the  upper  bounds  of  the  tixm 

the  highest  accelerations.  The  physical  mechanism  for  the  scaoer  in  these 
data  could  either  be  different  wave  propagation  effects  (depth  of  burial, 
proaimity  to  interfaces,  different  source-receiver  structures)  or  source 
coupling  effects.  Large  amplitudes  for  the  source  below  the  water  table 
support  a  coupling  effect.  On  the  ether  hand,  large  scatter  in  the  single 
explosion  dau  (Figure  3b)  suggests  propagation  path  effects. 

Scatter  in  observational  data  is  frequency  dependent.  The 
accelerograms  used  for  the  peak  amplitude  snidy  (Fgure  3b)  are  Fourier 
transformed  and  long  period  sprcual  estimates  are  made  and  converted  to 
moment: 

Ml  ■  (1) 

where  p  is  density,  R  is  s'lant  range,  a  is  compressionaj  velocity,  and  Clg 
is  the  long  period  dupbeemem  spectral  leveL  Since  the  dau  is  from  the 
free  surface,  a  simple  'actor  of  two  is  introduced  into  the  Oo  estimate. 
Froblems  with  this  simple  interpretation  of  the  free  surfxe  data  include 
the  inability  to  separate  bodv  and  surfxe  waves  in  these  dau  sets. 
Isotropic  momenu  interpreted  Crcm  vertical  (2).  radial  (R),  and  transverse 
(T)  dau  are  given  to  Figure  3  (srictiy  speaking  only  Z  and  R  should  be 
used).  The  scaiirr  in  momenu  from  this  single  explosion  is  a  fxtor  of 
2-3,  redxcd  from  variations  in  peak  accelenticn.  The  fxt  that  Z.  R.  and 
T  momenu  all  ibilow  the  same  spatial  pattern  argues  that  the  variation 
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RANGE  (km) 


Fig.  3.  Scumie  motnems  from  ihe  Pahuie  Mea  nuclear  explosion 
(Fi{.  3)  were  deiemined  from  observaiions  at  ranges  between  2.1  to 
SJktn.  Moments  were  calculated  from  vertical  (Z),  rvlial  (K)  and 
Oansvene  (T)  observations.  The  scatter  in  the  moment  estimates  is  a 
(actor  of  2-3,  reduced  from  the  scatter  in  peak  accelerations  given  in 
Figs.  3b  and  4. 


may  be  a  propagation  effect.  Similar  variations  are  found  for  moments 
estimated  from  near-source  earthquake  data  (Archuleta  et  al.,  1982; 
Fletcher  et  al.,  1984J.  Following  the  methodology  of  Archuleta,  average 
log  moments  and  their  standard  deviations  are  computed  according  to; 


log  (Me) 


(I/ns)  J  log  (Moi) 
L  i-l 


(2) 


std  dev  (log  (Mo)) 


ns 


I 

i-l 


log  Moi -log  (Mo))^ 


1/2 


0) 


"nie  multiplicative  error  is  defined  according  to: 

EMo  »  antilog  fird  dev  (log  (Mq))]  (4) 

The  average  moment  for  the  vertical  data  in  Figure  6  is  3.8  x  1023  dyne- 
cm  with  a  multiplicative  error  of  1.49. 


Changes  in  frequency  content  with  propagation  distance  for  data 
summanzed  in  Figure  3  are  characienaed  by  picking  comer  frequencies 
from  raw  accelerograms.  Accelerauon  spectra  are  interpreted  in  terms  ot  a 
long  period  nse  that  is  used  in  the  moment  estimate,  an  intermediate 
frequency  band  in  which  aceeleiauon  spectra  are  fiat,  followed  by  high- 
frequency  decay  (Figure  6).  A  first  comer  frequency  is  chosen  at  the  point 
of  iransiuon  from  long  period  nse  u>  the  constant  level  at  intermediate 
frequencies.  A  second  comer  frequency  is  identified  at  the  high-frequency 
point  where  the  accelerauon  spectra  begin  lo  decay.  These  two  comer 
fiequencies  are  plotted  against  source-receiver  olTsei  for  a  single  explosion 
in  Figure  6.  The  first  comer  frequency  changes  little  with  range  and  is 
esumated  to  be  1  Hz.  This  spectral  charactensuc  is  taken  to  be  the  source 
comer  frequency  which  is  rebted  to  the  source  ehsuc  radius.  The  second 
comer  shows  a  systemauc  decrease  with  range  indieauve  of  a  propagation 
path  effects.  Simple  frequency  independent  auenuauon  openajis; 

e-«frA3o 


suggest  that  a  Q  between  20-30  replicates  this  decay  with  range. 

These  frequency  domain  measures  show  less  variation  than  the 
amplitude  information  for  this  particular  explosion.  The  fiat  acceleration 
spectrum  between  the  two  comer  frequencies  (Figure  5)  argues  for  a  f~2 
high-frequency  source  model  for  this  explosion.  As  source-receiver  offset 
increases,  this  iniermedbie  high-frequency  decay  becomes  obscured  by 
aoenuauon  effects. 

Yucca  Flats  (FSS  <i  FSE) 

Yucca  Flats  is  a  valley  filled  with  alluvium  at  shallow  depths 
underlain  by  dry  and  wet  tuffs.  The  basement  consists  of  Paleozoic  rocks. 
A  set  of  FSS  dau  from  the  COALORA  experiment  b  reproduced  in 
Figure  7.  The  depth  and  range  of  the  spall  phenomena  is  constrained  with 
this  data.  Spall  eaiends  in  this  case  to  a  free  surface  range  of  1  DOB 
while  the  de^  of  the  spall  zone  is  esumated  to  be  1,2  DOB.  SpaJ  dwell 
time  increases  near  GZ  as  noted  for  Pahute  Mesa  data.  Spall  zone 
iceeietauons  ate  used  to  estimaie  escape  velociiiea  (after  integration)  of 
the  spalled  mass.  Spall  volume  is  delineaied  by  the  spatial  distribuuon  of 
accelerograms.  These  dau  combined  with  the  equivaiem  body  force  spall 
model  of  Day  et  al.  (1933)  are  used  to  estimaie  an  equivalent  body  force 
lime  history  for  the  secondary  source.  The  model  predicts  a  peak  force  of 
l.l  X  10^^  dynea  and  a  lime  duration  of  0.6  s. 

Peak  acceleration  dau  from  this  explosion  which  mark  the  transition 
from  FSS  to  FSE  are  pkxied  againat  free  surface  range  in  Fgure  8.  As 
found  for  Pihuie  Mesa  data,  explosion  geometry  coniiois  the  gradual 
decay  of  the  dau  close  to  GZ  followed  by  more  rapid  decays.  Radial  and 
vertical  peak  motions  merge.  Scatter  in  the  peak  motions  b  between  a 
factor  of  3-5  for  the  Y ucca  Flats  single  explosion  data.  Thb  variation  b 
less  than  the  single  shot  Pahuu  Mesa  acceleniion  scaoer. 

Differences  in  wave  shapes  and  spectra  from  COALORA  are 
documented  with  observational  displacements  at  the  549  m  range 
dbpbyed  in  Fgure  9.  Peak  radial  and  verucal  displacements  vary  by  30^!i 
about  the  mean.  Radial  and  verucal  wave  shapes  are  very  similar  at  the 
three  azimuths.  These  characteristics  argue  for  cylindrical  or  spherical 
symmetry  in  the  source  funciioa  Transverse  components  in  comparison 
are  much  shorter  in  durauon,  delayed  in  lime,  and  exhibit  significant 
fluctuations  with  azimuth  including  changes  in  first  motion. 
Displacement  spectra  are  eatimaied  from  these  records  (Figure  9). 
Envelope  funcuons  are  fit  UJ  each  spectrum  which  include  a  long  period 
level,  corner  frequency,  and  high-frequency  decay  (Table  2).  These  rbu 
show  little  variation  in  radial  and  vetticai  long  period  levels  and  a  factor  of 
two  increase  in  corner  frequency  for  the  transverse  component 
xcompanied  by  reduced  long  peried  levels.  Source  spectral  inietpreution 
of  the  above  dcia  supplemented  by  additional  near-source  gauges  gives  an 
isotropic  moment  (R,  Z)  of  1.95  x  lO^*  dynes-cm  (multiplicauve  error  of 


14 


/ 


FREE-FIELD  AND  FREE  SURFACE  GROUND  MOTIONS 

10- 


8 


1 

□ 


■  Z  CORNER  FREQ  (1) 
□  Z  CORNER  FREQ  (2) 

KOL.  d 
t  froncu 


61 


. 

b: 
o 
o 


21 


S 


f 


— r— 

4 


8 


10 


RANGE  (km) 


Fig.  (S.  Comer  frequency  interpreuiion  from  acceleration  spectra  of  the  same  data  used  in  Fig.  5.  As  shown  in  the 
inset,  two  comer  Gtequencies  were  measured.  The  first  represents  the  transition  from  the  rise  in  acceleration  spectra  at  tow 
frequoKy  to  the  natiening  at  intermediate  frequencies.  The  secord  comer  frequency  marks  the  transidon  from  the  flat, 
intermediate  fixqueney  region  to  the  decay  at  high  frequencies.  The  ftrs'.  comer  frequency  (1  Hz)  which  is  attnbuted  to  the 
sooice  is  insensitive  to  range  while  the  second  comer  decreases  with  range  from  9  to  3  Hz.  The  flat  spectra  between  the 
two  comets  supporu  a  source  model  The  second  comer  frequency  is  aohbuted  to  attenuation. 


1J6)  and  a  deviatoric  moment  (T)  of  3.97  x  dynes>cm 

(multiplicative  enor  of  t.A4).  The  accompanying  comer  frequency 
estimates  are  1.82  for  the  isotropic  spectra  (1.12  multiplicative  error)  and 
3.23  for  the  deviatoric  spectra  (1.13  multiplicative  error).  As  found  for 
Pahute  Mesa  spectral  data,  comer  frequencies  show  less  variation  than 
moments. 

The  isotropic  comer  frequency  lies  between  a  Mueller-Murphy  model 
prediction  of  2.3  Hz  for  1  lu  and  1.7  Hz  for  10  kt  which  is  in  agreement 


with  the  announced  yield  of  less  than  20  kt.  The  high  frequencies  decay 
IS  f'^,  like  the  Pahute  Mesa  data  and  the  Mueller-Murphy  source  model. 

A  simple  Brune’s  model  is  used  for  comparative  purposes  to  Interpret 
the  transverse  spectra  in  terms  of  stress  drop  and  source  dimension.  These 
estimates  are  applicable  only  for  earthquake  sources  and  may  not  be 
saicily  applied  to  tectonically  driven  motions  triggered  by  the  explosion. 
The  parameters  predicted  by  this  model  are  given  for  comparison.  The 
mean  source  radius  is  124  m  with  a  stress  drop  of  89  bars  and  an  average 
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Verticai  crass-iecuon  displaying  accderognnu  from  downhole  and  surface  pugei  within  a  range  equal  to  two  depths  of 
burial  which  constrains  the  spall  zone.  Solid  circles  represent  vertical  accelerometer  records  that  show  characteristic  -Ig  dwell 
indicative  of  i-palL  Solid  and  long  dashed  lines  are  two  dimensional  bounds  on  the  spall  acne. 


Free  Surface  Range  (krn) 


Fig.  8.  Peak  vertical  (solid  squares)  and  radial  (open  squares)  accelenuoo 
picued  against  free  surface  range  (unsealed)  for  the  COALORA  eaplosion 
at  Yucca  Flats.  This  data  like  that  for  Pahute  Mesa  in  Fg.  2  spans  the 
transition  from  the  free  surface  spaJI  (FSS)  to  the  free  surface  elastic 
(FSE)  regions. 


dispiacemeatof  33  cnL  These  displactmeais  calculaied  horn  the  radiated 
wavefield  are  similar  ia  size  to  dispUcemenu  observed  on  fiiuts  and 
bedding  planes  found  upon  re-enusing  uautels  surrounding  nuclear 
explosions  Rainier  Mesa  (Kennedy,  198A].  The  equivalent  elastic  radii 
for  1  and  10  kt  explosions  n  predioed  to  be  between  133  and  202  m. 
The  deviaioric  source  radius  falls  near  the  lower  bound  of  the  equivalent 
elastic  axzce  radius. 

As  seismic  waves  propagaia  within  Yucca  Flat  valley,  differences 
between  radial,  vertical,  and  nansvene  waveforms  and  spectra  disappear. 
Displacement  records  and  spectra  at  the  5.16  km  distance  ate  given  in 
Figure  10.  The  duration  of  ground  motion  has  grown  from  1.2  s  at  549 
m  to  nearly  20  $  at  5.16  km.  The  tnrtsvene  motion  it  the  largest  of  the 
three  components  and  comparable  ia  frequency  to  the  other  two 
components.  Spectral  differences  found  a  the  closest  ranges  and  attributed 
®  deviatoric/tsotropic  source  proeexies  have  disappeared  a  this  distance. 

Rouurr  Mtja  (FFW  i  FS£) 

Data  from  the  free*(ield  can  be  used  ®  constrain  coupling  of  explosive 
energy  into  the  linear  regime.  It  is  need  in  the  calculation  of  the 
representauve  seijrtuc  soiaee  function  known  u  Use  reduced  displacement 
potential  (Murphy  and  Bennett.  1979).  As  noted  by  Murphy  (1989), 
exisung  free-fieJd  data  is  often  complicated  by  complex  geological  models 
between  explosion  atd  receiver.  Much  of  the  existing  free-field  ground 
motion  data  bass  is  ilso  issoeiaied  with  underground  structures. 
Insirumenu  are  pnmanly  fielded  to  document  survivability  levels  for 
these  structures  and  not  free.field  motiont  necessary  for  seismic  source 
chaiacatzatica 

An  experiment  was  designed  »  quantify  the  seismic  souree  of  a 
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Fig.  9.  Vertical  (2)  and  nansvene  (T)  displacemems  spectra  along  with  ihe  displacement  waveforms  from  !he  COALCRA 
explosion.  The  three  stadons  (GMT,  Cl  18.  CM9)  are  all  at  a  free  surface  range  of  549  m.  Peak  displacements  for  each  wavefcrm 
ate  given  in  cm.  A  simple  spectral  model  consisting  of  a  constant  long  period  level,  high  frequency  decay  (r").  and  a  cort'er 
bequency  is  fit  to  the  specaal  data. 

Rainier  Mesa  explosion  using  free-lield  (FFW)  and  !ree  surface  (FSS  £  (19751  for  wet  tulT.  R~^^^  (aeceieration).  Waveforms  show  ;tn  azimuthal 

FSE)  puger  (Figure  lla).  Free-neld  measurements  were  made  with  effect  in  wave  shape  with  two  acceieromeiets  to  the  NE  having  longer 

three-component  accelerometers  locaierl  to  minimize  geological  effects  as  durations  than  gauges  to  the  NW.  Ouradon  of  motion  also  increases  with 

well  as  document  the  transiuon  of  the  modon  field  into  the  elasdc  range  as  noted  by  Murphy  [19391.  Another  measure  of  propagauon 

regime.  Free-surface  puges  were  fielded  so  ihat  comparisons  between  the  complexity  is  the  rado  of  radial  us  vertical  peak  acceleration  (free-neld) 

two  dau  sets  ..mid  be  made.  The  radial,  free-field  dau  is  reproduced  in  which  is  5.1  for  the  closest  gauge  (193  m)  and  decreases  to  l.S  for  the 

Figure  I  le  with  absolute  rmpliuides  compared  to  the  free-surface  dau  in  farthest  gaup  (38?  m). 

Figure  lib.  Amplitude  dau  indicates  a  smoothly  decaying  xcelerauon  Companson  of  free-neld.  ndial.  peak  acceletadon  (open  square)  and 
and  velocity  Held  which  follows  model  piedicdons  by  Ferret  and  Bass  free-surface,  peak  xcelerauon  (solid  square)  indicates  increased  scauer  for 
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Table  2.  COALORA  SptotI  Inierpreouon 


Station 

Range 

(m) 

DC 

(em-$) 

<c 

(Ha) 

Slope 

CM7Z 

612 

0.30 

1.8 

2 

CM7R 

612 

0J6 

2.1 

2* 

GMTT 

612 

0.09 

4.0 

4 

CMSZ 

614 

0.30 

1.6 

2 

CM8R 

614 

0J6 

1.8 

X 

CMSr 

614 

0.10 

3.3 

2* 

GM9Z 

612 

0.34 

1.8 

r 

GM9T 

612 

0.30 

3.3 

3* 

Z2.9S 


Fit  10-  Venieai  (2).  ndial  (R)  and  nnsvene  (T)  ditplacementa  and 
displacefnem  ipectn  at  iJie  j.  16  km  lource-recaver  range  for  COALORA 
are  displayed.  The  Z  and  R  spectra  are  aiditnniy  scaled  by  factonoriOO 
and  10  so  that  all  three  can  be  compared.  Where  the  Z  and  T  spectra 
showed  signiricantly  differeni  spectral  levels  (Qz  >  Or)  and  corner 
Gtequency  (fCT  »  fCZ)  “  5^9  m  (Fig.  9),  at  M6  km  both  spectral  levels 
and  comer  fraptenctea  for  ail  contpooents  are  idenucaL 


free-stufacedaia.  The  smgle  awree  scatter  in  free-surface  dau  is  similar  to 
that  from  Pahuie  Mesa.  Free-neld  scatter  ts  much  reduced  which  might  be 
a  reflecuon  of  the  fact  that  receivers  have  been  moved  away  from  the 
weathered  zone.  At  approximately  SOO-900  m.  free -field  and  free  .surface 
data  merge  with  no  factor  of  two  ampUncaiton  at  the  free-surface.  This 
observation  is  only  an  apparent  discrepancy  since  free  surface  data 
displayed  in  Figure  I  lb  are  gauges  from  outside  the  spall  zone  (mangles 
in  Figure  I3ai  and  as  such  represent  oolique  raiher  than  normal  free 
surface  incidcm  angles. 

Wavefieid  Modeling  and  Inversion 


FSS  to  FSE  tranituon 

The  uansiuon  from  FSS  to  FSE  is  idemined  by  a  region  of  gradual 
peak  amplitude  decay  to  a  range  near  100  m/kt  followed  by  more  rapid 
decay  at  greaier  ranges  for  both  Pahuic  Mesa  (Figure  2)  and  Yucca  Flats 
(Figure  8).  These  data  argue  that  explosion  geometry  has  a  dominant 
effect  on  ground  motions.  An  elastic  numerical  modeling  exercise  is 
undertaken  for  this  nnsition  region  at  NTS  in  an  attempt  to  quantify 
wave  propagation  effects  and  develop  an  understanding  of  free-surface 
motion  decay  rates.  The  velocity  model,  HOLE,  used  in  this  exercise  is 
given  in  Figure  I2a  ami  is  developed  from  emplacement  hole  data  at 
Pahuie  Mesa  (Table  3).  A  second  avenge  velocity  model  developed  by 
Leonard  and  Jonnson  [1987],  U,  is  given  in  the  figure  for  comparison. 
The  explosion  source  is  placed  at  a  depth  of  616  m  just  above  a  poorly 
welded  luff  that  shows  reduced  P  wave  velocity. 

Full  wave  synthetics  using  the  modiriad  renectivity  methodology 
[Muller,  198S]  are  computed.  Vertical  and  radial  velocity  waveforms 
between  0.1  and  3.0  km  are  reproduced  in  Figure  12b.  Waveforms  are 
multiplied  by  r^  to  balance  ampliiudes  for  viewmg.  Figure  13  displays 
the  syntheuc  velocities  plotted  as  a  futKuon  of  scaled  range  where 
the  explosion  yield  is  ufcen  to  be  ISO  kt  for  the  source  depth  of  616  m. 
Vertical  waveforms  show  liiUe  decay  to  approximately  lOO  m  (scaled 
range)  with  inoeased  decay  beyoitd  this  bee  surfxe  range.  Least  squares 
fit  to  Ksspliuide  data  resuit  in  a  spatial  decay  of  r~9-l^  u>  a  scaled  range  of 
100  mAi*^  while  the  decay  rate  at  greater  ranges  is  These 

theoretical  results  compsre  to  Pahuie  Mesa  observauons  (Figure  2)  that 
decay  as  r~9-i3  and  r*i-22  and  indicaia  that  elastic  wave  propagation 
effects  which  irwlude  auenuatioo  match  the  obser/ation  tian.*iuon.  Data 
and  models  imply  that  spall  momentum  comnbuuons  from  the  region 
inside  of  100  m/kt’^  range  may  be  dominant  due  u>  increased  spatial 
decay  at  greaier  ranges. 

Spall  Ballhae  Medal 

Vidocities  and  displacements  are  derived  tram  acceletation  observations 
within  FSS  region  for  the  Pahuie  Mesa  data.  These  measurements  as 
well  as  spall  dwell  time  are  inienelaied  for  a  gravity  controlled  process.  A 
gravity  driven  or  ballisuc  model  of  (he  spall  process  predicts  the 
following  relation  will  hold  between  peak  velocity  at  failure  and  peak 
displacemenc 


D- 


vJ_ 

2( 


(6) 


Peak  displacement  is  plotied  againsx  velocity  at  failure  (spall  initiation)  in 
Fgure  14.  The  leUuonship  predicied  by  equation  6  is  tksgnaied  by  open 
squares.  Observational  dau  follow  the  slope  of  predictions  but 
displacements  are  aa  much  as  a  factor  of  two  bigger  than  predictions 
throughout  the  dau  set 

Although  large  displacements  may  be  a  result  of  bias  introduced  in  the 
integration  procedure  m  deriving  these  values,  the  consistency  of  large 
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(b) 


SLAhn*  RANGE  (m> 


«.M  I.IM  U9  Ml  (lM  I.U  (.I* 


Ul  t.11  US  Ut  (.<•  U*  L7S 


IW  (s) 

Fig.  11.  (a)  Plan  view  of  ihe  inamuneei  anay  for  a  combined  free-rield  (solid  circles)  and  free  suface  (open  circles  and  ail 
triangles)  ground  motion  espenmeni  at  Rainier  Mesa,  NTS.  (b)  Compahson  of  the  peak  ndial  accelerations  in  the  free-Held 
lopen  squares)  and  at  the  free  surface  (solid  squares).  A  marked  increase  in  data  scatter  is  noted  (or  the  free  surfKe  data. 

(e)  Normalized  free-fieid  radial  velocity  waveforms  from  the  experimenL 

displacements  for  a  number  of  shots  at  i  variety  of  motion  levels  argues  scalier  in  the  ne3r>tource  region  is  a  result  of  lateral  variations  in  the 

against  this  interpretation.  Alternatively,  departure  from  a  purely  geological  structure  at  NTS.  There  are  strong  vertical  variations  in 

gravitational  spall  model  may  be  explained  by  viewing  spall  as  a  velocity  in  addition  to  lateral  variations.  Figure  12a  displays  the  velocity 

continuous  rather  than  a  simple,  discontinuous  prxxess.  Even  though  log  developed  from  an  emplacement  hole  at  Pahute  Mesa  (HOLE)  and 

peak  velocity  indicates  failure,  spalled  material  may  still  receive  long  contrasts  it  with  an  average  model  developed  by  Leonard  and  Johnson 

period  input  iiom  material  below  thus  boosting  displacements  over  those  (1987)  from  inversion  of  near>source  travel-time  dau  (LI).  The  site 

predicted  by  a  simple  ballistic  model.  More  dau  supplemented  with  specific  velocity  model  shows  strong  low  velocity  zones  (LVZ) 

numerical  modeling  is  needed  to  investigate  this  point.  represenutive  of  poorly  welded  tuff  uniu.  Differences  in  synthetics 

developed  from  the  average  and  site  specine  velocity  models  are 
Vtladij  Model  Effeca/Data  SeaiUf  investipied  to  determine  how  much  detail  must  be  included  in  the 

velocity  modeL  A  set  of  numerical  trials  will  be  dlmiww  using  the  site 
Sealed  acceleration  and  velocity  dau  from  Pahute  Mesa  (Figure  S)  speciAc  velocity  model  with  sources  just  above,  in,  and  below  a  LVZ 

show  considerable  scatter,  A  number  of  authors  have  argued  that  such  such  as  that  which  extsu  below  616  m  (HOLE)  in  Figure  12a.  This 
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(a) 

PAHUTE  MESA  P  VELOCITIES 


DEPTH  ten) 


OFFSET  (0,l-3.0tim) 


OFFSET  (0.1-3.01I*) 


Fig.  12.  (i)  Th«  P  wtv«  ydodty  aodel  inietpnted  from  empixemem  Me  data  M  PaAuie  Meaa  (jolid  iqoaies. 
Table  3)  cofltnuted  against  an  avoige  Pahute  Mesa  veiociiy  model  (open  circles)  developed  by  Leonatd  and 
Johnson  (1987).  (b)  Verncal  and  radial  syniheue  veiociiy  recotds  de^-eioped  (torn  the  site  specific  velocity  model 
using  ilK*  extended  leneetiviry  modeling  lechmque.  The  symheucs  caimlaifd  for  ihe  Cree  surface  langcs  of  0.1  lo 
33)  Im  have  all  been  scaled  by  free  surface  range  (r*)  lo  balance  anplioides  for  display. 


numerical  experiment  will  help  determine  the  efTect  of  oeierul  prepenies 
around  the  explosion  on  near'ssuree  wavefields. 

A  number  of  researchers  p/idale  and  Helmberger.  1987;  Johnson, 
1988;  Slump  and  Johnson.  1984]  have  noted  the  impottance  of  both  body 
and  surface  waves  in  near-source  wivefonns.  Although  distances  over 
which  observations  are  made  in  this  study  are  tdativefy  short  (telO  bn), 
the  shallow  source  depths  (2(X)-I0(X}  m)  lead  to  these  effects.  The 
modified  reflectivity  method  (Muller,  1983]  is  employed  in  seismogram 
synthesis  to  include  both  body  and  surface  wave  contributions.  Jovreei 


in  these  calculations  have  a  constant  moment  with  a  source  lime  function 
appropriaie  for  a  130-b  explosion.  A  simple  separation  of  dominant 
body  and  surface  wave  contnbutiona  ia  found  by  low-past  fillenng  the 
data  at  I  Hi  n  emphasize  siafaoe  waves  and  high-paaa  fillenng  ihe  data  at 
1  Hx  10  focus  on  body  wavs.  Peak  velociiiea  ai«  ptdred  from  lyntheucj 
for  compansen  belweea  the  varioua  models. 

The  Site  specific  velocity  model  (HOLE)  is  flist  investigated  with 
sources  above  ihe  LVZ  (616  m),  in  the  LVZ  (665  m),  and  below  the  LVZ 
(750  m)  (Figure  12a).  The  body  wave  dau  (Figure  15)  shows  strong 
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Fif.  13.  Peak  venical  velodciea  from  Fig.  12b  plooed  against  scaled  free 
surface  range.  Source  depth  for  the  synthetics  in  Fig.  I2b  was  appropriate 
for  an  explosion  of  ISO  xi  and  so  the  free  surface  ranges  were  scal^  by 
this  yield.  As  done  for  the  Pahute  Mesa  data  displayed  in  Fig.  2.  power 
law  decay  models  were  fit  to  the  synthetic  amplitudes  for  ranges  less  than 
and  greater  than  100  m/kt^. 


Velocity  (m/s) 

Fig.  14.  Peak  displacements  and  velocities  from  a  number  of  spall  zone 
observations  of  different  explosions  ate  plotted  against  one  another.  Fora 
gravitationally  controlled  process  peak  displacement  is  rehted  u>  velocity 
according  to  equation  6.  This  relation  is  displayed  as  open  squares  in  the 
figure.  On  average  for  a  given  peak  velocity  the  observed  peak 
dispbrements  are  greater  by  a  factor  of  two  than  that  predicted  by  the 
gnvitadonal  model 


Fig.  15.  Peak  vertical  body  wave  velocides  (f  >  1  Hz)  in  the  FSE  region 
for  extended  tefiectiviiy  synthetic  seismograms  utilizing  the  HOLE 
velocity  model  (Fig.  12a)  for  scuice  depths  of  616  m  (shove  a  low 
velocity  zone.  LVE).  663  m  (in  LVZ)  and  730  m  (below  LVZ).  Also 
included  are  peak  venical  body  wave  velocities  (f  >  1  Hz.  U  616  m) 
generated  from  the  average  Pahute  Mesa  model  of  Leonard  and  Johnson 
(1987). 

amplification  (factor  of  2)  for  body  waves  directly  above  a  source  when 
the  explosioa  is  above  the  LVZ.  Fbr  a  ISO-kt  shot  the  spall  zone  extends 
to  an  approximate  lateral  range  equal  to  2  DOB  which  includes  this  region 
of  enhanced  amplification.  For  the  souice  in  or  below'the  LVZ  little 
variaoon  is  observed.  Differences  between  the  three  models  disappear  as 
source-receiver  olbet  increases  and  downgoing  energy  becomes  more 
importanl  Surface  wave  synthetics  below  1  Hz  show  small  differences 
between  the  three  models. 

Comparison  of  peak  body  wave  amplitudes  for  a  616-m-Jeep  explosion 
in  the  average  Pahute  velocity  model  of  Leonard  and  Johnson  [1988]  are 
also  given  in  Figure  15.  Differences  between  the  two  models  are  less 
disi  those  observed  for  sources  above  and  below  a  LVZ. 

Momtia  Timor  Inveniam 

Inverse  modeling  of  observed  esplosion  waveforms  offers  another 
method  of  making  source  estimates.  Moment  tensor  inversions  of 
COALORA  data  (Yucca  Flats)  are  completed  in  an  attempt  to  utilize  the 
complete  FSE  data  in  a  manner  which  gives  stable  source  estimates  in 
spile  of  demonsuuied  ampliuide  variability.  The  following  representation 
is  used  in  these  inversions: 

Un(0-G„ij(0Mij(0  (7) 

where  UnfO  nre  the  frequency  domain  reptesenotion  of  the  observations. 
Gnij(0  are  the  Green’s  functions  or  wave  propagation  effects  accounted 
for  with  extended  reficctivity  solutions,  and  Mij(0  are  the  componen's  of 
the  moment  Umsor. 

Many  different  stations  are  utilized  simultaneously  in  these  studies  in 
an  attempt  to  minimize  the  effect  of  data  scatter.  Both  absolute 
amplitudes  and  wave  shapes  are  modeled  with  fits  completed  in  the 
frequerKy  domain.  Characteristic  comparisons  between  observed  and 
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modeled  pound  mouons  pve  eotrebuon  coelTicienu  liui  range  from  0.7 
to  0.9  in  the  ume  donum.  Tlie  moment  tenior  that  wu  determined  in 
this  invetiion  is  given  in  Figum  16.  The  source  is  dominated  by  me 
isotropic  component  The  initial  pulse  is  symmemc  and  is  followed  by  a 
secondary,  long  penod  conmibuuon,  largest  on  the  M33  component,  me 
vetucai  dipole.  The  size  of  mis  secondary  force  moment  is  in  general 
sgreement  wim  me  equivalent  elastic  force  model  of  spall  developed  from 
data  wimin  me  spall  zone  (Figure  7).  Off.di3gonal  moment  tensor 
elements  are  a  factor  of  3-10  smaller  man  diagonal  elements  supporung 
the  small  scalar  deviatoric  source  estimates  discussed  earlier.  Peak 
isotropic  moment  is  7.9  x  10^^  which  is  nearly  a  facur  of  3  smaller  man 
t.  /iivf  moments  estimated  from  spectra  at  549  m.  The  scalar  moment 
interpretation  of  the  549  m  spectra  assumes  a  simple  propagauon  pam 
effect  wim  no  way  to  return  downgoing  energy  to  me  free  surface  while 
inversions  uulize  Green's  functrons  which  return  downgomg  energy  to  me 
surface  and  lead  to  reduced  moment  estrmaics. 

Consuained  explosion  source  estimates  are  made  in  an  attempt  to 
investigate  bias  from  single  sution  measurements.  The  source  is 
assumed  isotropic  and  a  single,  ihree<omponent  set  of  observations  are 
inverted  for  the  lime  history,  T(0,  and  strengm.  ,Mu,  of  me  source: 

Un(0-G„kjt(0M|ckT(0  (8) 

Each  instnuneni  location  around  the  explosion  yields  a  source  estimate 
much  like  standard  scalar  estimate  aimough  equation  g  allows  for 
inclusion  of  mote  complex  propagation  pam  effects.  A  moment  for  each 
station  was  determined  wim  mis  methodology  for  COAX.ORA.  The 
average  of  all  these  single  sution  momenu  is  6  j  x  10^  dyne-cm,  close 
to  the  value  of  full  moment  tensor  inversions  and  signiTtcantly  smaller 
than  simple  scalar  moment  estimates.  Using  equations  3  and  4  10 
quantify  scalier  in  mese  estimates,  me  multiplicative  error  is  1.6.  close 
to  that  of  the  scalar  motrenu  determined  from  displacement  specn. 

Conclusiona 

A  number  of  near-source  data  sets  from  explosions  have  been  used  to 
constrain  the  explosion  source  function.  The  study  has  included  fice-fleid 
data  collected  at  shot  level,  ftee  surface  data  from  within  the  spall  zone, 
and  free  sur&ce  data  from  outside  me  spall  zone.  Data  from  explosions  at 
Rainier  Mesa.  Pahute  Mesa  and  Yucca  Flats  NTS  were  included. 

The  Rainier  Mesa  data  illusiraie  me  increased  variabiliiy  of  free  surface 
observations  compared  to  mose  made  in  the  free-Held  (at  shot  depm). 
Utile  free-Field  dau  (FFW)  scatter  was  found  about  the  spodal  decay  rate 
predicted  by  Perret  and  Bass  (19751  for  wet  uiff.  The  free  surface  data 
(FSE)  .'rom  me  same  explosion  had  a  factor  of  6-g  scatter.  This 
observational  result  illtisuaies  me  impact  of  weamered.  near-surface  layers 
on  ground  motions  and  indicates  difficulties  that  might  be  encountered  by 
making  source  esbmaies  from  limited  amounts  of  bee  surface  data. 

Observational  dau  from  wimin  the  spall  zone  (F5S)  transitioning  to 
the  elastic  region  (FSE)  was  analyzed  for  z  number  of  Pahute  Mesa 
explosions.  These  dau  illustrate  me  strong  impact  source  geomeuy  fas 
on  peak  motions  and,  possibly,  spall  as  a  secondary  seismic  source. 
From  gnnind  zero,  directly  above  me  explosion,  to  a  lateral  free  surface 
range  equal  to  1  DOB  peak  velocity  dau  shows  little  deeav, 

Beyond  this  range  me  dau  decays  much  more  rapidly,  f Simple 
elastic,  reneeiivity  calculations  for  a  layered  Pahute  Mesa  structure  wim 
appropriate  Q  (Tabic  3)  replicate  mis  decay  pauenv  The  dau  and  these 
numerical  results  suggest  that  momentum  wimin  the  spall  zone,  which  is 
available  for  a  secondary  seismic  source,  may  be  most  affected  by  spall 
processes  in  me  lateral  range  extending  to  I  DOB  for  sandard  conuinmem 
depms. 

Observational  peak  acceleration  and  velocity  dau  from  wimin  the  FSE 
region  for  Pahute  Mesa  show  a  scatter  as  large  as  a  fxtor  of  6-g  for  bom 
a  single  explosion  and  a  number  of  scaled  explosions.  These  observations 
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Ftg.  16.  Moment  rate  (M)  and  moment  (M)  lenson  from  an  inversion  of 
me  COALORA  observa'joral  data.  The  1  and  2  direcuons  are  in  me 
honzomal  plane  while  me  3  dircetion  is  vcnical. 


Table  3.  Velocity  Model  for  Hole^ 


Depm 

fkmi 

P  Vel 
(Vm/$> 

S  Vel 
fkmj'sf 

Densny 

femfccf 

Qp 

QS 

0.000 

1.00 

048 

1.70 

50.00 

2ZJ2 

0.100 

240 

147 

2.10 

50.00 

r42 

0450 

3.00 

1.73 

240 

30.00 

ZL22 

0450 

1.00 

048 

1.70 

30.00 

2242 

0.400 

1.00 

048 

1.70 

50.00 

2242 

0.450 

1.00 

048 

L70 

50.00 

2242 

0.300 

340 

2.02 

2.40 

30.00 

2242 

0.616*’ 

340 

i02 

140 

50.00 

2242 

0.630 

340 

2.02 

140 

50.00 

2242 

0,630 

240 

1.44 

110 

50.00 

2242 

0.700 

240 

1.44 

110 

50.00 

2242 

0.700 

342 

1.99 

140 

100.0 

44.44 

1.000 

3.70 

2.03 

140 

100.0 

44.44 

2.000 

4.37 

249 

140 

1000 

44.44 

4  000 

5.25 

3.03 

140 
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44.44 

4  Linear  interpolation  used  between  depth  points  except  in  case  uf 
repeated  depths  which  signifies  a  discondnuiiy  in  velocity. 

Souree  depth. 


support  the  conclusion  fiom  the  Rainier  Mesa  dau  that  the  weamered 
layer  at  me  earths  surface  can  have  significani  impact  on  near-source 
observiuons.  There  is  on  indication  in  the  scaled  dau  that  Ihe  one  shot 
from  below  me  water  ubie  has  higher  accelerations  man  me  other 
explosions.  Synthetic  near-source  reflectivity  calculations  were 
investigated  for  realistic  plane-layered  Pahute  Mesa  structures.  Velocity 
contrasts  at  depth  near  me  explosion  pve  body  wave  amplitude 
nuctuations  as  large  as  a  factor  of  2  directly  above  die  explosion.  These 
diffetences  decreased  wim  increasing  soiuce-reeeiver  offsets.  Despite 
strong  velocity  contrasts  near  me  source  the  ground  motion  fluciuauans 
did  not  approach  the  factor  of  6-g  observed  in  the  data.  The  observational 
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dau  argue*  that  lateral  vahaiiona  in  the  velocity  itnieture  (weathered 
xone).  not  taken  into  account  by  the  synthetic*,  nuy  have  a  strong 
concribution  to  the  ground  motions. 

The  scatter  in  the  free  surface  dau  from  Pahute  Mesa  is  apparently 
frequency  dependent  as  exemplified  by  the  determination  of  long-period 
momenL  Spectral  moment  estimates  sJxiwed  greatly  reduced  variation  for 
a  single  explosion,  a  factor  of  2-3.  Free-surface  dau  document  the 
subility  of  the  source  comer  frequency  followed  by  a  second  higher 
frequency  comer  which  decreases  in  frequency  with  increasing  source- 
receiver  offset.  This  second  comer  is  aonbuied  to  attenuation  and  gives 
an  average  Q  value  for  Pahute  Mesa  between  20  and  30.  At  small,  near¬ 
source  ranges  whem  the  source  and  attenuation  comers  separate  the  dau 
indicates  a  f~^  high-frequency  source  model.  At  greater  source-receiver 
ranges  source  and  attenuation  comers  are  diiTicult  to  separate  and  cloud  the 
f~^  source  intetpreuuon. 

Dau  from  Yucca  Flats  were  used  to  illusme  the  relative  deviatoric  and 
isotropic  source  contributions.  Simple  spectral  interpreutions  of 
observations  at  source-receiver  distances  of  2  km  or  less  give  deviatoric 
source  estimates  horn  transverse  motions  that  are  5-10  limes  smaller  than 
isotropic  source  estimates  from  vertical  and  radial  motions.  Transverse 
spectra  interpreted  as  a  Brune  type  model  give  a  deviatoric  source  radius 
which  is  at  the  lower  bound  of  the  equivalent  elastic  source  radius  for  the 
explosion.  As  found  for  the  Pahute  Mesa  data,  spectral  differences 
between  the  radial,  vertical  and  transverse  spectra  decrease  with  range.  At 
5  km  all  three  comoonenis  of  motion  give  spectra  that  are  identical  in 
shape  and  strength,  indicating  the  strong  effeeu  of  scattering  and 
attenuation. 

Full  wave  modeling  via  moment  tensor  inversion  are  compared  to 
simple  spectral  inierpreuiion  of  free  surface,  near-source  data.  These 
.  resulu  inrlicaie  that  the  scalar  spectral  inierpretatitm  for  isotropic  moment 
may  be  biased  high  by  as  much  as  a  factor  of  3.  This  bias  is  a  result  of 
the  free  surface  interaction  and  the  generadon  of  surface  waves  not  taken 
into  account  in  the  scalar  moment  esdmates. 

Dau  &oin  the  dee-field  and  free  surface  have  been  used  to  constrain  the 
explosion  source  function.  Contributions  from  the  isotropic  explosion, 
tectonic  stress  release  or  driven  modons  and  spall  have  been  documented. 
Dau  from  Rainier  Mesa  indicate  the  superiority  of  free-ficld  dau  in 
making  source  esdmates.  Spall  zone  motions  offer  opportunity  for 
constraining  this  secondary,  seismic  source  and  indicate  a  strong 
geometrical  effect.  At  near-source  distances  strong  trade-offs  can  develop 
between  source  esdmates  and  the  effects  of  attenuation  and  scattering. 
The  advanuge  of  near-source  dau  sets  is  the  ease  at  which  experiments 
can  be  designed  to  separate  these  source  and  propagauon  effects  as 
documented  by  this  review. 
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VARIABILITY  OF  NEAR  SOURCE  WAVEFORMS 
FROM  CONTAINED  EXPLOSIONS 


Pahute  Mesa:  A  Case  Example 
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ABSTRACT 


The  importance  of  simple  elastic  (with  attenuation)  wave 
propagation  effects  on  near  source  observations  (0.5-10  km)  of 
seismic  waves  generated  by  underground  explosions  is  considered. 
Specifically,  the  study  is  designed  through  a  numerical  exercise  to 
answer  the  question  of  how  well  the  geological/geophysical 
structure  of  a  test  site  must  be  known  in  order  to  make  good 
estimates  of  the  explosion  source  function.  This  study  has 
ultimate  implications  on  the  energy  radiated  to  regional  and 
teleseismic  distances.  Modeling  of  these  wavefields  is  completed 
using  the  extended  reflectivity  method  for  plane  layered  geologies. 
This  method  is  utilized  because  of  the  necessity  to  include  both 
body  and  surface  wave  effects.  Full  wave  velocity  synthetics  are 
generated  and  quantified  by  peak  amplitude  and  a  new  energy 
measure  which  is  defined  as  the  integral  of  the  squared  velocity. 
Four  specific  effects  are  quantified  in  the  study:  (l)Source- 
receiver  offset  effects:  (2)Source  depth  of  burial  effects;  (3)Near 
source  interface  effects;  and  (4)  ’Average’  vs  'site  specific' 
velocity  model  effects.  The  offset  effects  were  found  to  be  quite 
strong  in  these  models  with  spatial  decays  of  both  body  and 
surface  waves  exceeding  spherical  and  cylindrical  decay 
respectively.  Body  waves  are  found  to  decay  between  and  r 
■*•9  and  surface  waves  between  and  ri  *^.  Body  waves  with 
frequencies  greater  than  1  Hz  showed  strong  variations  which  did 
not  fit  these  power  laws  in  the  0-2  km  range.  Spatial  changes  in 
energy  over  the  0.5  to  10  km  range  were  1000-2000  for  the  body 
waves  and  75-160  for  the  surface  waves.  Sources  above  and  below 
interfaces  were  found  to  have  the  second  strongest  effect  on 
amplitudes.  High  frequency  body  waves  exhibited  factors  of  4 
change  in  amplitude  as  sources  moved  into  and  below  a  low 
velocity  zone.  These  differences  disappeared  as  source  to  receiver 
offsets  increased.  Depth  of  burial  effects  were  large  for  sources 
at  200  and  616m.  These  differences  decreased  with  source  to 
receiver  offset  also.  For  reasonable  source  depth  estimates  it  was 
determined  that  this  effect  was  small  as  long  as  the  source 
remained  in  a  smoothly  varying  media.  The  comparison  between  an 
average  and  site  specific  velocity  model  with  a  constant  source 
depth  showed  small  variations  in  peak  amplitude  and  total  energy. 
The  amplitude  and  energy  variations  quantified  in  these 
calculations  are  smaller  than  the  variations  seen  in  near  source 


data  sets.  Lateral  variations  in  the  gei 
coupling  effects  may  account  for  additi 
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MOTIVATION 


Near  source  ground  motions  from  underground  nuclear  explosions 
offer  a  diagnostic  that  can  be  utilized  in  the  constraint  of  physical 
processes  in  the  seismic  source  region.  These  source  processes 
include  but  are  not  limited  to  the  nonlinear  response  of  the 
material  in  which  the  explosion  is  detonated,  the  release  of  stored 
tectonic  stress  in  and  around  the  explosion,  and  the  interaction  of 
the  explosive  wavefield  with  the  free  surface  and  other  interfaces. 
Waveforms  are  affected  by  all  of  these  processes  as  well  as  the 
geological  media  through  which  the  seismic  waves  propagate.  In 
any  study  attempting  to  either  resolve  the  explosive  source 
function,  discriminate  the  resulting  wavefield  from  that  developed 
by  a  naturally  occurring  earthquake,  or  constrain  the  yield  or  size 
of  the  explosion  one  must  be  able  to  separate  the  wave  propagation 
effects  from  the  source  itself. 

The  purpose  of  this  study  is  to  begin  to  quantify  the  role  simple 
near  source  wave  propagation  effects  play  in  the  observed 
wavefields  separate  from  source  issues  such  as  coupling.  We  shall 
focus  upon  wave  propagation  at  Pahute  Mesa  in  the  0-10  km  range 
and  quantify  the  variability  of  amplitudes  that  might  be 
experienced  within  the  constraints  of  current  geological  and 
geophysical  models  of  this  area  of  the  Nevada  Test  Site.  Through 
this  numerical  exercise  we  would  like  to  answer  the  question  of 
how  well  the  geological  structure  must  be  known  in  order  to  make 
good  estimates  of  the  explosion  source  function.  A  secondary  but 
longer  term  goal  of  this  study  is  to  investigate  the  effects  that 
this  near  source  structure  has  on  energy  observed  at  regional  and 
teleseismic  distances.  This  second  goal  is  beyond  the  scope  of  the 
current  study. 

We  will  first  address  effects  of  source-receiver  offset  in  the 
radiated  wavefield.  The  purpose  of  this  part  of  the  study  will  be  to 
identify  distance  ranges  in  which  the  wavefield  exhibits 
substantial  changes  and  those  in  which  it  does  not.  Strong 
gradients  in  the  radiated  energy  will  translate  into  large 
variations  of  amplitudes  and  possible  sources  of  error  in  analysis. 

The  effect  of  source  depth  of  burial  will  be  addressed.  Relative 
coupling  of  energy  into  body  and  surface  waves  rust  be  quantified 


before  the  adequacy  of  either  or  both  phases  as  diagnostics  for 
source  strength  or  character  can  be  judged. 

Pahute  Mesa  is  characterized  by  rhyolitic  lava,  bedded  tuffs,  and 
ash  flow  tuffs  ranging  from  partially  to  densely  welded.  This 
general  structure  leads  to  interfaces  which  can  in  some  cases  be 
represented  by  relatively  steep  increases  or  decreases  in  velocity. 
In  view  of  this  geological  model  we  will  investigate  the  wave 
propagation  effects  of  interfaces  near  explosive  sources. 

It  is  hoped  that  this  study  of  the  'geometrical'  effects  of  wave 
propagation  at  Pahute  Mesa  (by  'geometrical'  we  mean  elastic  wave 
propagation)  will  help  in  ascertaining  the  usefulness  of  near 
source  data  in  constraining  the  source  physics  from  explosive 
sources  and  possibly  the  yield  of  the  nuclear  device.  In  light  of 
this  last  goal  a  comparison  of  an  average  velocity  model  for  Pahute 
Mesa  will  be  made  against  a  site  specific  model  developed  through 
downhole  exploration  studies. 

The  strongest  motivation  for  a  numerical  study  such  as  this  one  is 
the  observational  data  itself.  Figure  1  is  a  compendium  of  scaled 
peak  accelerations  observed  from  a  number  of  Pahute  Mesa  nuclear 
explosions  plotted  against  scaled  range.  There  is  an  approximate 
factor  of  ten  scatter  in  peak  accelerations  at  any  given  range.  The 
physical  reason  for  this  scatter  could  either  be  the  effect  of 
different  wave  propagation  (depth  of  burial,  proximity  to 
interfaces,  or  receiver  structure)  or  source  coupling.  Sources 
below  the  water  table  yield  the  highest  amplitudes  at  most  ranges 
(solid  squares  in  Figure  1)  suggesting  a  coupling  effect.  This  data 
can  be  contrasted  to  the  variation  of  peak  accelerations  from  a 
single  explosion  given  in  Figure  2.  Although  the  data  density  is 
reduced,  there  is  still  a  large  amount  of  scatter  in  the  data 
indicative  of  wave  propagation  differences  from  this  single  source. 
The  scatter  in  the  data  decreases  at  longer  periods.  From  the  same 
data  as  used  in  Figure  2  acceleration  spectra  were  determined  and 
a  long  period  spectral  level  chosen.  This  level  was  then  converted 
to  moment.  The  average  moment  in  this  case  is  3.8  x  1023  dyne-cm 
with  a  multiplicative  error  of  1.49.  Scatter  in  moment  values  is 
reduced  from  that  of  peak  accelerations  (Figure  3).  Even  in  this 
data  set  there  is  an  indication  of  local  receiver  effects  as  stations 
with  high  moments  from  vertical  waveforms  also  yield  high 
moments  interpreted  from  the  radial  and  transverse  motions 
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PAHUTE  MESA  SCALED  ACCEL 
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(Figure  4).  The  multiplicative  errors  from  the  radial  and 
transverse  data  were  1.58  and  1.79  respectively. 

The  approach  taken  in  this  study  is  to  bound  the  variability  in  near 
source  waveforms  expected  from;  (1)  Offset  effects;  (2)  Depth  of 
burial  effects;  (3)  Interface  effects;  and  (4)  'Average'  versus  'site 
specific'  geological/geophysical  structure. 


Z-MOMENT  (10**23) 


SINGLE  SHOT  MOMENTS 
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METHODOLOGY  AND  MODELS 


A  number  of  researchers  (Vidale  and  Helmberger,  1988;  Johnson, 
1989;  Stump  and  Johnson,  1984)  have  noted  the  importance  of  both 
body  and  surface  waves  in  near  source  waveforms.  Although  the 
distances  over  which  observations  are  made  are  relatively  short 
(0-10  km,  this  study),  the  shallow  source  depths(200-1000m)  lead 
to  these  effects.  It  is  important  to  develop  a  computational 
technique  that  includes  both  body  and  surface  waves.  Focus  in  the 
initial  part  of  this  study  will  be  on  offset  effects,  depth  effects, 
and  interface  effects  all  of  which  can  be  addressed  with  a  plane 
layered  structure. 

With  the  listed  constraints  and  goals  the  extended  reflectivity 
method  was  chosen  for  seismogram  synthesis  (Muller,  1985).  This 
methodology  allows  for  the  determination  of  the  complete 
response  of  plane  layered  structures  and  includes  effects  of  both 
body  and  surface  waves.  In  some  instances  where  waveforms  are 
dominated  by  vertically  propagating  energy  the  reflectivity  method 
can  be  extended  to  include  separate  source  and  receiver  structures 
(Kind,  1985). 

I 

The  choice  of  appropriate  geological/geophysical  model  for  Pahute 
Mesa  lies  at  the  heart  of  this  exercise.  Two  different  models  were 
chosen  in  order  to  replicate  variability  in  waveforms  one  might 
expect  from  the  range  of  data  types  available  to  a  seismologist  as 
he  approaches  a  new  testing  area.  The  first  model  was  developed 
from  downhole  logs  recovered  from  an  emplacement  hole.  This 
model  is  represented  by  solid  lines  in  Figures  5a  and  b  (HOLE  ). 
Near  surface  velocities  begin  at  1.0  km/s  increase  to  a  depth  of 
approximately  350m  (3.0  km/s)  where  a  bedded  tuff  is  enco’intered 
and  P  velocities  drop  to  1.0  km/s.  Below  500m  a  more  competent 
material  is  encountered  and  velocities  rise  (3.5  km/s)  until  630m 
where  a  nonwelded  tuff  is  encountered  and  velocities  drop  again 
for  approximately  70  m  to  2.5  km/s.  This  detailed  velocity  model 
is  contrasted  with  a  relatively  smooth  model  (dashed  lines  Figures 
5a  and  b)  developed  from  the  linearized  damped  least-squares 
travel-time  inversion  of  Leonard  and  Johnson,  1987.  Both  velocity 
models  show  a  similar  increase  in  velocity  with  depth  but  low 
velocity  zones  identified  in  downhole  acoustic  logs  are  missing  in 
the  Leonard  and  Johnson  model.  Differences  in  these  models  will 
allow  us  to  quantify  the  effect  of  detail  structure  on  radiated 
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wavefields  and  identify  the  importance  of  interfaces  close  to  an 
explosive  source.  At  a  depth  scaling  of  122m/kti/3^  a  -150  ^ton 
explosion  would  be  emplaced  at  a  depth  of  648  m  close  to  the 
deeper  velocity  decrease  in  the  HOLE  model. 

This  investigation  has  several  shortcomings  and  is  intended  as  the 
start  of  a  more  exhaustive  numerical  exercise  to  quantify 
propagation  effects.  The  most  apparent  problem  is  the  plane 
layered  structure  with  similar  source  and  receiver  effects.  The 
separate  source  and  receiver  algorithm  of  Kind  (1985)  will  allow 
the  quantification  of  some  of  these  effects.  Two  and  three 
dimensional  models  which  include  the  effects  of  scattering  may  be 
called  for  at  a  latter  date.  No  attempt  to  replicate  source  depth 
effects  such  as  those  proposed  by  Mueller  and  Murphy  (1973)  have 
been  included.  Source  coupling  differences  are  ignored  when  there 
is  obvious  data  (Figure  1)  that  supports  its  importance. 

A  characteristic  set  of  waveforms  at  0.5  km  increments  from  the 
HOLE  velocity  model  is  given  in  Figures  6a-d.  The  numerical  data  ' 
are  sampled  at  32  samples/s  and  the  waveforms  are  relative 
velocity  for  a  source  of  constant  moment.  In  Rgures  6a  and  b  one 
can  identify  the  early  arriving  body  waves  followed  by  latter, 
longer  period  surface  waves.  In  order  to  crudely  separate  body  and 
surface  waves  for  this  model  the  data  were  high  and  low  pass 
filtered  at  1  Hz.  The  low  pass  filtered  (4  pole)  vertical 
components  are  given  in  Figure  6c  while  the  high  pass  data  are 
replicated  in  Figure  6d.  In  the  amplitude  analysis  that  follows 
these  two  representations  will  be  used  to  quantify  the  effects  of 
body  and  surface  waves. 

In  order  to  contrast  the  effects  of  the  various  propagation  models 
peak  velocities  were  picked.  Careful  analysis  of  Figures  6a-d 
reveals  rapid  changes  in  waveshapes  and  amplitudes,  particularly 
for  the  body  waves.  As  a  second,  more  stable  estimate  of  radiated 
seismic  energy,  the  integral  of  squared  velocity  was  chosen  as  a 
measure  of  source  strength.  Since  the  waveforms  involve  both 
body  and  surface  waves,  it  cannot  be  taken  as  an  absolute  energy 
estimate.  The  measure  is  defined  as: 

IV2»  f  v2(x,t)dt 


A  representative  set  of  energy  calculations  as  a  function  of  time 
are  given  in  Figures  7a  and  b.  In  this  case  the  source  is  at  a  depth 
of  665m  in  the  low  velocity  tuffs  of  model  HOLE  .  The  limit  as  to 
approaches  infinity  is  taken  as  the  measure  which  will  be  applied 
to  the  observational  data.  The  dispersion  in  the  propagation  path 
effects  can  be  easily  identified  as  at  1  km  the  energy  plateau  is 
reached  in  approximately  2  seconds  while  at  10  km  it  takes  nearly 
18  seconds  for  the  final  static  value  to  be  reached.  These  plots  are 
for  waveforms  which  have  been  low  pass  filtered  at  1  Hz  and  as 
such  they  are  dominated  by  surface  waves.  Plots  have  been  scaled 
by  r"!  to  account  for  cylindrical  divergence  in  surface  waves.  Peak 
values  of  energy  in  both  radial  and  vertical  plots  exhibit  smooth 
decays  as  a  function  of  range.  This  well  behaved  decay  is 
interpreted  as  a  stable  estimate  of  source  strength  with  no  ranges 
at  which  there  are  erratic  variations. 

Similar  energy  plots  but  now  for  high  pass  filtered  waveforms 
which  emphasize  the  body  waves  are  given  in  Figures  8a  and  b. 
These  waveforms  have  been  scaled  by  r2  to  replicate  spherical 
divergence  of  body  waves.  Dominance  of  body  waves  is  illustrated 
by  the  fact  that  peak  energy  values  are  reached  at  early  times  at 
all  range,  at  1 .5s  for  1  km  and  4s  for  the  1 0km  range.  Energy 
values  show  significant  variation  in  the  0.5  to  2.0  km  range.,  Th? 
radial  value  at  the  1.0  km  range  is  a  factor  of  two  larger  than  the 
0.5  km  value  (after  scaling).  These  wave  propagation  variations 
indicate  that  for  short  source-receiver  offsets  one  would  expect 
strong  variations  in  the  radiated  wavefield  from  the  HOLE 
velocity  model.  A  more  smoothly  spatial  varying  wavefield  is 
found  at  the  greater  source-receiver  offsets. 

Comparing  the  body(>1Hz)  and  surface{<1Hz)  wave  results,  it  is 
apparent  that  high  frequencies  or  body  waves  exhibit  strong 
variations  in  the  0-2  km  range  with  surface  waves  or  long  periods 
showing  more  constant  decay  patterns. 
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OFFSET  EFFECTS:  Effects  of  offset  on  amplitude  and  energy  in  a 
sense  quantify  transition  from  upgoing  to  horizontally  propagating 
to  turning  rays  on  seismograms.  Waveforms  displayed  in  Figures 
6a-d  illustrate  the  dominance  of  body  waves  in  the  0. 5-2.5  km 
range  with  development  of  a  strong  surface  wave  contribution  at 
more  distant  ranges.  The  body  wave  packet  at  closest  ranges  has  a 
duration  of  1-2s  while  it  extends  to  4s  at  10  km.  Surface  waves  in 
contrast  extend  beyond  16s  at  the  10  km  range  for  the  HOLE 
(Figure  5)  velocity  model.  As  discussed  previously,  energy 
estimates  support  a  similar  picture  in  terms  of  ground  motion 
duration  with  surface  waves  dispersing  energy  over  a  window 
nearly  four  times  the  duration  of  that  from  the  body  waves  at  10 
km. 

Peak  velocity  estimates  as  a  function  of  offset  for  the  HOLE 
velocity  model  are  summarized  in  Figures  9a  and  b  (>1Hz)  and 
Figures  10a  and  b(<1Hz)  while  equivalent  energy  plots  are  given  in 
Rgures  11  and  12.  These  plots  include  data  from  sources  at  a 
variety  of  depths  which  will  be  discussed  in  detail  latter.  For 
purposes  of  quantifying  offset  effects  numerical  data  for  a  source 
at  616m  (solid  squares  in  all  plots)  will  be  discussed  as  a 
characteristic  data  set. 

High  pass  filtered  velocity  data  have  been  scaled  by  ri  to  account 
for  body  wave  decay  while  low  pass  filtered  data  have  been  scaled 
bv  to  correct  for  surface  decay.  Similarly  energy  plots  are 

scaled  by  r^  for  high  pass  filtered  data  and  r"*  for  low  pass  filtered 
data. 

Both  vertical(Z)  and  radial(R)  high  pass  filtered  velocities  show 
strong  geometrical  effects  surpassing  standard  ri  elastic  decay. 
Beyond  1  km  vertical  velocities  decay  as  ri-5  while  radial 
velocities  decay  more  quickly  at  These  are  elastic  effects 

introduced  by  the  layered  velocity  structure  reproduced  in  Figure  5. 
The  decay  in  amplitudes  is  not  a  smooth  function  as  suggested  by 
the  simple  decay  laws  but  a  good  deal  of  complexity  is  introduced 
by  geological  structure  with  focusing  and  defocusing.  Any  one  peak 
amplitude  estimate  can  diverge  from  the  simple  decay  laws  by  50- 
100%. 
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Long  period  peak  amplitudes  (Figures  10a  and  b)  which  are 
dominated  by  surface  wave  contributions  also  show  spatial  decays 
in  excess  of  simple  rO-5  elastic  half  space  results.  -Beyond  1-2  km 
decay  rates  are  pi -3  for  vertical  surface  waves  and  for  radial 
motions.  Scatter  about  these  decay  laws  is  much  less  than  that 
observed  for  body  waves.  As  indicated  in  the  figures,  variation  of 
any  one  measurement  from  the  decay  law  is  at  most  10-20%.  This 
result  indicates  the  more  consistent  nature  of  long  periods  from 
the  HOLE  structure. 

Energy  estimates  for  both  high  pass  (Rgures  11a  and  b)  and  low 
pass  (Figures  12a  and  b)  filtered  data  are  more  smoothly  varying  as 
a  function  of  offset  than  the  peak  amplitudes.  This  conclusion  is 
particularly  true  for  body  wave  energy  above  1  Hz.  For  both  high 
and  low  pass  filtered  data  there  is  a  change  in  decay  rate  at 
approximately  2  km.  Beyond  this  range  data  decays  at  a  constant 
rate.  At  shorter  ranges  data  decays  in  some  instances  and 
increases  in  others.  There  are  significant  variations  in  energy 
estimates  at  distance  ranges  of  less  than  three  depths  of  burial. 
This  result  was  found  in  the  time  series  displayed  in  Figures  8a 
and  b. 

DEPTH  OF  BURIAL  EFFECTS:  Source  depth  effects  are  best 
illustrated  using  the  smoothly  varying  velocity  model  of  Leonard 
and  Johnson  (Figures  5a  and  b).  Since  peak  amplitude  results 
showed  significant  station  to  station  variability  and  energy 
estimates  utilizing  the  entire  waveform  are  more  consistent 
between  individual  stations,  depth  of  burial  results  will  be 
quantified  with  the  integral  of  velocity  squared.  Frequency  effects 
were  determined  by  high  and  low  pass  filtering  of  the  data  at  1  Hz. 

Major  depth  changes  were  considered  with  sources  at  200m,  400m, 
and  616m.  Low  pass  filtered  vertical  and  radial  peak  energy 
estimates  are  reproduced  in  Figures  13a  and  b  while  high  pass 
filtered  data  are  given  in  Figures  14a  and  b.  Shallow  sources  yield 
larger  energy  estimates  at  the  free  surface  as  determined  from  the 
time  integral  of  ground  velocity  squared.  As  much  as  a  factor  of 
ten  increase  in  energy  existed  for  a  source  buried  at  200m 
compared  to  the  one  at  616m.  Differences  in  energy  estimates 
decreased  with  offset,  reaching  values  close  to  2  at  10km.  These 
source  energy  differences  can  be  compared  to  the  geometrical 
decay  of  energy  over  the  observation  range  of  0.5  to  10  km.  For 
high  frequencies  the  energy  differences  (0.5  to  10  km)  are  between 
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1000-2000  while  for  the  low  frequencies  they  are  between  75- 
160.  In  both  the  high  and  low  pass  data  sets  the  spatial  decay  is 
between  1  and  2  orders  of  magnitude  greater  than  the  energy 
differences  at  a  constant  range  for  sources  at  200,400  and  616m 
depths  These  'energy'  values  are  for  the  integral  of  velocity 
squared.  Taking  the  square  root  of  results  gives  numbers 
comparable  to  velocity  estimates  given  earlier.  Variations 
between  the  400  and  616m  sources  in  the  low  pass  data  then 
become  1.4  at  0.5  km  and  1.3  at  10  km. 

Errors  in  depth  estimation  should  be  much  smaller  than  variations 
introduced  in  these  numerical  trials.  It  is  concluded  that  errors  in 
waveforms  introduced  through  depth  errors  are  small.  In  contrast 
spatial  decay  of  these  effects  is  strong  enough  that  accurate 
spatial  locations  of  receivers  is  important.  Data  in  the  5-10  km 
range  exhibit  the  smallest  variation  with  source  depth.  Data  below 
1  Hz  were  found  to  give  more  stable  estimates. 

INTERFACE  EFFECTS:  The  third  goal  of  this  numerical  exercise  was 
to  quantify  effects  near  source  interfaces  have  on  radiated  seismic 
energy.  The  HOLE  velocity  model  was  chosen  to  explore  this 
variation.  As  indicated  in  Figure  5,  there  is  an  approximate  30% 
decrease  in  P  velocity  at  630m  depth  extending  to  700m.  This 
depth  range  is  close  to  the  scaled  depth  for  a  150  kton  explosion. 
The  velocity  decrease  is  representative  of  velocity  changes  found 
at  Pahute  Mesa  of  the  the  Nevada  Test  Site.  They  are 
representative  of  variations  introduced  by  the  amount  of  welding 
in  the  tuffs  and  transitions  to  rhyolitic  lavas.  In  order  to  quantify 
these  geological  effects  on  seismic  radiation  a  set  of  synthetic 
calculations  were  com^pleted  with  sources  at  616m  (above  low 
velocity  zone),  665m  (in  low  velocity  zone),  and  750m  (below  low 
velocity  zone).  Waveforms  for  the  616m  source  have  been 
reproduced  in  Figures  6a  and  b.  Radial  and  vertical  velocity 
records  for  the  665m  source  are  in  Figures  15a  and  b  while  those 
for  the  750m  source  are  in  Figures  16a  and  b.  Comparison  of  these 
three  sets  of  time  series  identifies  a  dramatic  decrease  in 
amplitudes  for  the  665  and  750m  sources  relative  to  that  at  616m. 
This  decrease  in  amplitude  for  the  two  deeper  sources  is  greatest 
at  closest  offsets  and  is  balanced  at  farther  ranges  by  slower 
spatial  decay  rates  for  deeper  sources. 

These  conclusions  can  be  -further  illustrated  by  comparing  peak 
velocities.  Figures  9a  and  b  reproduce  the  peak  amplitudes  from 
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the  high  pass  filtered  data  {>1Hz)  which  are  dominated  by  body 
waves.  Over  the  0.5-3km  range,  the  616m  source  yields  peak 
amplitudes  that  are  a  factor  of  two  or  more  greater  than  those  for 
the  two  deeper  sources.  Although  there  is  a  great  deal  of  scatter 
in  these  peak  amplitude  estimates,  data  from  ail  three  depths  of 
burial  merge  near  10  km. 

Energy  estimates  were  made  for  the  same  high  pass  filtered  data 
as  well  as  for  a  source  at  200m  DOB  in  the  HOLE  velocity 
structure.  Comparisons  of  these  estimates  are  given  in  Figures 
11a  and  b.  As  with  peak  velocity  data,  energy  estimates  indicate  a 
wide  difference  for  data  from  above  (616  and  200m)  and  below{665 
and  750m)  the  interface  at  the  shortest  offsets.  This  difference 
decreases  with  range.  The  major  effect  in  the  near  source 
waveforms  occurs  as  the  source  moves  below  the  first  interface  of 
the  low  velocity  zone.  The  665m  source  which  is  in  the  low 
velocity  zone  gives  results  which  are  nearly  identical  to  those  for 
the  source  at  750m  which  is  completely  below  the  low  velocity 
zone. 

Long  period(>1Hz)  peak  amplitude  data  are  given  in  Figures  10  a  and 
b.  Differences  between  sources  above  and  below  the  interface 
exhibited  by  the  body  waves  are  greatly  reduced  for  near  source 
long  period  (<1Hz)  surface  waves.  There  is  a  20%  difference 
between  peak  amplitudes  for  all  three  source  depths.  Long  period 
energy  estimates  are  reproduced  in  Figures  12a  and  b  along  with 
estimates  for  the  200  m  source.  There  is  less  than  a  factor  of  2 
difference  in  total  energy  for  the  three  source  depths  with  this 
difference  again  decreasing  with  range.  At  iOkm  this  energy 
difference  is  a  factor  of  1.8  for  the  vertical  motions  ana  1.2  for 
the  radial.  In  units  of  velocity  this  reduces  to  factors  of  1.3  and 
1.1  respectively.  These  small  differences  are  contrasted  by  the 
large  differences  already  discussed  for  the  200m  DOB  source. 

Interface  effects  are  strongest  in  high  frequency  body  waves.  In 
all  cases  the  effect  of  moving  a  source  above  or  below  an  interface 
near  the  scaled  depth  for  a  150  kton  explosion  decreases  with 
range.  The  smallest  variations  in  surface  seismic  radiation  are 
found  for  long  period  motions  (<1Hz)  which  are  dominated  by 
surface  waves.  At  10  km  range  the  square  root  of  the  integral  of 
squared  velocity  is  found  to  vary  by  a  factor  of  1. 3-1.1  as  sources 
are  moved  above  and  below  a  realistic  low  velocity  zone  near  shot 
depth. 
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VELOCITY  MODEL  COMPARISONS:  The  final  issue  that  was  addressee 
in  these  synthetic  tests  was  the  importance  of  a  site  specific 
velocity  model  (such  as  that  developed  from  a  downhole  log  at  the 
test  site)  in  contrast  to  an  average  model  for  a  test  site.  In  order 
to  quantify  this  discussion  synthetic  velocity  records  were 
developed  for  the  'average'  Leonard  and  Johnson  velocity  model  and 
the  'site  specific’  HOLE\  velocity  log.  Sources  for  this  final 
comparison  were  placed  at  a  depth  of  616m.  Vertical  and  radial 
synthetics  were  high  pass  filtered  at  1Hz  (Figures  17a  and  b)  and 
low  pass  filtered  at  1  Hz  (Figures  18a  and  b)  for  comparison.  Peak 
velocity  differences  in  these  models  are  remarkably  small.  The 
biggest  differences  occurs  at  the  smallest  offsets  Radial  data  in 
Figure  17b  for  the  two  velocity  models  nearly  match  one  another 
beyond  the  3  km  range.  The  one  difference  from  this  general  trend 
is  in  the  low  frequency  radial  data  where  the  two  velocity  models 
lead  to  a  diverging  set  of  curves  with  increasing  range. 

Of  the  effects  considered  which  include  source-receiver  offset, 
source  depth,  proximity  to  an  interface,  and  average/site  specific 
velocity  model,  the  velocity  model  seems  to  have  the  least  effect 
on  peak  amplitudes  for  the  target  depth  of  616m. 
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CONCLUSiONS/RSCOMMENDATIONS 


A  program  designed  to  quantify  geometricai  effects  of  near  source 
seismic  wave  propagation  has  begun.  Analysis  is  motivated  by 
observational  data  from  Pahute  Mesa  of  the  Nevada  Test  Site  and  as 
a  result  velocity  models  developed  in  this  study  are  applicable  to 
this  area.  Extended  reflectivity  modeling  is  used  to  synthesize 
seismograms  from  plane  layered  structures  in  the  0.5  -  10.0  km 
range.  Calculations  show  strong  contributions  from  body  and 
surface  waves  over  this  near  source  region.  Surface  waves  are 
particularly  strong  because  of  shallow  burial  depths  (200-750m). 

Data  were  crudely  separated  into  body  and  surface  wave 
contributions  using  high  and  low  pass  filters  centered  at  1  Hz.  A 
general  observation  is  that  body  waves  show  stronger  amplitude 
fluctuations  than  surface  waves  when  peak  velocities  are 
considered.  In  order  to  develop  a  more  stable  estimate  of  source 
strength  from  these  observations  the  time  integral  of  the  squared 
velocity  was  introduced  as  a  new  waveform  measure.  Preliminary 
analysis  of  this  function(roughly  called  an  energy  measure)  with 
the  synthetic  data  indicates  that  it  is  a  more  robust  estimate  of 
source  strength. 

Four  specific  questions  were  addressed  in  this  initial  waveform 
modeling  exercise  using  the  measures  discussed.  These  questions 
dealt  with  the  quantification  of  the  following  effects:  (l)Source' 
receiver  offset;  (2)Proximity  of  interfaces  to  the  source;  (3) 

Source  depth  of  burial;  and  (4)’Average'  vs  'Site  Specific'  velocity 
models. 

The  effect  of  offset  on  the  radiated  wavefield  was  found  to  be  the 
strangest  at  horizontal  distances  of  0  to  3  depths  of  burial.  The  P 
waves  (as  quantified  by  data  above  1  Hz)  were  most  variable  over 
this  range.  Adjacent  stations  separated  by  only  0.5  km  were  found 
to  change  in  amplitude  by  a  factor  of  two.  These  changes  were 
seen  as  both  velocity  increases  and  decreases.  At  greater 
distances  body  waves  decayed  according  to  a  power  law  which 
always  exceeded  simple  geometrical  spreading.  Radial  velocity 
decayed  as  pi-s  while  vertical  velocities  decayed  as  r'l-S. 

Surface  wave  peak  amplitudes  also  showed  some  variation  in 
amplitudes  at  close  source  to  receiver  offsets.  At  ranges  greater 
than  1-2  km  surface  waves  fit  power  law  decay  patterns.  These 
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decays  exceeded  cylindrical  spreading  with  values  of  for 
vertical  components  and  r-i-^  for  radial  components.  Variations  in 
Individual  surface  wave  observations  from  this  decay  were  10- 
20%,  much  smaller  than  that  observed  for  body  waves.  Energy 
estimates  (integral  of  the  square  of  velocity)  exhibited 
dramatically  reduced  scatter  from  station  to  station  when 
compared  to  peak  amplitude  estimates. 

Source  proximity  to  an  interface  was  replicated  with  a  series  of 
sources  at  depths  of  616,  665  and  750  m  in  the  HOLE  velocity 
model  (Figure  5).  The  616m  source  is  above  a  low  velocity  zone, 
the  665m  source  is  in  the  low  velocity  zone,  and  the  750m  source 
is  below  the  zone.  Strong  differences  in  the  body  waves  generated 
by  the  source  above  the  low  velocity  zone  and  those  in  and  below 
the  zone  were  obser/ed.  Over  the  0.5  to  3.0  km  range,  the  616m 
source  yielded  body  waves  above  1  Hz  that  were  a  factor  of  two  to 
four  larger  than  those  for  either  of  the  two  deeper  sources. 
Differences  between  the  data  from  the  three  source  depths 
decreased  with  increasing  offset.  Long  period  waveforms 
represented  by  surface  wave  contributions  were  much  less 
affected  by  source  depth  with  differences  in  the  range  of  20% 
between  individual  peak  amplitudes  for  the  three  source  depths. 

Source  depth  effects  were  quantified  by  considering  sources  at 
200,  400,  and  616m.  Energy  estimates  indicated  that  as  much  as  a 
factor  of  ten  increase  in  surface  energy  for  the  200m  source 
compared  to  the  616m  source.  These  differences  were  greatest  at 
short  source-receiver  offsets  and  decreased  to  a  factor  of  2  at  10 
km.  Source  depth  effects  were  1  to  2  orders  of  magnitude  less 
than  the  change  in  energy  due  solely  to  geometrical  spreading  over 
the  0.5  to  10  km  range.  Variations  between  the  400  and  616m 
sources  were  even  less.  The  conclusion  drawn  from  these  trials 
was  that  for  reasonable  constraints  on  source  depth  variations  in 
radiated  energy  will  be  small  at  large  (3-1 0km)  offsets  for  the 
velocity  model  considered. 

The  final  set  of  trials  focused  on  two  sets  of  synthetics,  one  from 
an  average  Pahute  Mesa  ve'ocity  structure  and  the  second  from  a 
specific  model  developed  with  downhole  information.  Sources 
were  placed  at  616m  in  both  structures.  Differences  between 
synthetics  for  the  two  models  were  small  compared  to  other 
effects  identified  in  these  trials.  Again  the  biggest  differences 
occurred  at  small  offsets. 
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These  trials  have  identified  that  peak  velocity  estimates  in  the 
near  source  region  may  exhibit  strong  lateral  variations  at  close 
offsets  and  high  frequencies  which  are  dominated  by  body  waves. 
These  differences  degrade  when  greater  offsets  or  longer  periods 
(surface  waves)  are  considered.  The  integral  of  squared  velocity  is 
found  to  be  a  stable  source  estimate.  These  new  criteria  may  be 
useful  in  analyzing  observational  data  sets.  Effects  of  interfaces 
and  low  velocity  layers  were  found  to  be  strongest  for  body  waves 
at  short  offsets.  Finally,  with  reasonable  estimates  of  source 
depth  this  error  should  not  be  dominant  in  observational  data  as 
long  as  the  source  is  not  near  an  interface. 

Purely  geometrical  wave  propagation  effects  have  been  identified 
In  these  numerical  exercises  with  no  consideration  of  coupling 
effects.  Strong  effects,  as  big  as  a  factor  of  4,  have  been  found  for 
high  frequency  body  waves  at  short  offsets.  Reduced  effects  are 
found  for  longer  period  surface  waves  at  great  source  to  receiver 
offsets.  These  factors  are  less  than  the  factor  of  10  found  in  the 
observational  data  set.  Additional  scatter  must  be  introduced  by 
lateral  variations  in  the  geological  structure  and  coupling  effects 
from  explosion  to  explosion.  The  fact  that  scatter  in  peak 
amplitude  data  does  not  decay  with  range  as  suggested  by  these 
numerical  trials  also  calls  for  one  of  these  other  mechanisms. 

Many  things  have  not  been  considered  by  these  trials.  Our  next  step 
is  to  complete  ray  tracing  in  the  structure  to  come  to  a  better 
understanding  of  the  reasons  for  the  large  variations  in  body  waves 
at  short  offsets.  The  stability  of  long  period,  long  offset  data 
motivates  a  new  analysis  of  tl  a  observational  data  reproduced  in 
Figures  1-4.  The  reduced  scatter  in  long  period  moment  estimates 
may  be  a  reflection  of  physical  effects  explored  in  these  trials. 
Implementation  of  energy  estimation  to  the  observational  data  set 
is  now  warranted. 

Consideration  of  separate  source  and  receiver  structures  on  near 
source  waveforms  must  also  be  completed.  This  modeling  is  the 
first  step  in  investigating  lateral  variations  in  propagation  path. 
We  intend  to  begin  this  analysis  with  the  method  of  Kind  (1985). 
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Abstract 

Near-source  data  from  the  nuclear  explosion  Coalora  detonated  at  Yucca  Flats, 
Nevada  Test  Site  are  utilized  to  constrain  the  seismic  source  function.  The 
equivalent  seismic  source  is  interpreted  in  terms  of  physical  processes  in  the 
source  region  with  the  aid  of  data  from  within  the  explosion’s  nonlinear 
region.  The  isotropic,  deviatoric,  and  cylindrical  spall  contributions  are 
separated  and  quantified.  Standard  spectral  interpretations  of  the  radiated 
wavefield  for  source  resolution  are  contrasted  against  complete  waveform 
modeling  with  moment  tensor  determination.  Individual  waveform  spectra 
(source-receiver  offsets  <  2  km)  can  be  interpreted  in  terms  of  an  isotropic 
source  model  which  is  in  agreement  with  a  Mueller-Murphy  model, 
including  f  high  frequency  decay  and  source  comer  frequency  of  1.8  Hz. 
Moment  tensor  inversion  produces  an  isotropic  source  strength  of  8x10^0 
dyne  cm  while  scalar  moments  from  the  spectral  interpretation  are  a  factor  of 
2.5  larger.  This  difference  is  attributed  to  the  application  of  simple 
propagation  path  corrections  in  the  spectral  interpretation.  The  deviatoric 
component  of  the  moment  tensor  is  a  factor  of  5-10  times  smaller  than  the 
isotropic  component.  Deviatoric  source  radius,  as  estimated  from  the  spectral 
data,  is  124  m,  smaller  than  the  equivalent  elastic  source  radius  which  is 
boimded  between  133-202  m.  Stress  drop  estimated  with  the  Brune  source 
model  is  84  bars  with  an  average  slip  of  33  cm.  The  spall  source  contribution 
is  longer  period  and  delayed  in  time  fr  om  the  initial  explosion.  Its 
contribution  to  the  diagonal  elements  of  the  moment  tensor  is  dominant  on 
the  Mzz  component  at  least  a  factor  of  3  larger  than  the  Myy  and  h^xx 
components.  Spall  source  strength  from  waveform  inversion  is  within  a 
factor  of  2  of  forward  spall  models  developed  from  acceleration  data  within 
the  spall  zone  but  is  longer  in  duration  than  the  forward  model,  reflecting  the 
effect  of  a  quasi-point  source  assumption  in  the  forward  models.  Complex 
propagation  effects  dominate  the  data  beyond  2  km.  This  complexity  is 
exemplified  by  wavetrains  at  5.1  km  that  extend  to  beyond  20  s  in  duration 
and  similarity  of  radial,  vertical,  and  transverse  acceleration  spectra  in 
contrast  to  short  waveform  durations  and  strong  spectral  dif.'erances  between 
transverse  and  radial-vertical  spectra  at  distances  less  than  2  km.  These 
propagation  path  effects  suggest  that  source  biases  can  develop  at  ranges  as 
dose  as  2-5  km. 


Introduction 


Seismic  source  studies  of  underground  explosions  are  motivated  by  the  need 
to  separate  physical  processes  that  generate  seismic  waves  at  near-source, 
regional,  and  teleseismic  distances.  These  physical  processes  include  but  are 
not  limited  to  direct  coupling  of  the  spherical  explosion,  repartitioning  of 
energy  near  the  free  surface  by  tensile  failure  and  driven  motions  along 
planes  of  weakness,  and  tectonic  stress  release.  Application  of  such  physically 
based  models  to  seismic  data  is  necessary  for  a  complete  imderstanding  of 
discrimination  issues  and  yield  determination  for  the  purposes  of  nuclear  test 
ban  treaty  verification.  A  physically  based  model  of  the  source  provides  the 
methodology  to  extrapolate  from  one  geological  environment  to  another  and 
therefore  can  be  tested  against  new  data  sets. 

The  simplest  characterization  of  the  seismic  source  consists  of  a  spherically 
symmetric  or  isotropic  equivalent  elastic  source  function.  A  variety  of  such 
models  have  been  developed  with  different  parameterizations  and 
constrained  by  not  only  near-source  but  regional  and  teleseismic  data 
(Haskell,  1967;  Mueller  and  Murphy,  1971;  von  Seggem  and  Blandford,  1972; 
Wallace  and  Burger,  1987).  These  characterizations  have  been  extended  to 
chemical  explosions  (Stump  and  Reinke,  1987;  Stump,  1987). 

Additional  phenomena  accompany  the  detonation  of  an  underground 
explosion.  One  of  these  secondary  processes  is  near-surface-layer  tensile 
failure,  spalKEsler  and  Qrilton,  1964;  Eisler  et  al,  1966).  As  the  initial 
spherical  compressive  wave  from  the  explosion  encounters  the  free  surface,  it 
reflects  as  a  tensile  wave  propagating  away  from  the  free  surface.  This  tensile 
wave  may  fail  near-surface  layers.  Momentum  trapped  in  the  layer  results  in 
the  failed  material  following  a  ballistic  brajectory.  This  nonlinear  process  may 
be  responsible  for  repartitioning  part  of  the  explosion  energy  as  seen  in  the 
radiated  seismic  wavefield.  The  importance  of  this  process  as  a  soiirce  of 
seismic  waves  was  addressed  by  Viecelli  (1973).  He  argued  that  spall  had  a 
significant  contribution  to  surface  waves  at  10-20  km.  Day  et  al  (1983) 
presented  a  spall  model  which  conserved  momentum  and  found  that  the 
process  had  no  contribution  to  20  s  surface  waves  at  teleseismic  distances. 


Stump  (1985)  argued  from  data  within  the  spall  zone  that  the  process  made  a 
significant  contribution  to  near-source  smrface  waves  from  chemical 
explosions  although  at  a  much  shorter  period  (.1-1  s)  than  Day's  surface  wave 
analysis.  Patton  (1990)  has  analyzed  near-source  obser/ations  in  the  spall 
zone  at  Pahute  Mesa,  NTS  in  order  to  develop  scaling  relations  for  spall  mass 
and  momentum.  Day  and  McLaughlin  (1991)  show  the  equivalence  of  point 
force  and  moment  tensor  representations  of  the  equivalent  spall  source 
model  which  could  be  developed  from  such  measurements.  Other  authors 
have  argued  that  the  spall  secondary  source  can  have  important  contributions 
to  regional  waves  (Patton,  1988;  Taylor  and  Randall,  1989)  and  teleseismic  P 
waves  (Schlittenhardt,  1991). 

Transverse  shear  motions  and  Love  waves  have  been  observed  from  all  types 
of  explosions  (Kisslinger  et  al,  1961;  Ohta,  Y.  and  A.  Kubotera,  1968;  Aki  and 
Tsai,  1972;  Wallace  et  al,  1985;  Burger  et  al,  1986).  Transverse  motions  have 
been  attributed  to  mode  conversions  and  scattering  (Gupta  and  Blandford, 
1983),  tectonic  strain  release  (Akr  and  Tsai,  1972;  Wallace  et  al,  1985;  Burger  et 
al,  1986),  or  driven  motion  along  planes  of  weakness  (Bache  et  al,  1979). 
Quantification  of  these  complex  propagation  effects  and  deviatoric  source 
contributiorrs  are  important  for  verification  applications  as  they  make  the 
explosions  look  more  earthquake-like  and  may  bias  seismic  yield  estimates. 

Trade-offs  between  the  seismic  source  (it's  strength  and  physical 
characterization)  and  propagation  path  effects  are  possible  in  any  study  of  the 
equivalent  linear  source  representation.  The  focus  of  this  study  is  on  near¬ 
source  data  (<  6  km)  where  propagation  path  effects  can  be  minimized.  Even 
at  these  dose  ranges  a  number  of  authors  have  presented  evidence  for 
propagation  path  complications.  McLaughlin  et  al  (1983)  utilizing  dosely 
spaced  arrays  found  a  degradation  in  coherence  above  5  Hz  for  observations  at 
6  km  from  nuclear  explosions.  Vidale  and  Heimberger(1987)  give  evidence 
for  two  dimensional  structural  effects  on  near-source  waveforms  from 
nuclear  explosions  at  the  Nevada  Test  Site  (NTS). 

The  seismic  experiment  on  the  Coalora  event  (Yucca  Flats,  NTS)  which  we 
report  was  designed  to  quantify  both  the  deviatoric  and  isotropic  parts  of  the 
equivalent  elastic  nudear  source  in  partially  saturated  tuffs  and  alluvium. 
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primary  materials  in  which  the  United  States  tests  nuclear  devices  (US,  Dept 
of  Energy,  1987).  Several  investigators  have  studied  near-source  ground 
motion  from  explosions  in  these  materials  (Werth  and  Herbst,  1962;  Ferret 
and  Bass,  1976).  The  Yucca  Flats  test  site  was  chosen  because  of  the  vast 
quantity  of  geological  and  geophysical  site  characterization  information. 
Ferguson  et  al  (1988)  studied  the  basin  structure  under  Yucca  Flats  with  high 
resolution  seismic  reflection  and  gravity  data.  The  geologic  structure  close  to 
the  experiment  is  well  constrained  by  many  downhole  geological  and 
geophysical  measurements  (Figure  2).  As  a  result,  structural  effects  can  be 
taken  into  account  close  to  the  source  (<  2  km). 
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Experimental  Design 


The  Coalora  experiment  had  an  announced  yield  of  "less  than  20  kilotons 
(kt)"  (US  Dept  of  Energy,  1987).  The  working  point  (WP)  depth  of  the 
explosion  was  274  m  in  Yucca  Eats  (Figure  1).  A  great  variety  of  geological 
and  geophysical  information  was  available  to  constrain  the  wave  propagation 
model.  Figure  2  summarizes  these  data  which  include  the  surface  expression 
of  Yucca  Fault,  location  of  a  seismic  reflection  survey,  downhole  acoustic  logs, 
and  locations  of  available  geologic  cross  sections  such  as  that  developed 
around  the  Coalora  emplacement  hole  (Figure  1). 

The  Yucca  Eats  structure  summarized  in  the  inset  of  Figure  1  is  based  upon 
the  seismic  reflection  survey  and  inversion  of  available  gravity  data 
(Ferguson  et  al,  1988).  The  model  consists  of  layers  of  alluvium,  tuff,  wet  tuff, 
and  Paleozoic  sedimentary  rocks.  Paleozoics  are  near  the  surface  on  either 
side  of  the  10  km  wide  valley.  At  the  center  of  the  valley,  where  the  shot  was 
detonated,  depth  to  the  Paleozoics  is  1  km.  The  Coalora  cross  section  in 
Figure  1  gives  more  shot  region  detail.  Major  geologic  characteristics  include 
Yucca  fault  with  a  maximum  150  m  of  displacement  within  300  m  of  the 
working  point,  interbedded  tuffs,  and  500  m  water  table.  Large  velocity  and 
density  contrasts  e.xist  within  the  interbedded  tuffs.  A  downhole  acoustic  log 
333  m  NW  of  Coalora  (Figure  3a)  indicates  near-surface  velocity  of  1.5  km/s 
with  a  positive  gradient  to  2.0  km/s  at  400  m  .  Below  400m,  the  vitrified 
Rainier  Mesa  tuff  reaches  a  velocity  of  3.0-3.5  km/s.  This  particular  log  does 
not  penetrate  below  this  material.  A  second  log,  U31a,  (Figure  3b),  1.5  km  N  of 
Coalora,  shows  similar  effects.  At  this  location,  two  vitrified  tuff  units,  the 
Ammonia  Tanks  between  210-275  m  and  the  Rainier  Mesa  between  335-40o 
m,  are  apparent. 

The  instrumentation  array  within  two  depths  of  burial  of  the  Coalora  ground 
zero  is  illustrated  in  Figure  4.  These  gages  are  used  to  constrain  the  depth  and 
range  of  spall  and  imp.ove  resolution  of  spall  contributions  to  near-source 
waveforms.  Included  in  this  array  are  two  downhole  packages  in  the 
emplacement  hole  for  the  purposes  of  constraining  the  depth  of  spall.  The 
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Figure  1:  Geological  aoss  section  passing  through  the  emplacement  hole  for 
the  Coalora  explosion.  The  general  Yucca  Flats  velocity  model  as  determined 
by  Ferguson  et  al,  1988  is  given  in  figures  lower  left  hand  comer. 
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Figure  2:  Map  summarizing  geological  information  available  for  constraining 
structure  in  and  around  the  Coalora  test  site.  The  SMU  cross  section  is  based 
upon  a  reflection  survey,  stars  represent  geologic  cross  sections  based  upon 
emplacement  or  exploratory  holes,  and  DHAL  are  dry  hole  acoustic  logs. 
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Figure  3b:  Dry  hole  acoustic  log  located  ISOCm  N  of  Coalora. 
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abundance  of  data  at  different  azimuths  and  ranges  from  the  source  provides 
opporttinity  for  improved  resolution  of  the  deviatoric  source  contribution. 

All  13  close-in  accelerometers  were  hardwired  to  a  central  recording  site.  Data 
were  recorded  on  high  speed  analog  tape  recorders  and  digitized  post-test. 

The  second  part  of  the  instrumentation  array  consisted  of  three-component 
force-balance  accelerometers  with  digital  event  recorders  placed  between  0.86 
and  5.10  km  (Figure  5).  A  total  of  20,  three-component  stations  were  installed. 
The  data  were  sampled  at  200  samples  per  second.  Data  words  were  each  12 
bits. 

T.  -  two  =»ccelerometer  arrays  were  designed  to  provide  data  over  a  range  of 
e...  -uths  and  distances  such  that  contributions  from  the  spherical  explosion, 
th  a  cylindrical  spall,  and  the  deviatoric  source  components  could  be 
separately  resolved. 
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Kgiure  5:  Three-component  force-balance  accelerometers  recorded  with  digit: 
event  recorders.  Stations  span  the  boxmds  of  Yucca  Hats  with  closest  station* 
0.86  km  and  most  distant  5.16  km. 


Observational  Data 


Vertical  accelerograms  from  within  the  spall  zone  are  reproduced  in  Figure  6. 
Records  are  characterized  by  arrival  of  the  compressive  wave,  minus-one-g- 
dweU  during  free-fall,  and  impulsive  rejoin.  The  closest  downhole  gauge  at  a 
depth  of  114  m  shows  evidence  of  multiple  spall  while  the  most  distant  free- 
surface  gauge  that  showed  a  spall  signature  was  at  a  horizontal  range  of  274  m 
(a  radial  distance  equal  to  one  depth  of  burial).  Spall  e.xtends  to  just  beyond 
one  depth  of  burial  in  the  Coalora  experiment.  Two  delineations  of  the  total 
spall  zone  are  given  by  the  solid  and  dashed  lines  in  Figure  6.  Limited  data 
existed  for  defining  the  depth  of  spall  directing  below  ground  zero  (GZ)  and 
no  data  existed  for  spall  depth  beyond  GZ.  The  depth  of  spall  below  GZ  was 
placed  conservatively  just  below  the  deepest  gage  that  spalled.  The  boundary 
of  the  spall  zone  at  other  distances  was  taken  in  its  simplest  form  as  a  linear 
line  connecting  the  deepest  spalled  gage  and  the  first  free  surface  range  where 
spall  was  not  observed.  Based  on  more  complete  spall  data  sets  from 
chemical  explosions  (Stump,  1985),  spall  zone  shape  is  found  to  be  more  bowl 
like,  thus  the  more  extensive  dashed  line  model.  The  exact  shape  of  the  spall 
zone  is  dependent  upon  the  material  properties  around  the  explosion,  source 
depth  and  size.  Due  to  a  lack  of  complete  azimuthal  coverage  in  the  spall 
zone  the  process  was  assumed  cylindrically  symmetric  Records  from 
geophones  at  a  horizontal  distance  of  93  m  and  four  azimuths  separated  by 
90°  yielded  arrival  times  which  were  identical  within  3  ms.  This  scatter  is 
within  measurement  error  and  so  supports  at  least  an  initially  cylindrical 
pulse  from  the  explosion.  These  gages  were  clipped  after  the  first  arrival  and 
cannot  be  used  to  estimate  symmetry  from  latter  arriving  phases. 

Observational  data  within  two  depths  of  burial  are  summarized  in  Figure  7. 
Vertical  (Z),  radial  (R),  and  transverse  (T)  displacement  waveforms  at  three 
549  m  range  stations  are  displayed.  Displacements  result  from  twice 
integrating  the  accelerograms.  Waveshapes  for  the  radial  and  vertical 
displacements  at  the  different  azimuths  are  very  similar  although  peak 
amplitudes  vary  by  30%  about  the  mean.  Relative  to  the  vertical  and  radial, 
transverse  displacements  are  shorter  in  duration,  exhibit  larger  variations  in 
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COALORA  DISPLACEMENT 
AT  549m  RANGE 


Figure  7:  Displacemer.t  waveforms  (accelerations  integrated  twice)  at  three 
azimuths  from  somrce  and  at  a  range  of  549  m.  Vertical(Z),  radiaKR),  and 
transverseCT)  displacements  are  given.  Peak  displacements  in  cm. 


amplitude,  and  change  first  motion  as  a  function  of  azimutli.  Large 
amplitude,  relatively  simple  transverse  motions  are  observed  at  ranges  equal 
to  two  depths  of  burial.  Displacement  spectra  were  estimated  for  the  vertical 
and  transverse  records  at  the  549  m  range  (Figure  8).  Envelope  functions 
were  fit  to  each  spectrum  with  a  long  period  level,  comer  frequency,  and  high 
frequency  decay.  Transverse  spectra  (549m  range)  exhibit  higher  comer 
frequencies,  more  rapid  decay  of  high  frequencies,  and  lower  long  period 
levels  than  the  vertical  spectra.  These  spectral  parameters  will  be  used  as 
source  model  constraints. 

Accelerograms  recorded  at  more  distant  ranges  (0.36-5.10  km)  are  summarized 
in  Figures  9a(Z),  9bCR),  9c(r).  Four  representative  ranges  (0.86, 1.66, 3.38,  and 
5.10  km)  have  been  chosen  from  the  20  instruments  fielded.  The  effect  of  the 
two-  and  three-dimensional  structure  of  Yucca  Flats  (Figure  1)  becomes 
apparent  in  these  waveforms.  At  5.10  km  range,  near  the  valley  edge,  all 
three  components  of  motion  exhibit  a  complex  wave  train  with  duration  of 
over  20  s.  Displacement  records  and  spectra  for  the  5.10  km  data  are  given  in 
Figure  10.  Unlike  the  simple  waveforms  at  549  m  with  obvious  spectral 
differences  between  vertical  and  transverse  components,  all  three 
components  at  this  range  have  nearly  identical  spectral  shapes  and 
amplitudes.  The  resulting  source  comer  frequency  one  would  interpret  from 
the  data  is  reduced  from  that  observed  at  the  closer  ranges,  especially  on  the 
transverse  components.  The  5.10  km  waveforms  and  spectra  are  evidence 
that  more  distant  observations  are  strongly  influenced  by  complex  geological 
structure  underneath  Yucca  Flats.  Simple  spectral  source  interpretations 
would  be  severely  biased  at  this  more  distant  data. 

The  effect  of  the  valley  structure  around  Coalora  was  quantified  with  two  and 
three  dimensional  finite  element  calculations  by  Wojdk  and  Vaugn  (1984). 
Results  of  these  calculations  show  asymmetries  in  the  radiated  wavefield 
beyond  1-2  km.  Transverse  motions  in  these  calculations  attributable  to 
diffraction  by  the  Yucca  Fault  were  20-30%  the  size  of  radial  and  vertical 
motions.  Since  the  main  purpose  of  this  study  is  to  constrain  the  physical 
parameters  of  the  seismic  source  function,  primary  emphasis  vrill  be  placed 
on  observational  data  well  away  from  boundaries  of  the  valley.  Modeling 
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COALORA  DISPLACEMENT  SPECTRA 
AT  549ni  RANGE. 


Figure  8:  Spectral  estimates  from  vertical(Z)  and  transverse(T)  waveforms  in 
Figure  7.  Original  time  series  were  windowed  with  a  10%  cosine  taper.  A 
simple  spectral  model  consisting  of  a  constant  long  period  level,  a  high 
frequency  decay  (f”),  and  a  comer  frequency  (f^  was  fit  to  the  data.  Light  line 
in  each  spectrum  is  this  model  with  arrow  representative  of  f^. 


COALORA  VERTICAL  ACCELERATION 


Figure  9a,b,c  Vertical(a),  radial(b),  and  transverse(c)  accelerograms  spanning 
the  0.86-5.10  km  range. 
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COALORA  TRANSVERSE  ACCELERATION 


COALORA  DISPLACEMENTS 
AT  5.16  km 


Fsgiire  10:  VexticaKZ),  radialQl),  and  fcransverseCT)  displacements  and 
displacement  spectra  from  data  at  5.16  km  range.  Z  and  R  spectra  have  been 
arbitrarily  scaled  by  factors  of  100  and  10  so  that  all  three  could  be  compared 
on  a  single  plot 
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and  interpretation  will  be  restricted  to  observations  within  2  km  of  the 
explosion. 
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Standard  Temporal  and  Spectral  Source  Interpretations 

The  simplest  interpretation  of  near-source  seismic  data  involves  utilization 
of  spectral  measures  such  as  long  period  level  for  seismic  moment  and  comer 
frequency  for  source  dimension.  In  order  to  test  these  measures  with  the 
Coalora  data  set,  spectra  such  as  those  displayed  in  Figvre  8  were  used  to 
estimate  long  period  spectral  level,  comer  frequency,  and  high  frequency 
decay  (Table  1).  Radial  and  vertical  components  were  assumed  to  be 
dominated  by  the  isotropic  source  component.  Transverse  motions  were 
assigned  to  the  deviatoric  source.  Based  on  the  valley  size  beneath  Yucca 
Flats,  increasing  waveform  complexity  with  range,  growing  homogeneity  of 
RT2  spectra  with  range,  and  2/3D  simulations  of  Wojdk  and  Vaugn  (1984) 
source  analysis  was  limited  to  data  within  2200  m  of  the  explosion  working 
point.  Moment  estimates  for  isotropic  and  deviatoric  sotarce  contribution  are 
given  in  Figure  11; 

Mi  a  4jfpRa2Qo 

Md  =  4jtpRj3^!;2o  (15) 

Density  (p)  was  taken  from  emplacement  hole  data  at  shot  depth  as  1.8 
gm/cm3.,  compressional  velocity  (a)  as  1.9  km/s,  and  shear  velocity  (p)  as  1.16 
km/s  which  implies  a  Poisson's  ratio  (v)  of  0.2.  Do  is  long  period 
displacement  spectral  level.  The  high  shear  velocity  was  chosen  so  that 
deviatoric  moments  would  be  an  upper  bound  relative  to  isotropic  moments. 
Following  the  methodology  of  Archuieta  et  al.  (1987)  average  log  moments 
and  their  standard  deviations  were  computed  according  to: 

log  (Mq)  =  [(1/ns)  ^  log  (Moi)] 

M  (2a) 


ns 

std  dev  aog  (Mo»  =  [(-!-)  X  {log  Moi  -  log  {Mo))1 
ns-1 


2,V2 


(2b) 
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Table  la:  Deviatoric  Estimates 

Range  DC  Slope  Momenl 

(m)  (m-s)  (Hz)  (dyne-cm) 


GM2T 

274 

0.37 

GM4T 

365 

0.34 

GM6T 

389 

0.32 

GM7T 

612 

0.09 

GM8T 

614 

0.10 

GM9T 

612 

0.30 

GMIOT 

307 

0.49 

GMllT 

495 

0.24 

SIT 

893 

0.10 

S3T 

940 

0.10 

S5T 

1692 

0.07 

S7T 

1663 

0.06 

S8T 

1722 

0.10 

SllT 

2117 

0.10 

3.4 

3 

3.6x10^° 

3.3 

2"^ 

4.4x10^° 

3.5 

3 

4Jxl0^° 

4.0 

4 

Z0xl0^° 

3.5 

2"^ 

2.1x10^° 

3.5 

3^ 

6.4x10^° 

3.5 

3' 

5.3x10^° 

3.7 

3*^ 

43x10^° 

3.5 

4 

3.1x10^° 

2.8 

3 

3.3x10^® 

3.4 

4 

43x10^° 

2.6 

2“^ 

3.8x10^° 

3.0 

4 

6.1x10^° 

2.2 

3 

73x10^° 

Table  lb:  Isotropic  Estimates 


Station 

Range 

(m) 

DC 

(cm-s) 

(Hz) 

Slope 

GM7Z 

612 

0.30 

1.8 

2 

GM7R 

612 

0,36 

Z1 

2* 

GM8Z 

614 

0.30 

1.6 

2 

GM8R 

614 

0.36 

1.8 

2 

GM9Z 

612 

0.34 

1.8 

2^ 

GMllZ 

495 

0.34 

ZO 

2 

GMllR 

495 

0.90 

1.4 

2 

SIZ 

893 

0.20 

1.9 

2 

SIR 

893 

0.38 

1.6 

2-^ 

S3Z 

940 

0.24 

2.0 

3"^ 

S3R 

940 

0.40 

ZO 

3 

S5Z 

1692 

0.21 

1.4 

2 

S5R 

1692 

0.20 

1.7 

2+ 

S7Z 

1663 

0.16 

1.6 

2* 

S7R 

1663 

0.26 

1.7 

2"^ 

S8Z 

1722 

0.11 

1.7 

2'*' 

SllZ 

2117 

0.12 

1.9 

3 

SllR 

2117 

0.10 

Z1 

3 

Moment 

(dyne-cm) 

2.8x10^^ 

3.4x10^^ 

2.8x10^^ 

3.4x10^^ 

32x10^^ 

2.6x10^^ 

6.9x10^^ 

2.8x10^^ 

53x10^^ 

33x10^^ 

5.8x10^^ 

53x10^^ 

53x10^^ 

4.1x10^^ 

6.7x10^^ 

2.9x10^^ 

3.6x10^^ 

3.3x10^^ 
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Rgure  11:  Isotropic(radiai  and  vertical)  and  deviatoric{transverse)  moment 
estimates  for  Coalora. 
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Mean  isotropic  and  deviatoric  logarithmic  moments  from  observational  data 

were  21.29  (standard  deviation  of  mean  0.14)  and  20.60  (s.d.m.  0.16) 
respectively.  These  estimates  for  Mj(isotropic)  and  Mj(deviatoric)  moments 

are  1.95  x  10^^  and  3.97  x  10^^  dyne-cm.  The  isotropic  scalar  moment  is  a 
factor  4.91  larger  than  the  deviatoric  scalar  moment.  The  multiplicative  error 
term  is  detined  according  to: 

EMo  =  antilog  [std  dev  (log  (Mo))]  (3) 

The  multiplicative  error  term  for  the  isotropic  and  deviatoric  moments  are 
1.36  and  1.44.  The  standard  deviations  (std  dev)  of  the  two  means  are  nearly 
equal. 

Comer  frequencies  were  estimated  by  fitting  an  envelope  function  to  the  data 
(Figure  8).  These  comer  frequencies  are  displayed  in  Figure  12  as  a  function  of 
source-receiver  separation.  Transverse  comer  frequencies  are  significantly 
larger  than  radial  or  vertical  estimates  at  close  range.  This  difference 
decreases  with  increasing  range.  The  slow  decay  of  the  transverse  comer 
relative  to  the  more  constant  radial/vertical  values  may  reflect  the 
importance  of  the  attenuation  model  in  interpreting  high  frequency  source 
comers  beyond  1200-1400  nu  The  lower  isotropic  comer  frequency  appears 
xmaffected  to  2200  m.  Mean  isotropic  comer  frequency  is  1.82  Hz  with  a 
multiplicative  error  term  of  1.12.  Mean  deviatoric  corner  frequency  is  3.28  Hz 
with  a  multiplicative  error  term  of  1.13. 

In  order  to  interpret  the  data  Mueller-Murphy  source  models  were  calculated 
for  detonations  of  1, 10/  and  100  kt  in  alluvium.  A  scaled  depth  of  122 
m/kt  was  used.  Velocity  and  density  structure  was  chosen  appropriate  for 
the  Coalora  site  and  consistent  with  values  used  to  interpret  the  scalar 
moments.  Envelop  functions  were  fit  to  the  Mueller-Murphy  source 
spectrum  to  replicate  analysis  applied  to  the  observations.  High  frequency 
decay  in  Mueller-Murphy  source  model  is  close  to  that  found  in  the  R/Z  data, 
(Table  la).  The  observed  isotropic  comer  frequency  of  1.82  falls  between 
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model  values  of  2.5  Hz  for  1  kt  and  1.7  Hz  for  10  kt,  consistent  with  the 
announced  yield  of  "less  than  20  kt." 


Spall  is  a  secondary  sou-ce  process  that  occurs  near  the  free  surface. 
Accelerometer  data  from  within  the  spall  zone  (Figure  6)  were  used  to 
develop  an  eqiiivalent  body  force  mode*,  for  this  process.  The  model 
represents  spall  as  a  set  of  vertical  point  forces  that  conserve  momenftim  (Day 
et  ah  1983)  which  can  be  shown  to  be  equivalent  to  the  standard  moment 
tensor  representation  (Day  and  McLaughlin,  1991).  Escape  velocity  of  near 
sixrface  material  is  constrained  by  integrals  of  the  accelerograms  while  total 
spall  mass  (mj)  is  estimated  from  the  data  (Figure  6).  Boimds  placed  on  the 
spall  process  were  discussed  in  the  earlier  observational  data  section.  Spall 
mass  estimates  (solid  and  dashed  lines.  Figure  6)  for  this  explosion  were  1.2- 
4.4x10^^  kg.  FoUcvving  Stump  (1985),  temporal  and  spatial  finiteness  of  the 
spall  zoxie  is  replicated  with  a  source  rise  time,  Tsr.  Initiation  time  of  the 
failure  process  across  the  test  bed  is  used  as  an  estimate  of  the  secondary 
source  lise  time.  Physically  this  rise  time  reflects  the  fact  that  spall  initiation 
does  not  occur  simultaneously  over  the  e*.  :re  test  bed  but  is  dependent  upon 
tejisile  strength  of  the  material,  geometry  of  the  explosion,  and  propagation 
veloci^  of  the  media.  In  addition  to  spall  rise  time  ,  an  average  spall  dwell 
time,  Ts,  is  needed  in  the  three  vertical  forces  that  make  up  the  total 
eqxuvalent  .seismic  source.  This  value  is  also  estimated  from  the  spall  zone 
data.  With  these;  parameters  the  three  vertical  forces  that  represent  (1)  spall 
initiation,  (2)  material  relaxation,  and  (3)  rejoin  become: 


fl  (t)  =  [(2*1^^  -  ^  mt  Vo  (H(t)  -  H(t  -  Tsr))] 

TiR  T|r  tIr 


(4a) 


f2(t)*-mTg 


+  (H(t-TsR)-H(t-Ts)j 
[[^  (6(t-TslP  15(t-T,)^  ^  lO(t-Ts)^ 
^  I  TIr  T^r  tIr 
\  •(H(t-Ts)-H(t-Ts-TsR)) 


(4b) 
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f3(t)  = 


(4c) 


i  tIr  T^r  tIr  I 

-H{t-Ts-TsR)) 


Spall  rise  time  for  Coalora  was  taken  as  300ms.  Dwell  time  varied  across  the 
test  bedy  a  value  of  300  ms  with  a  consistent  escape  velocity  of  1.47  m/s  was 
chosen.  The  approximate  equivalence  of  spall  rise  time  and  dwell  time  has 
been  observed  in  chemical  explosion  tests  (Stump,  1985).  The  smoothness  of 
the  time  function  deaeases  the  high  frequency  spall  signal.  A  lower  bound 
source  contribution  was  calculated  by  assuming  a  spall  mass  of  1.2x10^0  kg. 

The  resulting  source  time  function  and  its  spectrum  are  given  in  Figure  13. 
The  time  domain  results  emphasise  the  initial  downward  force  due  to 
material  failure  (1),  the  upward  relaxation  of  the  unspalled  body  (2),  and  the 
downward  force  associated  with  spaU  rejoin  (3).  The  model  predicts  a  peak 
force  of  1.1x10  dynes  and  a  time  duration  of  approximately  0.6  s.  The 
peaked  spectrum  results  in  little  contribution  from  spall  at  the  long  periods  or 
high  frequencies,  but  strong  effects  around  the  dominate  period  of  0.6  s.  This 
'dominant'  period  is  below  the  characteristic  comer  frequency  of  the 
explosion  model  (1. 7-2.5  Hz)  but  close  enough  in  frequency  to  bias  comer 
frequency  or  long  period  spectral  levels  if  it  is  strong.  The  peaked  nature  of 
the  secondary  source  could  in  some  cases  be  misinterpreted  in  terms  of 
explosion  overshoot.  The  tinie  integral  of  the  spall  force  model  gives  a 
maximum  value  of  7.2x10'^  dyne-cm  which  will  be  used  latter  to  compare  to 
the  time  derivative  of  the  moment  tensor  determined  by  inversion  of  the 
near-source  data. 

Full  use  of  transverse  waveforms  and  spectra  at  the  closest  ranges  (<2200m)  is 
dependent  upon  a  deviatoric  soxuce  model.  One  possible  representation  is 
that  attributed  to  Brune  (1970)  in  which  a  shear  stress  is  applied 
instantaneously  to  the  fault  surface  and  the  resulting  displacement  field 
determined.  Strictly  speaking  this  model  is  only  appropriate  for  earthquake 
sources.  One  could  imagine  that  this  shear  stress  is  a  result  of  the  initial 
compressive  waves  from  the  source  being  resolved  along  planes  of  weakness 
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Fjgiire  6:  Vertical  section  displaying  accelerograms  from  dovmhole  and 
siirface  gages  within  two  depths  of  burial  which  constrain  the  spall  zone. 

Solid  circles  represent  vertical  accelerometer  records  that  show  characteristic 
-1  g  dwell  indicative  of  spall.  Solid  and  heavy  dashed  lines  are  bounds  on  the 
depth  of  the  spall  zone. 
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Figure  13:  Spall  source  model  developed  from  observational  data  ^Sthin  spall 
zone.  Time  function  of  vertical  force  along  with  its  spectrum  are  given 
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MOMENT  RATE  TENSOT  (M) 


Figure  19:  Moment  rate  (Mjj)  and  moment  CM.j)  tensors  from  an  inversion  of 
the  Coalora  observational  data. 
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or  bedding  planes  in  the  geological  material  surrounding  the  explosion.  For 
such  a  deviatoric  source  model  the  shear  wave  (T)  spectrum  can  be 
interpreted  in  terms  of  moment,  stress  drop,  fault  radius,  and  average  slip. 
Application  of  this  model  to  the  Coalora  transverse  spectral  data  (Figure  8  and 
Table  la)  results  in  the  source  properties  given  in  Table  2.  Comer  frequency 
decreases  with  increasing  range  beyond  about  1  km  (Figure  12).  Since  the 
Brune  model  interpretation  assumes  the  spectral  corner  frequency  is 
characteristic  of  the  source,  data  within  1  km  of  the  explosion  were  used  in 
the  deviatoric  source  analysis.  Source  radii  estimated  from  these  data  range 
from  110-144  m  (124  m  mean,  8  85  standard  deviation  of  the  mean  -s.d.m). 
Stress  drops  varied  from  47-149  bars  (89  bars  mean,  34  s.d.m).  Average 
displacements  on  the  fault  were  18-56  cm  (33  cm  mean,  12  s.d.m).  Source 
displacements  calculated  from  the  radiated  wavefield  are  similar  in  size  to 
displacements  on  faults  and  bedding  planes  found  upon  re-entering  tunnels 
surroimding  nuclear  explosions  at  Rainier  Mesa  (Kennedy,  1984). 

Using  Coalora  shot  medium  properties,  equivalent  elastic  radii  for  a  1  and  10 
kt  explosion  were  calodated  to  be  133  m  and  202  m  respectively.  The  average 
deviatoric  source  radius  from  this  study,  124  m,  falls  inside  the  equivalent 
elastic  radii.  Bache  et  al  (1979)  studied  regional  records  (131-368  km)  from  the 
nuclear  explosions  Mighty  Epic  and  Diablo  Hawk  at  Rainier  Mesa,  NTS  to 
constrain  block  motions.  They  estimated  deviatoric  moments  for  the  two 
shots  to  be  1.3x10  dyne-cm  and  2.6x10  dyne-cm  respectively.  Since  comer 
frequencies  could  not  be  unambiguously  determined  by  Bache  et  al.  (1979), 
trade-offs  existed  for  the  range  of  possible  source  parameters.  Source  radii 
between  168-385  m  coupled  with  displacements  between  11-46  cm  and  stress 
drops  of  20-120  bars  bounded  the  possibilities.  Equivalent  elastic  radii  for  the 
two  events  were  estimated  at  160  m.  Deviatoric  source  radii  in  Bache’s 
analysis  appear  to  be  larger  than  elastic  radii  although  these  estimates  are  not 
directly  obtained  from  spectral  data.  Motions  along  faults  and  bedding  planes 
were  measured  following  the  Mighty  Epic  and  Diablo  Hawk  tests 
(Kcimedy,1984).  Fault  motions  from  Mighty  Epic  were  between  0.24  m  and 
1.68  m  at  ranges  of  74-124  m.  The  radius  which  encompasses  these 
displacements  is  slightly  less  than  the  elastic  radius  as  found  for  the  Coalora 
experiment.  These  displacements  represent  motions  on  isolated  planes 


Table  2:  Bnme  Model 


Station 

Range 

(m) 

Radius 

(m) 

Stress  Drop 
(bars) 

GM2 

274 

127 

76 

GM4 

365 

131 

87 

GM6 

389 

123 

100 

GM7 

612 

no 

66 

GM8 

614 

123 

49 

GM9 

612 

123 

149 

GMIO 

307 

123 

124 

GMll 

495 

117 

120 

SI 

850 

123 

70 

S3 

900 

144 

47 

S5 

1670 

149 

62 

S7 

1640 

144 

51 

S8 

1700 

144 

87 

Sll 

2100 

172 

54 

‘  Displacement 
(cm) 

29 

34 

37 
21 
18 
56 
46 
42 
26 
20 
28 
22 

38 
28 


in 


which  when  summed  would  yield  the  Coalora  near-source  transverse 
observations. 
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Wavef  onn  Modeling  and  Inversion 


Synthetic  Green's  functions  for  Yucca  Flats  around  the  Coalora  explosion 
were  developed  so  that  waveform  modeling  for  source  determinadon  could 
be  completed.  Geological  and  geophysical  structure  data  utilized  in  these 
calculations  were  presented  in  the  experimental  section.  The  velodty-depth 
model  is  given  in  Figure  14.  A  linear  velocity  gradient  in  the  top  0.3  km  is 
followed  by  a  jump  in  the  P  and  S  velocity  at  the  vitrified  tuffs.  A  slight  jump 
in  P  velocity  at  the  water  table  (500  m)  is  followed  by  a  positive  velocity 

gradient  to  the  Paleozoics  where  another  large  step  in  velocity  is  recorded.  A 
simple  attenuation  model  was  assumed  with  Qp  equal  to  100  and  equal  to 

44.  Although  important,  little  information  exists  for  constraining  shear  wave 
velocities.  Poisson's  ratio  was  taken  as  0.25.  Synthetics  were  calculated 
utilizing  the  extended  reflectivity  methodology  (Muller,  1985)  so  that  the 
effects  of  both  body  and  surface  waves  could  be  included.  Green's  functions 
for  a  pure  explosion  source  are  displayed  in  Figure  15  along  with  observations 
at  411/850  m  (solid  line-radial,  dotted  line- vertical). 


Rather  than  attempt  a  trial  and  error  waveform  matching  procedure,  Green’s 
functions  were  used  in  a  constrained  inversion  of  observational  data.  The 
first  trial  assumed  that  the  waveforms  could  be  modeled  with  just  an 
isotropic  source.  In  this  case  the  only  imknown  is  the  source  strength  and  its 
time /frequency  fimction: 

(5) 

where  are  the  radial  and  vertical  motions,  G^j^  ^  are  the  Green's 
functions  of  the  isotropic  source,  Mj^  are  the  absolute  source  strength,  and 

T(f)  is  the  source  spectral  shape.  In  the  constrained  inversions  the  source 
strength  (Mjck)  and  shape  (T(f))  are  determined  by  matrix  inversion  given  a 
set  of  Green's  functions  (Gnk,k)  and  observational  data  (Un).  This  procedure 
can  be  applied  to  single  or  multiple  source  observations.  Figure  16  is  the 
resulting  source  spectrum  and  time  functions  from  application  of  this 
procedure  to  one  set  of  vertical  and  radial  waveforms  (transverse  unfit 
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COALORA  VELOCITY  MODEL 


Figure  14:  Velocity  model  developed  for  Coalora  explosion  in  Yucca  Flats. 
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DEPTH  (km) 


Figure  15:  Isotropic  Green's  functions  for  velocity  model  in  Figure  16  and 
selected  observations  in  the  411-850  m  range.  Radial  components  are 
designated  by  solid  lines  while  vertical  components  are  dotted  lines. 
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ISOTROPIC  DECONVOLUTION 
0.411  km 


W0-3NAa 


Rgure  16:  Result  of  a  single  station  isotropic  sc 
was  at  41  Im.  Source  spectrum,  its  time  historj 
given. 


because  of  isotropic  assumption)  from  the  Coalora  explosion  at  a  range  of  411 
m.  Isotropic  moments  determined  by  this  procedure  for  all  single  station  data 
between  411  and  2100  m  are  given  in  Figure  17,  The  mean  isotropic  moment 
estimate  is  6.9x10  dyne-cm,  smaller  than  the  scalar  source  estimates  since 
Green's  functions  used  in  the  inversion  focus  more  energy  to  the  free  surface 
than  the  simple  propagation  path  corrections  used  in  scalar  source 
interpretation.  Although  the  constrained  isotropic  moments  have  a  factor  of 
4  scatter,  they  indicate  no  systematic  increase  or  decrease  in  moment  with 
range.  This  result  indicates  that  propagation  path  effects  within  this  range 
were  properly  taken  into  account  in  a  relative  way. 

These  constrained  inversions  may  be  biased  if  spall  or  deviatoric  source 
contributions  are  important.  Full  moment  tensor  inversions  were  conducted 
to  assess  these  secondary  source  contributions: 

In  this  case  all  six  elements  of  the  moment  tensor,  Mij,  are  determined  at  each 
frequency  value  given  the  matrix  of  Green's  functions,  Gni,j,  and  the 
observational  data,  Un,  in  the  frequency  domain.  The  observations  within  2 
km  of  tlie  explosion  working  point  were  used  in  these  inversions,  a  total  of 
eleven,  three-component  observations.  The  basis  for  this  data  choice  was 
discussed  earlier. 

Summary  of  one  source  inversion  is  given  in  Figure  18,  where  observed  and 
calculated  vertical,  radial,  and  transverse  velocity  records  are  reproduced. 

One  measure  of  the  adequacy  of  the  source  inversion  is  how  well  the 
resulting  source  predicts  the  observations.  Correlation  coefficients  for 
individual  records  range  from  0.71  to  0.94.  Transverse  motion?  are  as  well 
modeled  as  vertical  and  radial  components. 

A  second  criteria  in  judging  the  inversions  is  the  size  of  the  condition 
number  (ratio  of  largest  to  smallest  eigenvalue)  of  the  Gni,j  matrix  as  a 
fimcdon  of  frequency.  These  values  give  a  measure  of  resolution  of  model 
parameters  by  the  data.  The  condition  numbers  are  between  10  and  20  across 
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SINGLE  STATION  CONSTRAINED 
ISOTROPIC  MOMENT 


Figure  17:  Summary  of  peak  time  domain  moment  estimates  from  single 
station  isotropic  source  inversions  such  as  the  one  in  Figure  18. 


STATION  (km) 


c 


Figure  18:  Observed  and  calculated  seismograms  from  moment  tensor 
inversions  at  four  stations.  Correlation  coefficient  for  each  fit  is  noted  next  to 
waveform  pairs. 
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the  entire  frequency  band  for  the  array.  Experience  with  synthetic  inversions 
with  realistic  noise  indicate  that  these  values  are  more  than  acceptable. 

The  time  derivative  of  the  moment  tensor  and  its  integral  for  the  Coalora 
source  inversion  are  given  in  Figure  19.  The  source  is  dominated  by  diagonal 

elements  of  the  moment  tensor.  The  initial  pulse  is  symmetric  and  is 
followed  by  a  secondary  long  period  contribution,  largest  on  the  M22 

component,  the  vertical  dipole.  This  secondary  source  found  on  the  diagonal 
elements  of  the  moment  tensor  is  consistent  with  the  spall  forward  model 
discussed  earlier.  As  Day  and  McLaughlin  (1991)  show,  the  moment  tensor 
representation  of  the  spall  source  function  is  proportional  to  the  doubly 
integrated  body  force  representation  of  spall  scaled  by  a^/h,  where  o  is  the  P 
velocity  (1.9  km/s)  and  h  is  the  source  depth  (274iu;  assumed  in  the  moment 
tensor  inversions.  The  integral  of  the  body  force  spall  model  developed 
earlier  (7.15x10^'^  dyne-s)  can  be  compared  to  the  strength  of  the  secondary 
long  period  arrival  on  the  Mzz  component  of  the  moment  rate  tensor  in 
Figure  19.  The  source  strength  of  this  secondary  pulse  from  the  moment 
tensor  inversion  is  1.7x10^^  dyne-cm/s  which  compares  to  the  forward  spall 
model  estimate  of  9.4x10^^  dyne-om/s.  The  spall  forward  model  assumed  the 
lower  bound  of  spall  mass  (1.2  -  4.4x10^®  kg)so  with  an  intermediate  mass 
estimate  the  forward  model  from  the  spall  zone  data  and  the  inverse  model 
from  moment  tensor  inversion  would  be  in  very  dose  agreement  for  spall 
source  strength.  The  secondary  spall  pulse  in  Figxire  19  is  1.5  s  in  duration 
which  is  a  factor  of  two  longer  than  the  forward  spall  model.  This  inaeased 
duration  may  reflect  the  approximate  way  in  which  spall  spatial  finiteness  is 
taken  into  accouiit  in  the  forward  model. 

Off-diagonal  elements  of  the  moment  tensor  are  a  factor  of  3-10  smaller  than 
the  diagonal  elements  supporting  the  small  deviatoric  scalar  moments 
detennined  earlier.  Separation  of  the  deviatoric  components  of  the  moment 
tensor  is  complicated  by  the  apparent  spall  contribution.  No  attempt  was 
made  at  further  analysis  of  the  deviatoric  component. 
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The  peak  in  the  diagonal  elements  of  the  moment  tensor  is  7.9x10^®  dyne-an. 
This  moment  estimate  compares  with  the  mean  scalar  value  of  1.95x10^^  and 
the  average  constrained  value  of  6.9x1020 
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Implications 


The  Coalora  experiment  offered  the  opportunity  to  constrain  important 
source  processes  contributing  to  near-source  ground  motion  from  a  nuclear 
explosion.  Data  identified  important  contributions  from  the  isotropic  source, 
spall,  and  other  deviatoric  sources.  The  instrumental  array  included  gages 
from  within  the  spall  zone  out  to  ranges  comparable  to  the  scale  of  the  major 
geological  feature,  Yucca  Flats.  Detailed  site  characterization  data  from 
throughout  the  valley  were  available  for  constraining  propagation  path 
effects.  At  ranges  of  2  km  or  less  the  waveforms  were  simple  with  distinct 
differencet  between  the  radial /vertical  (R/Z)  and  transverse  (T)  motions.  As 
propagation  distance  increased,  all  three  components  of  motion  (RTZ) 
became  similar  in  appearance  in  both  the  time  and  frequency  domains.  At  5.1 
km  range  ground  motions  had  durations  in  excess  of  20  s  with  similar  comer 
frequencies  and  long  period  levels  for  all  (RTZ)  components.  Three 
dimensional  wave  propagation  calculations  for  the  site  indicate  that  the 
Yucca  Flats  two  and  three  dimensional  structure  become  important  at  ranges 
beyond  2  km.  Strong  similarity  in  spectra  and  long  duration  of  motion  at  5.1 
km  are  attributed  to  attenuation  and  scattering  effects  within  the  basin 
underlying  Yucca  Flats.  Effects  of  attenuation  on  the  deviatoric  source  can  be 
seen  as  a  decrease  in  apparent  comer  frequency  for  transverse  spectra  with 
range.  Source  studies  were  constrained  to  data  within  2100  m  of  the  source. 
Use  of  more  distant  data  could  residt  in  a  biased  source  estimate  at  Yucca 
Flats. 

The  isotropic  part  of  the  source  dominates  moment  tensor  inversions  with  a 
total  source  strength  near  8  x  10^^  dyne-cm.  Observed  isotropic  comer 
frequencies  (1.8  Hz)  fall  within  the  bounds  of  plausible  values  predicted  by  a 
Mueller-Murphy  source  model  (1  kt  -  2.5  Hz;  10  kt  - 1.7  Hz).  High  frequency 
decay  of  T  is  observed  also  in  support  of  the  Mueller-Murphy  source  model. 
Scalar  source  estimates  from  obseived  spectra  are  approximately  2.5  times 
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larger  than  complete  moment  tensor  estimates.  This  difference  is  attributable 
to  the  different  way  the  two  estimates  treat  downgoing  energy  and  surface 
waves  with  the  scalar  moments  being  biased  high.  The  multiplicative  error 
factor  for  the  scalar  isotropic  moment  was  1.36  representative  of  the  spatial 
variability  of  the  seismic  energy  at  these  near-sotirce  distances. 

The  deviatoric  source  has  higher  comer  frequencies  although  the  absolute 
moment  is  5-10  times  smaller  than  the  isotropic.  Simple  Brune  type  models 
lead  to  an  equivalent  source  radius  of  124  m  which  is  slightly  smaller  than 
the  elastic  radius,  a  moderate  stress  drop  of  89  bars,  and  an  average 
displacement  of  33  cm.  Thq  source  radius  and  average  slip  are  comparable  to 
observations  in  the  tuimels  of  Rainier  Mesa  following  explosior^s.  These 
near-source  observations  when  not  dominated  by  attenuation  or  scattering 
effects  give  a  consistent  characterization  of  the  deviatoric  source  contribution. 
In  the  case  of  data  from  Yucca  Flats  these  estimates  may  become  biased  when 
data  beyond  2  km  are  used  because  of  scattering  and  attenuation  effects. 

Forward  calculations  of  the  size  of  the  spall  contribution  are  comparable  to  a 
late  time,  longer  period  arrival  seen  in  the  moment  tensor.  This  source 
contribution  is  dominant  on  the  Mzz  component  and  agrees  with  similar 
results  for  chemical  explosions  (Stiimp,  1987).  The  spall  source  strength 
determined  from  forward  models  constrained  by  spall  zone  data  and  moment 
tensor  inversions  are  consistent  with  one  another.  The  duration  of  the  spall 
contribution  from  the  inversions  is  longer  than  that  determined  by  the 
forward  model  and  may  reflect  the  approximate  way  in  which  source 
fiiuteness  is  taken  into  account.  The  fact  that  the  spall  source  function  is 
delayed  in  time  and  longer  period  in  duration  than  the  spherically  symmetric 
explosion  function  means  that  it  could  lead  to  a  biased  estimate  of  the 
explosion  source  time  function  if  not  taken  into  account  Its  close  proximity 
in  time  to  the  explosion  time  function  means  that  it  could  be  interpreted  as 
overshoot  in  the  source  time  function.  Isotropic  source  deconvolutions  of 
the  Coalora  data  (Figure  16)  do  indicate  significant  overshoot  in  the  source 
time  function. 

Constrained  isotropic  source  inversions  offer  opportunity  for  single  station 
explosion  source  estimates.  Applying  this  procedure  to  data  spanning  a  wide 
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range  of  distances,  we  illustrated  that  there  is  no  trend  in  moment  (with 
respect  to  data  scatter)  over  the  range,  549  -  2200m.  This  observation  argues 
that  attenuation  has  been  properly  taken  into  account  over  these  distances. 

The  Coalora  experiment  .Jlustrates  the  importance  of  combining  motion 
measurements  from  within  the  source  region  (spall  zone)  and  outside  it  so 
that  a  complete  unambiguous  picture  of  the  equivalent  seismic  source  can  be 
interpreted.  Such  studies  could  provide  the  basis  for  developing  physically 
supported  scaling  relations  for  each  of  the  individual  source  contributions. 


124 


References 


Aki,  K-  and  Y.  B.  Tsai  (1972).  Mechamsms  of  Love-wave  excitation  by 
explosive  sources,  J.  Geophys.  Res.  77, 1452-1475. 

Bache,  T.  C,  W.  E  Farrell,  and  D.  G.  Lambert  (1979).  Block  motion  estimates 
from  seismological  observations  of  Mighty  Epic  and  Diablo  Hawk,  Defense 
Nuclear  Agency  Report  5007F,  Washington,  D.  C.  20305,  184  pages. 

Brune,  J.  N.  (1970).  Tectonic  stress  and  the  spectra  of  seismic  shear  waves  from 
earthquakes,  J.  Geophys.  Res.  75, 4997-5009. 

Burger,  R.  W.,  T,  Lay,  T.  C.  Wallace,  and  L.  J.  Burdick  (1986).  Evidence  of 
tectonic  release  in  long-period  S  waves  from  underground  nuclear  explosions 
at  the  Novaya  Zemlya  test  sites,  BuU.  Seism.  Soc  Am.  76,  733-755. 

Day,  S.  D.,  N.  Rimer,  J.  T.  Cherry  (1983).  Surface  waves  from  underground 
explosions  with  spall:  analysis  of  elastic  and  nonlinear  source  models.  Bull. 
Seism.  Soc  Am.  73, 247-264. 

Day,  S.  D.  and  K.  L.  McLaughlin  (1991).  Seismic  source  representations  for 
spall,  BnlL  Seism.  Soc.  Am.  81, 191-201. 

Eisler,  J.  and  F.  Chilton  (1964).  Spalling  of  the  earth’s  surface  caused  by 
underground  nuclear  explosions,  J.  Geophys.  Res.  69,  5285-5293. 

Eisler,  J.,  F.  Chilton,  and  F.  Sauer  (1966).  Multiple  subsurface  spalling  by 
underground  nuclear  explosions,  J.  Geophys.  Res.  71,  3923-3927. 

Ferguson,  J.  F.,  R.  N.  Felch,  C.  L.  V.  Aiken,  J.  S.  Oldow,  and  H.  Dockery  (1988). 
Models  of  the  Bouguer  Gravity  and  Geologic  Structure  at  Yucca  Flat,  Nevada, 
Geophysics,  53, 231-244. 

Gupta,  I.  N.  and  R.  E  Blandford  (1983).  A  mechanism  for  generation  of  short- 
period  transverse  motion  from  explosions,  BulL  Seism.  Soc  Am.  73,  571-591. 


125 


Haskell,  N.  A.  (1967),  Analytic  approximation  for  the  elastic  radiation  from 
contained  underground  explosion,  J.  Geophys.  Res.  72,  2583. 

Kennedy,  R.  P.  (1982).  Mighty  Epic/Diablo  Hawk  block  motion  program,  in 
DARPA/AFOSR  Symposium  on  the  Physics  of  Nonisotropic  Source  Effects 
from  Underground  Nuclear  Explosions,  DARPA-GSD-8203/AFOSR-NP-8201, 
Defense  Advanced  Research  Projects  Agency,  Arlington,  Va  22209. 

Kisslinger,  C.,  E.  J.  Mateker,  Jr.,  T.  V.  McEvilly  (1961).  SH  motion  from 
explosions  in  soil,  J.  Geophys.  Res.  66, 3487-3496. 

McLaughlin,  K.  L.,  L.  R.  Johnson,  and  T.  V.  McEvilly  (1983).  Two-dimensional 
array  measurements  of  near-source  ground  acceleration.  Bull.  Seism.  Soc. 

Am.  67,  349-375. 


Mueller,  R.  A.  and  J.  R.  Murphy  (1971).  Seismic  characteristics  of 
tmderground  nuclear  detonations.  Part  I,  Seismic  scaling  law  of  underground 
detonations.  Bull.  Seism.  Soc  Am.  61,  1675-1692. 

Ohta,  Y.  and  A.  Kubotera  (1968).  An  example  of  Love  waves  generated  by 
small  explosions,  Zisin  (J.  Seismol,  Soc  Japan),  Ser.  2, 21, 109-120, 

Patton,  H.  J.  (1988).  Source  models  of  the  Harzer  explosion  from  regional 
observations  of  fundamental-mode  and  higher  order  mode  surface  waves. 
Bull.  Seism.  Soc  Am.  78, 1133-1157. 

Patton,  H.  J,  (1990).  Characterization  of  spall  from  observed  strong  ground 
motions  on  Pahute  Mesa,  Bull.  Seism-  Soc  Am.  80,  1326-1345 

Perret,  W,  R.  and  R.  C.  Bass  (1974).  Free  field  ground  motion  induced  by 
underground  explosioris,  SAND74-0252,  Sandia  National  Laboratories, 
Albuquerque,  New  Mexico. 


126 


Schlittenhardt,  J.  (1991).  The  effects  of  spall  on  teleseismic  P-waves:  an 
investigation  with  theoretical  seismograms,  AGU  Monograph  on  Explosion 
Soxirce  Phenomenology  (in  press). 

Sobel,  P.  A.  (1978).  The  effect  of  spall  on  mj^  and  M^,  Teledyne  Geotech  report 
SDAC-TR-77-12,  Dallas,  Texas. 

Stump,  B.  W.  (1985).  Constraints  on  explosive  sources  with  spall  from  near¬ 
source  waveforms.  Bull.  Seism.  Soc  Am.  75,  361-377. 

Stump,  B.  W.  and  R.  E.  Reinke  (1987).  Experimental  seismology:  In  Situ 
source  experiments.  Bull.  Seism.  Soc  Am.  77,  1295-1311. 

Stump,  B.  W.  (1987).  Mathematical  representation  and  physical  interpretation 
of  a  contained  chemical  explosion  in  alluvium,  BiiU.  Seism.  Soc.  Am.  77, 
1312-1325. 

Taylor,  S.  R  and  G.  E.  Randall  (1989).  The  effects  of  spall  on  regional 
seismograms,  Geophys.  Res.  Lett.  16, 211-214. 

Vidale,  J.  E.  and  D.  V.  Helmberger  (1987).  Path  effects  in  strong  motion 
seismology,  in  Seismic  Strong  Motion  Synthetics  edited  by  Bruce  A.  Bolt, 
Academic  Press,  Orlando,  FI. 

Viecelli,  J.  A.  (1973).  Spallation  and  the  generation  of  surface  waves  by  an 
imderground  explosion,  J.  Geophys.  Res.  78, 2475-2487. 

von  Seggem,  D.  H.  and  R.  R.  Blandford  (1972).  Source  time  functions  and 
spectra  for  underground  nuclear  explosions,  Geophys.  J.  31, 83-97. 

Wallace,  T.  C.,  D.  V.  Helmberger,  and  G.  R  Engen  (1985).  Evidence  of  tectonic 
release  from  underground  nuclear  explosions  in  long-period  S  waves.  Bull. 
Seism.  Soc  Am.  75, 157-174. 


127 


( 


Wallace,  T.  C  and  J.  F.  Barker  (1985).  Modeling  near-field  observations  from 
nuclear  explosions  in  dry  tuff  and  alluvium,  AFGL-TR-85-0321,  Air  Force 
Geophysics  Laboratory,  Hanscom  AFB,  Ma  01731.  ADAI65227 

Werth,  G.  C.  and  R.  F.  Herbst  (1963).  Comparison  of  amplitudes  of  seismic 
waves  from  nuclear  explosions  in  four  mediums,  J.  Geophys,  Res.  68,  1463- 
1474. 

Wojdk,  G.  L.  and  D.  K.  Vaugn  (1984).  Large-scale  analysis  of  three- 
dimensional  seismic  waves,  WA-CA  R  8403,  Weidlinger  Associates,  Menlo 
Park,  California,  94025. 


128 


/ 


Near-Source  Characterization  of 
the  Seismic  Wavefield  Radiated  from  Quarry  Blasts 


Sharon  K.  Reamer* 
Klaus-G.  Hinzen** 
Brian  W.  Stump* 


March  1992 


*both  at  Department  of  Geological  Sciences,  Southern  Methodist 
University,  Dallas  Texas  75275 

•*Bundesanstait  fur  Geowissenschaftcn  und  Rohstoffe,  Stilleweg  2, 
3000  Haimover  51,  Federal  Republic  of  Germany 


keywords:  temporal  finiteness,  spatial  finiteness,  chemical  explosion, 

quarry  blast,  cylindrical  seismic  source,  linear  superpositioning 
techniques 


129 


i 


SUMMARY 


A  scries  of  controlled  seismic  experiments  perfonned  in  a  limestone  quarry  demonstrate 
the  utility  of  high  precision  electronic  detonators  in  studying  source  characteristics  of  multiple 
explosive  arrays.  At  near-source  ranges  (80-130  m),  where  source  dimensions  are  on  the  same 
order  as  source-receiver  distances,  the  influence  of  the  difference  in  travel  path  length  among 
individual  e.xplosions  on  the  seismograms  is  significant.  Focusing  of  the  seismic  energy  is 
observed  as  a  function  of  station  location  with  respect  to  the  source  array  and  is  attributed  to  the 
extended  source  length  (68-94  m)  and  firing  time  of  the  source  (380-544  ms). 

We  examine  two  methods  for  modeling  ripple-fired  explosions  at  near-source  ranges 
using  the  principles  of  superpositioning.  The  first  method  is  based  primarily  on  acquisition  of  an 
adequate  single  shot  signal  and  requires  well-constrained  shot  times.  Amplitude  variations 
which  result  from  travel  path  differences  are  not  modeled,  which  restricts  use  of  this  technique  for 
purposes  of  blast  vibration  reduction  to  larger  distances  (>2-3  source  dimensions)  where  the 
spatial  finiteness  effects  of  the  source  begin  to  diminish.  For  near-source  distances  (<  2  source 
dimensions),  we  successfully  model  multiple-source  seismograms  by  convolving  a  synthetic 
seismic  source  signal  for  a  single  explosion  with  individual  half-space  Green’s  functions 
calculated  for  each  explosion  in  the  array.  Our  single-source  model  for  a  cylindrically-shaped 
single  charge  (borehole  length  of  17.5  m  and  diameter  of  90  mm)  of  68  kg  consists  of  a  modified 
Mueller-Murphy  approximation  which  utilizes  source  parameter  estimates  taken  from  chemical 
explosion  study  results.  Model  parameters  include  a  final  cavity  radius  of  0.25  m  and  an  elastic 


radius  of  18  m.  The  final  model  is  obtained  by  convolving  the  simulated  single-source  time  series 


with  half-space  Green’s  fiinctions  calculated  for  several  source  depths  and  superposed  to 
approximate  the  spatial  extent  of  the  borehole.  The  relative  amplinide  and  phase  characteristics 
of  the  observed  single-source  signal  at  the  some  distance  (8U.6  m)  are  reproduced  by  this  model. 

Multiple-source  synthetic  seismograms  contain  individual  Green’s  functions  for  each 
source-receiver  distance  but  utilize  identical  sources  for  the  explosive  array.  Focusing  effects  are 
shown  to  be  due  to  the  effect  of  propagation  path  differences  between  individual  explosions  in 
agreement  with  the  results  of  Anderson  and  Stump  (1989)  in  simulating  multiple-source 
seismograms.  Good  fits  to  the  measured  production  shot  amplinide  spectra  are  obtained  with  the 
synthetic  spectra.  Spectral  peaks  are  well  matched  due  to  precision  of  the  firing  times  which  were 
controlled  by  electronic  detonators.  Our  example  of  delay  time  variances  for  32  ms  production 
shot  (Appendix)  argues  for  better  constraint  of  firing  times  for  controlled  seismic  experiments. 
Such  constraint  requires  a  1%  error  or  less  in  cap  firing  times  which  can  be  realized  by  the  use  of 
firing  systems  with  an  order  of  magnitude  increase  in  precision  compared  to  pyrotechnic 
detonators. 


INTRODUCTION 

The  present  study  focuses  on  the  generation  and  propagation  of  seismic  waves  from 
multiple  explosion  arrays.  We  analyze  seismic  data  collected  from  three  production  blasts  and  a 
single  explosion  in  a  limestone  quarry  to  assess  the  influence  of  an  explosive  source  extended  in 
space  and  time  on  the  seismic  wavefield.  We  also  explore  the  efficacy  of  linear  superpositioning 
techniques  and  equivalent  source  models  for  chemical  explosions. 
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Previously,  seismological  studies  of  ripple-fired  explosions  were  limited  and  dispers  i 
somewhat  unevenly  in  the  literamre.  Interest  in  recent  years  has  increased  mainly  due  to  the 
necessity  of  discriminating  large  commercial  blasts  from  small  nuclear  explosions.  Several 
recent  studies  have  addressed  the  discrimination  problem  at  regional  distances  (Baumgardt  and 
Seglcr,  1988;  Hedlin  et  al,  1989;  and  Su  ct  al,  1991),  although  in  most  cases  with  limited 
knowledge  of  or  constraints  on  the  seismic  source  used  in  deriving  the  proposed  discriminants. 
Smith  (1989)  examined  spectral  characteristics  of  multiple-delay,  multiple-row  explosions  and 
identified  a  significant  tradeoff  between  attenuation  due  to  scattering  and  source  characteristics 
at  regional  distances. 

Anderson  and  Stump  (1989)  present  results  of  modeling  near-source  seismograms 
recorded  from  multiple-row,  multiple-delay  production  shots  in  a  granite  quarry  in 
Massachusetts.  Individual  Green’s  functions  for  each  source-receiver  distance  in  the  array  were 
calculated  and  convolved  with  a  .simulated  single-source.  Extended  source  time  duration 
(temporal  finiteness)  and  propagation  path  differences  due  to  extended  source  spatial  dimensions 
(spatial  finiteness)  were  well  reproduced  despite  a  limited  knowledge  of  the  actual  firing  times  or 
constraints  for  the  single-source  model.  Hinzen  (1988)  used  superpositioning  with  weighted 
amplitudes  to  model  a  series  of  five  explosions  fired  in  a  single  row.  The  recorded  seismograms 
from  a  single  explosion  fired  separately  from  the  row  shot  but  located  next  to  the  multiple  array 
were  used  for  the  supeipositioning.  By  using  an  in  situ  single  shot,  the  seismic  response  of  the 
individual  sources  and  extended  time  duration  of  the  multiple  explosions  were  well  modeled, 
although  amplitude  dififerences  due  propagation  path  differences  from  the  extended  source 
dimensions  were  not  reproduced.  It  was  thought  that  amplitude  variations  could  be  attributed  to 
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coupling  differences  between  individual  charges.  Some  studies  have  indicated  that  nonlinear, 
dynamic  source  effects  may  be  observable  in  the  elastic  region  (Minster  and  Day,  1986). 
However,  experimental  confirmation  of  supeipositioning  for  two  charges  fired  simultaneously 
with  spatial  separations  of  different  lengths  (Stump  and  Reinke,  1988)  a^ues  against  dynami'; 
interaction  of  the  sources  being  observable  in  far-fleld  seismic  data. 

As  noted  by  Hinzen  (1988),  knowledge  of  the  exact  firing  times  is  essential  to  any  study  of 
ripple-fired  explosions.  The  data  in  the  present  study  are  unique  in  that  firing  times  are  controlled 
with  electronic  detonators.  The  electronic  firing  system  reduces  scatter  in  firing  times  and  allows 
precise  measurement  of  actual  detonation  times.  This  provides  us  the  opportunity  to  test  and 
compare  observed  production  shot  seismograms  to  time  series  calculated  from  two 
superpositioning  methods:  (1)  empirical  superpositioning  with  an  in  sini  single  shot,  and 
(2)  linear  superpositioning  with  calculated  source  and  travel  path  functions. 

Additionally,  the  single  shot  data  provides  constraints  on  our  single-source  model.  As 
discussed  by  Henmarui  et  al  (1989),  there  is  a  need  for  refinement  of  source  models  for 
small-yield  chemical  explosions.  We  examine  the  Mueller-Murphy  (Mueller  and  Murphy, 
1971)  source  model  with  parameter  estimation  adjusted  for  our  quarry  e.xplosions.  Scaling  of 
source  parameters  for  small  chemical  explosions  at  depth  using  the  Mueller-Murphy  model  is 
explored  in  detail  by  Grant  (1988)  using  both  a  forward  modeling  approach  and  moment  tensor 
inversions.  We  pay  particular  attention  in  the  synthetic  calculations  to  modeling  the  single 
explosion  in  the  forward  case  and  thereby  extend  the  work  of  Anderson  and  Stump  (1989)  in 
modeling  the  temporal  and  spatial  finiteness  of  the  multiple  explosion  seismograms. 
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DESCRIPTION  OF  EXPERIMENTS 


SlaslJkslim 

The  experiments  consisted  of  three  multiple  explosions  (production  shots)  and  one  in  situ 
single  shot  fired  in  a  limestone  quarry  in  nonhcm  Italy.  Figure  1  shows  the  plan  of  the  quarry 
bench  and  shot  locations  for  three  production  shots  (SV13,  SV15  and  SV16)  and  the  single  shot 
(SVI4)  examined  in  the  present  study.  Each  production  shot  was  fired  in  a  single  row  at 
neighboring  parts  of  the  1 5  m  highwail.  The  firing  direction  for  each  production  sliot  is  indicated 
by  the  arrows  in  front  of  the  highwail  in  Figure  1 .  The  single  shot,  SVI4  was  also  located  on  the 
quarry  bench  and  recorded  at  similar  distances  to  the  production  .shot  experiments.  Shot  SV13 
produced  14.4  kilotons  (kt)  of  material  with  a  total  charge  weight  of  1.26  tons  (t)  of  explosives, 
SVI5  produced  16.0  kt  with  1.35  f,  and  SVI6  produced  20.7  kt  with  1.73  t. 

The  single  shot,  SVI4,  consisted  of  68  kg  of  explosives  including  11  leg  of  high  energy 
explosives  in  the  bottom  of  the  hole  and  37  kg  of  ammonium  nitrate  fuel  oil  (ANFO).  Tlie 
c-xplosives  were  initiated  from  the  top  of  the  hole.  The  charges  were  stemmed  with  3  m  of  drill 
cuttings,  fhe  borehole  was  drilled  at  an  inclination  of  30’  to  the  venical  with  a  borehole  lengdi  of 
17.5  m,  diameter  of  90  mm.  and  burden  of  3  m.  Loading  of  the  boreholes  for  the  production  shots 
was  identical  to  the  single  s.hot  for  most  but  not  ail  cases.  For  example,  16  out  of  20  holes  forSVI3 
were  loaded  as  dc.scribed  for  SVI4,  Tlic  other  four  holes  contained  high  energy  mixtures  in  the 
middle  of  the  charge  column  instead  of  at  the  bottom.  However,  total  charge  weights  per  hole 
were  kept  con.stant  at  63  kg  for  oil  production  .shon  consi<iered  here.  A  crosscut  view  of  the 
typical  borehole  is  shown  in  Figure  2  in  addition  to  a  pl.ui  view  of  SV15  and  The  20 
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boreholes  for  SVI5  had  spacing  and  burden  of  4  m  and  3  m,  respectively,  for  a  total  array  length  of 
76  ra.  A  constant  delay  time  of  20  ms  was  used  which  gives  a  total  desired  firing  time  of  380  ms. 
Actual  shot  tin.es  for  S  VI3  were  not  recorded.  For  the  26-hoIe  array,  SVId,  borehole  spacing  and 
burden  were  the  same  as  for  S  VI3  except  for  holes  22-26  where  the  quarry  wall  turns.  Delay 
intervals  for  SVI6  alternated  between  18  and  27  ms  for  a  total  desired  firing  time  of  540  ms 
{Figure  2).  Actual  shot  times  recorded  for  SVI6  (Table  1)  vary  from  desired  firing  times  by  values 
ranging  from  0.02%  to  1.5%.  Shot  times  for  SV13  (not  shown  in  Figure  2),  with  18  holes  (burden 
and  spacing  at  4  and  3  m,  respectively)  and  a  constant  32  ms  delay  (544  ms  total  desired  firing 
time),  were  also  recorded  with  a  maximum  deviation  from  desired  times  of  1.4%.  This  small 
variation  between  desired  and  actual  firing  times  was  achieved  by  the  use  of  electronic 
detonators. 

Seisniic  Station  Locatiuns  and  Data  Aenuisition 

Station  locations  for  each  production  shot  formed  a  stretched  semicircle  around  the  row  of 
explosions.  The  distance  between  each  of  these  stations  and  the  closest  borehole  was  80  m. 
Stations  A  through  E,  as  shown  in  Figure  1,  are  considered  in  this  study.  Station  A  was  located 
80  m  behind  the  last  borehole  (hole  1 8)  of  S  V13  (measured  perpendicular  to  the  explosive  array), 
while  stations  B  and  C  were  positioned  80  m  behind  the  first  and  last  boreholes  of  SVI3, 
respectively  .  Station  D  was  located  80  m  behind  the  first  borehole  of  SVI6,  and  station  E  was 
positioned  80.6  m  distance  from  the  single  shot,  SVI4.  Each  station  consisted  of  a 
three-component  4.5  Hz  geophone  with  62%  critical  damping.  All  stations  were  calibwted  with 
a  shaking  table.  Signals  were  digitally  recorded  with  a  central  recording  system  ufilizir<g  64 
channels.  The  data  were  sampled  at  the  rate  of  1 5.625K  samples  per  second  with  a  dynamic  range 
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of  12  bits  including  polarity.  Length  of  the  recorded  signal  was  1  s  for  the  production  shot  and 
s  for  the  single  shot  e-ipeviment.  The  recording  system  is  described  in  detail  by  Hinzen  (198b). 

The  production  shots  were  initiated  by  a  computer  driven  electronic  firing  system.  TTiis 
system  offers  a  maximum  of  60  time  steps  of  which  18, 20,  and  26  have  been  used  in  S  VI3,  S  VI5, 
and  SVI6.  The  detonators  contain  an  integrated  circuit  which  interacts  with  the  blast  computer. 
Each  detonator  is  programmed  at  the  factory  for  one  of  the  possible  time  steps.  In  the  field,  the 
time  length  of  the  delay  step  is  preset  at  the  blast  computer  to  values  between  1  and  100  ms.  This 
preset  delay  time  multiplied  by  the  time  step  of  the  detonator  gives  the  absolute  firing  time.  By 
dropping  time  steps  in  the  row  of  the  detonators,  unequal  delay  time  intervals,  as  were  used  for 
SVI6,  can  be  realized.  A  detailed  description  of  the  firing  system  is  given  by  Hinzen  et  al  (1987). 

EXPERIMENTAL  DATA 

Figure  3  shows  velocity  seismograms  from  stations  A-D  for  the  SVI3,  SV15,  and  SVT6. 
The  X  and  y  components  arc  the  two  perpendicular  horizontal  components  where  the  x  direction  is 
parallel  to  the  explosive  array  (Figure  1).  The  z  component  is  vertical.  All  seismograms  are 
normalized  to  maximum  trace  amplitudes,  indicated  in  mm/s  at  the  end  of  each  seismogram.  The 
firing  direction  for  SVI5  and  SVI6  recedes  from  statioas  B  and  D  and  for  S  VT3  and  SVI5  comes 
closer  to  stations  A  and  C  (Figure  1). 

The  time  duration  of  the  seismic  record  is  proportional  to  the  total  duration  of  the  sources. 
Duration  of  the  seismograms  is  approximately  1.2-1. 3  times  the  total  firing  time  at  stations  B  and 
D  and  1.3-1. 4  for  stations  A  and  C.  Except  for  the  initial  compressive  wave  pulse,  individual 
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signal  characteristics  of  a  single  source  are  obscured  in  the  multiple  source  seismograms.  The 
most  notable  feature  is  the  constructive  and  destructive  interference  of  seismic  energy  as  seen  in 
the  seismograms  due  to  the  multiple  explosions.  Seismic  energy  peaks  in  the  beginning  of  the 
seismogram  for  stations  B  and  D  and  at  the  end  for  stations  A  and  C.  This  is  clear  from  the 
cumulative  seismic  trace  energy,  shown  below  each  set  of  station  seismograms  (Figure  3).  The 
dashed  lines  connect  the  first  arrival  of  seismic  energy  with  the  point  where  98%  of  the  total 
energy  per  trace  is  reached.  Cumulative  trace  energy  at  time  I-At  is  measured  by  the  quantity 

^eum* 

Ecuni(I)  = 

i  *1  j»il 

where  Vy  is  ground  velocity  in  mm/s.  The  index  i  runs  from  1  to  the  total  number  of  samples  and  j 
represents  the  two  horizontal  and  one  vertied  components  of  the  recording.  If  a  constant  amount 
of  seismic  energy  arrives  as  a  function  of  time,  the  cumulative  energy  would  follow  the  dashed 
lines.  The  convex  curvature  of  the  cumulative  energy  for  stations  B  and  D  illustrates  that  the 
energy  arriving  from  SV15  and  SVI6,  respectively,  is  greater  at  the  beginning  and  decreases 
gradually  as  the  shot  detonation  moves  farther  away  from  the  station  with  time.  The  energy  curve 
forms  a  concave  shape  for  stations  A  and  C,  where  the  arriving  energy  increases  with  time  as  the 
detonating  charges  progress  closer  to  the  stations. 

SUPERPOSITIONING 

Empirical  Model  Calculations 

Our  first  approach  is  to  reproduce  the  production  shot  seismograms  in  the  time  domain 
using  the  hybrid  modeling  technique  of  Hinzen  (1988).  This  empirical  modeling  method  consists 
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of  linear  superpositioning  of  an  observed  single-source  seismogram  with  appropriate  delay  times 
to  reproduce  a  multiple  shot  seismogram.  This  method  is  frequently  utilized  to  reduce  ground 
vibrations  due  to  blasting  (Crenwelege,  1991;  Hinzen  and  Reamer,  1991).  Propagation  effects 
and  source  propenies  are  contained  in  the  recorded  signal  from  a  single  explosion  used  to 
synthesize  production  blasts  at  the  same  or  a  nearby  location.  The  optimal  synthetic  production 
blast  signal  is  calculated  by  varying  firing  time  sequences  until  seismic  vibrations  are  adequately 
reduced  (usually  in  the  frequency  range  between  5  and  20  Hz).  The  method  requires  precise 
knowledge  of  the  firing  times.  Linear  interaction  of  seismic  waves  from  neighboring  explosions 
and  identical  seismic  source  functions  arc  assumed.  Since  the  single  shot  contains  only  the 
propagation  effects  along  a  single  source-receiver  path,  differences  in  propagation  due  to  the 
actual  path  are  not  modeled  directly;  however,  differences  in  path  length  are  included  by  the 
addition  of  an  assumed  compressional  wave  travel  time.  In  the  time  domain,  this  can  be 
expressed  as  a  series  of  convolutions  given  by: 


n 


V(x,t)  =  S(x’,t’)  ®  G(x,t;x’,t’)  ®  Z  a;  5(t-ii)  (2) 


i«l 

where,  V(x,t) 

=> 

particle  velocity 

G(x,t ;  .x’,t’)®S(x’,t‘) 

= 

representation  of  measured  seismic  signal 

X 

= 

spatial  coordinates  of  station 

x’ 

= 

spatial  coordinates  of  source 

ai 

= 

weighted  amplitude.s 

ti 

delay  and  travel  times,  and 

n 

= 

total  number  of  explosions. 

For  this  experiment,  the  measured  singl 

e-source  seismogram  at  station  E  (inset  in  Figure  4) 

is  convolved  with  a  time-delayed  sequence  of  unit  impulses.  Each  time  delay  includes  firing  time 
of  the  shot  and  travel  time  (with  an  assumed  P-wave  velocity  of  4.5  km/s)  from  source  to  receiver. 
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The  impulse  series  simulates  the  temporal  extent  of  the  multiple-source  array.  Results  of  the 
linear  superpositioning  for  station  B  of  SV15  and  station  A  of  SVI3  (Figure  4)  show  that  total 
source  duration  is  well  matched  by  the  linearly  superposed  seismograms.  However,  simple  linear 
superpositioning,  as  shown  in  the  upper  traces  in  Figure  4,  does  not  adequately  reproduce  the 
amplitudes  observed  in  the  measured  data.  Variations  in  amplitude  can  be  modeled  direcdy  only 
if  the  firing  times  are  well  known,  as  is  the  case  for  this  experiment.  Individual  amplitude  values 
from  the  observations  are  modeled  by  weighting  the  amplitudes  of  the  linearly-superposed 
seismograms  at  the  appropriate  delay  times  to  agree  with  the  amplitudes  of  the  observed  multiple 
source  event  (weighting  factors  are  determined  directly  from  the  measured  production  shot 
seismograms).  The  resulting  velocity  seismograms  (bottom  traces  in  Figure  4)  match  the 
observational  data  quite  well  demonstrating  the  effectiveness  of  this  superpositioning  technique 
if  shot  times  are  well  constrained.  This  technique  can  be  a  useful  tool  for  blast  vibration 
reductions;  however,  the  method  does  not  provide  us  with  any  physical  understanding  of  the 
mechanism  for  interference  effects  observed  in  the  production  shot  seismograms.  As  noted  by 
others  (Smith,  1989;  Anderson  and  Stump,  1989),  spatial  finiteness  of  the  production  shot  array 
and  temporal  finiteness  due  to  time  delay  blasting  combine  to  produce  a  seismogram  extended  in 
time  (proportional  to  total  source  firing  time)  and  exhibiting  characteristics  of  destructive  and 
constructive  wave  interference.  The  rest  of  this  paper  will  be  devoted  to  understanding  the  nature 
of  these  interference  effects,  as  observed  at  near-source  distances. 

Theoretical  Model  Calculations 

We  attribute  the  observed  amplitude  variations  in  the  production  shot  seismograms  to 
differences  in  propagation  path  among  individual  explosions  in  the  array.  Our  primary 


139 


motivation  is  to  quantify  the  contributions  to  the  multiple  shot  seismograms  from  the  seismic 
source  and  propagation  path.  We  therefore  simulate  the  production  shot  signals  using  a  forward 
modeling  approach.  We  calculate  individual  Green’s  functions  for  each  shot  of  the  array 
following  the  method  of  Anderson  and  Stump  ( 1 989).  Each  Green’s  function  is  convolved  with  a 
single-source  time  function  and  linearly  superposed  with  the  appropriate  delay  times  to 
reproduce  the  multiple  shot  seismograms.  This  modeling  approach  is  represented  by: 


V(x,t) 
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i=l 


S(Xi,t’)  ®  G(x,t;xi.t’)  0 
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n 


seismic  particle  velocity 
individual  Green’s  functions 
single-source  model 
spatial  coordinates  of  station 
spatial  coordinates  of  sources 
delay  times,  and 
the  total  number  of  explosions. 


(3) 


As  in  the  empirical  modeling  approach,  each  individual  shot  of  the  multiple  array  is 


assumed  to  have  an  identical  source  function.  An  elastic  half-space  with  a  compressional  wave 
velocity  f  4.5  km/s  is  assumed  for  the  Green’s  function  as  a  first-order  approximation  to  the 
quarry  stnicture  estimated  from  average  P-wave  arrival  times  for  me  measured  data.  Density  is 
assumed  tc  be  2600  kg/m^,  a  reasonable  value  for  the  competent  limestone  in  this  area.  Individual 
Green’s  functions  arc  calculated  using  the  method  of  Johnson  (1974)  for  each  source-receiver 
distance  b.  he  array  configuration.  Figure  5  shows  radid  and  vertical  component  Green’s 
functions  convolved  with  the  single-source  model  calculated  for  station  D.  The  details  of  the 


single  source  model  calculations  will  be  fully  discussed  in  the  next  section.  Here  we  only  wish  to 
illustrate  that  the  change  in  Green’s  function  from  the  first  to  last  explosion  is  not  dramatic 


140 


(Figure  5),  but  is  significant  when  considered  in  addition  to  the  delay  times  between  shots.  The 
source-receiver  distance  increases  from  80  to  126.7  m  between  the  first  and  last  explosion  in  the 
SVI6  array  (26  holes).  P  and  SV-Rayleigh  waves  are  clearly  separated  at  these  distances.  The 
amplitude  ratio  of  P  to  SV-Rayleigh  peaks  is  0.22  and  0.18  for  the  vertical  and  0.78  and  0.72  for 
the  radial  component  at  the  shortest  and  largest  distance,  respectively. 

The  single  source  is  estimated  from  a  parametric  seismic  model  for  explosions  originally 
derived  for  simulating  the  seismic  response  of  contained  nuclear  explosions  (Mueller  and 
Murphy,  1971).  The  model  proceeds  at  the  "boundary”  between  elastic  and  non-elastic  response 
in  a  continuum  to  the  application  of  a  sphericaily-symmetric  pressure  function.  The  utility  of  the 
model,  whose  basic  principles  were  first  derived  by  Sharpe  (1942),  resides  in  the 
parameterization  of  the  solution  in  terms  of  "measurable  characteristics”  and  including  the  effect 
of  source  depth  (Mueller  and  Murphy,  1971).  Stump  (1985)  and  Grant  (1988)  have  successfully 
used  this  method  to  model  small  chemical  explosions  (114,6  and  2.3  kg)  at  depth  in  alluvium. 
Details  of  the  mathematical  derivations  of  the  semi-analytical  model  will  not  be  reproduced  here 
as  they  are  well  documented  in  the  literature  (Mueller  and  Murphy,  1971;  Stump,  1985;  Grant, 
1988).  Calculation  of  the  model  in  the  frequency  domain  is  given  by  the  following  equation, 

p(C0) 

^  (4) 

4\l  (cOg^  -  j3co^  +  icOgCo) 


where,  \}/((d)  =  far-field  displacement  potential  (units  of  volume) 
rei  =  elastic  radius 
c  =  compressional  wave  velocity 
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H  =  shear  modulus 

O)  =  angular  frequency 

C0(,  =  theoretical  comer  frequency 

P  =  (X+2li)/4|i,  (k,  p.  are  Lame's  constants) 

p(co)  =  frequency  domain  pressure  function. 


The  pressure  function  as  specified  as  in  Stump  (1985)  and  in  Mueller  and  Murphy  (1971) 
represents  a  step  pulse  in  time  with  an  e.xponential  decay  and  is  parameterized  by  peak  pressure, 
static  pressure,  and  a  decay  constant  (proponional  to  tOo).  Estimation  of  the  model  parameters 
elastic  radius,  peak  pressure  and  static  pressure  bears  further  discussion.  Grant  (1988)  narrows 
the  acceptable  parameter  range  by  imposing  constraints  on  the  material  propenies,  panicularly 
the  shear  wave  velocity.  However,  in  the  present  study,  we  have  a  limited  knowledge  of  the 
material  properties  in  our  area,  and  we  must  constrain  the  range  of  reasonable  values  for  the 
pressure  function  and  elastic  radius  parameters  using  a  different  approach. 


Static  pressures  arc  calculated  according  to  the  proportionality  relation  from  Mueller  and 
Murphy  (1971)  for  competent  rock. 


P,  =  4/3n 


(5) 


where,  P, 

=  static  pressure 

It 

=  shear  modulus 

fc 

=  final  cavity  radius 

fel 

=  elastic  radius. 

Explicit  in  this  calculation  are  the  parameters  for  final  cavity  radius  and  clastic  radius.  We 
attempt  to  determine  reasonable  constraints  for  these  values  in  our  model  calculations  by  using 
results  from  other  chemical  explosion  studies. 
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First  we  consider  that  cavity  radius  for  chemical  explosions  scales  differently  than  for 
nuclear  explosions.  Grant  (1988)  scaled  cavity  radius  for  2.27  kg  explosions  by  altering  the 
medium-dependent  proportionality  constant  (k)  relating  cavity  radius  (rc)  in  meters  to  yield  (Y) 
in  kilotons  and  depth  of  burial  (h)  in  meters  given  by  : 

yO.29 

fc  =  k  •  (6) 

Using  reasonable  values  for  compressional  velocity  (4.5  km/s),  density  (2600  kg/m^),  and  shear 
modulus  (7.8x10*°  N/m^)  for  the  limestone  in  the  quarry,  and  a  proportionality  constant,  k,  of  25 
derived  using  a  relation  for  nuclear  explosions  (Mueller  and  Murphy,  1971),  equation  (4)  predicts 
cavity  radius  values  for  the  68  kg  single  explosion  in  our  study  of  between  1 .4  to  1 ,2  m  (between  3 
and  17.5  m  depth).  However,  based  on  both  theoretical  and  experimental  studies  of  chemical 
explosions  in  various  media,  most  predictions  for  cavity  radius  are  from  2  to  1 2  times  the  original 
borehole  radius  (Qiiappetta  et  al,  1987)  producing  a  range  of  values  between  0. 1  and  0.6  meters 
for  an  initial  cavity  radius  of  0.05  m. 

Elastic  radius  is  often  but  not  always  defined  as  the  ’Tsoundary"  around  the  explosion 
beyond  which  the  material  behaves  elastically  (Mueller  and  Murphy,  1971;  Sharpe,  1942).  For 
chemical  explosions,  it  can  also  be  related  to  the  region  of  fracture  growth  and  damage  to  the 
material  surrounding  the  borehole  (Kutter  and  Fairhurst,  1971)  and  is  a  more  transitional 
measure.  There  exist  many  published  relationships  between  explosive  yield  and  clastic 
radius/damage  zones  (e.g.,  Atchison  and  Toumay,  1959;  D'Andrea  et  al,  1970;  Siskind  and 
Fumanti,  1974)  for  chemical  explosions  in  hard  rock  from  studies  which  sample  small  explosive 
yields  (0.002  to  4  kg).  The  range  of  elastic  radii  predicted  by  these  relations  falls  between  0.5  and 


5.3  m  for  the  single  68  kg  charge  in  this  study.  Another  frequently  used  estimator  for  elastic 
radius  is  the  transition  zone  observ'ed  in  amplitude  decay  rates  as  a  function  of  scaled  range. 
Using  the  decay  rate  curves  for  nuclear  explosions  in  hard  rock,  the  range  of  scaled  elastic  radius 
predicted  for  68  kg  is  between  4  and  8  m  (Ferret  and  Bass,  1975).  Based  on  peak  velocity 
amplitude  decay  rates  for  chemical  explosions  (mainly  production  blasts)  in  hard  rock,  an  initial 
estimate  for  the  elastic  radius  is  between  16-18  m  (Ambraseys  and  Hendron,  1968). 
Interestingly,  this  measure  coincides  with  the  length  of  the  borehole  (17.5  m). 

One  possible  interpretation  of  the  variability  in  predicted  elastic  radius  values  is  that  a 
cylindrically-shaped  charge  scales  differently  than  a  charge  of  spherical  shape,  elastic  radius 
being  constrained  by  the  length  of  the  borehole  in  the  cylindrical  case.  This  is  supported 
experimentally  by  the  different  scaling  relationships  for  spherical  and  cylindrical  charges 
necessary  for  determining  the  radius  of  damage  for  explosions  in  plexiglass  and  hard  rock  (Kutter 
and  Fairhurst,  1971).  Many  Bureau  of  Mines  studies  were  also  carried  out  with  cylindrical 
charges.  Nicholls  and  Duvall  (1966)  present  results  which  support  a  volumetric  rather  than 
charge  weight  relationship  to  the  damage  zone.  Siskind  and  Fumanti  (1974)  also  allude  to  a  ’’rule 
of  thumb”  relation  of  between  1  and  2  times  the  length  of  the  charge  column. 

When  values  for  cavity  radius  obtained  from  the  nuclear  explosion  relations  are  substituted 
in  equation  (5),  static  pressures  are  1-2  orders  of  magnitude  greater  than  peak  pressures 
calculated  from  the  overburden  relation  given  by  Mueller  and  Murphy  (1971), 

Pp  =  1.5  pgh.  (7) 

> 
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where,  Pp  =  peak  pressure 

p  =  density 

g  =  acceleration  due  to  gravity 

h  =  overburden  depth. 

Reasonable  static  pressure  values  are  obtained  from  equation  (5)  using  cavity  radius  values 
(0. 1-0.6  m)  based  on  the  chemical  explosion  study  results.  We  depend  on  the  results  of  the  model 
calculations  to  better  constrain  the  elastic  radius  values. 

A  suite  of  theoretical  sources  in  the  frequency  domain  are  first  calculated  using 
equation  (4),  differentiated,  Fourier  synthesized,  and  finally  convolved  with  the  calculated 
Green’s  function  at  80  m  to  produce  velocity  seismograms.  Reasonable  approximations  to  the 
observed  single  source  seismogram  at  station  E  are  obtained  by  varying  cavity  radius  from  0.05  to 
0.3  m  ( 1-6  times  the  borehole  radius)  and  elastic  radius  from  10-20  m.  As  a  first  approximation, 
the  source  functions  are  calculated  assuming  an  overburden  depth  of  8  m  and  a  Green’s  function 
depth  of  8  m  (near  the  middle  of  the  borehole). 

Two  criteria  used  in  adjusting  the  model  parameters  arc  peak  amplitude,  as  measured  in  the 
time  domain,  and  comer  frequency,  as  measured  from  the  velocity  spectra.  While  not  the  only 
criteria  used  in  the  model  selection,  these  two  measures  are  important  because  they  can  be 
directly  estimated  from  and  compared  between  the  observed  and  calculated  seismograms. 
Comer  frequency  is  inversely  proportional  to  clastic  radius  (Mueller  and  Murphy,  1971)  and  is 
used  as  a  constraint  on  this  model  parameter.  We  estimate  comer  frequency  directly  from  the  data 
by  finding  the  change  in  slope  at  the  peak  spectral  amplitude  level.  This  is  consistent  with  a 
parametric  spectral  model  for  explosions  which  postulates  a  static  long-period  displacement 
amplitude  below  the  comer  frequency  and  a  high  frequency  spectral  decay  above  the  comer 
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(Reamer  and  Stump,  1992).  We  feel  that  the  comer  frequencies  estimated  in  this  way  ax- 
consistent  to  within  ±2  Hz. 

The  second  criterium,  peak  amplitude,  is  a  more  delicate  parameter  to  adjust.  As  shown  in 
equation  5,  as  cavity  radius  increases  and  elastic  radius  remains  constant,  static  pressure  increases 
(and  overshoot  decreases).  Peak  amplitude  also  increases  as  the  source  time  function  receives  a 
bigger  contribution  from  the  static  pressure  while  peak  pressure  remains  the  same.  In  the 
frequency  domain,  the  effect  is  an  increased  long-period  amplitude  level  (propononal  to  the 
increase  in  static  pressure).  Returning  again  to  equation  5,  if  elastic  radius  increases  and  cavity 
radius  remains  constant,  static  pressure  decreases  (and  overshoot  increase^).  In  this  case,  the 
source  time  function  peak  amplitudes  are  increased  due  to  the  increased  overshoot.  In  our  first 
series  of  tests,  we  experiment  with  combinations  of  elastic  and  cavity  radius  values  to  determine 
their  effect  on  peak  amplitude  and  comer  frequency  in  the  resulting  time  series  (Table  2).  Wr 
conclude  that  differences  between  calculated  and  observed  peak  amplitudes  and  comer 
frequencies  are  too  large,  even  though  synthetic  waveforms  and  spectral  shapes  are  in  fair 
agreement  with  the  observations. 

For  the  second  series  of  tests,  refinement  of  model  parameters  is  necessary.  If  elastic  radius 
and  cavity  radius  are  kept  constant  (constant  static  pressure),  peak  amplitude  can  still  be  adjusted 
by  varying  the  overburden  depth,  thereby  adjusting  the  peak  pressure  (equation  7).  Since  the 
boreholes  were  stemmed  with  drill  cuttings  from  the  surface  to  3  m  depth  and  the  charges  were 
initiated  from  the  top  of  the  charge  column,  an  overburden  of  3  m  seems  a  logical  choice. 
However,  after  some  experimentation,  we  found  that  an  overburden  depth  of  5  m  produced  better 


peak  amplitude  results.  A  physical  justification  for  this  result  may  be  found  by  reviewing  the 
explosive  process.  A  detonation  or  shock  front  fonns  behind  the  explosive  reaction  as  it 
propagates  through  the  borehole,  and  a  (measurable)  detonation  pressure  builds  directly  behind 
the  detonation  front.  For  cylindrical  explosions,  peak  borehole  pressure  can  be  expressed  as  a 
percentage  of  the  explosive  detonation  pressure  although  variations  can  be  quite  high  for  ANFO 
(30-70%)  (Chiappetta  et  aJ,  1987).  A  panial  explanation  for  the  difference  in  explosive  versus 
borehole  pressure  may  be  that  some  of  the  explosive  (shock)  energy  is  channeled  into  fracturing 
the  quarry  face,  thereby  decreasing  the  effective  source  strength. 

The  cylindrical  shape  of  the  charge  column  (14.5  m  length)  also  creates  a  dilemma  for  the 
Green’s  function  calculations  which  require  a  single  source  depth.  After  testing  several  source 
depths,  we  obtain  the  best  re.suits  for  the  single-shot  seismograms  by  linearly  superposing  the 
Green's  functions  at  different  depths.  Individual  Green's  functions  are  calculated  at  depths  of  3, 
5,7,9, 11, 13, 15  and  17  m.  A  downhole  detonation  velocity  of  4800  m/s  is  assumed,  giving  a  total 
detonation  time  of  3  ms  from  the  top  to  the  bottom  of  the  hole.  In  our  second  set  of  tests,  selected 
source  models  are  convolved  with  the  Green’s  functions  and  linearly  superposed  with  the 
appropriate  downhole  detonation  time.  The  seismic  moment  of  the  synthetic  source  is 
’’distributed”  over  the  source  depths  by  normalizing  the  linearly-superposed  seismogram  by  the 
number  of  source  depths  (8)  used  in  the  superpositioning.  Tabic  3  gives  parameters  and  results  of 
the  second  set  of  models  using  the  superposed  Green's  functions.  Variations  in  peak  amplitude 
and  comer  frequency  are  small  compared  to  the  first  model  tests.  Based  on  comparison  to  the 
peak  amplitude  and  comer  frequency  estimated  from  the  observed  seismogram  at  station  E  (top 
traces  in  Figure  6),  as  well  as  some  other  criteria  (discussion  follows),  we  chose  a  "best  fit”  model 
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from  this  scries  of  tests  (indicated  in  bold  in  Table  3)  with  a  final  cavity  radius  of  0.25  m  (5.0  times 


the  original  borehole  radius  of  0.05  m)  and  an  elastic  radius  of  18  m  (Figure  6). 

The  radial  and  vertical  component  observed  seismograms  at  station  E  (Figure  6)  consist 
mainly  of  a  compressional  wave  pulse  (15  ms  period)  followed  by  the  SV-Rayleigh  wave 
(20-25  ms  dominant  period).  The  major  phases  and  relative  peak  amplitudes  between  radial  and 
venical  components  arc  well  modeled  by  the  synthetics  (bottom  traces  in  Figure  6).  Ratio  of 
vertical  to  radial  peak  amplitudes  is  1.05  for  the  observed  seismograms  and  1.09  for  the 
synthetics.  For  the  radial  observed  seismograms,  the  peak  amplitude  of  the  Rayleigh  wave  is 
slightly  larger  than  the  P-wave  pulse  due  to  the  extended  source  depth.  We  are  able  to  match  this 
relative  amplitude  only  by  using  the  superposed  (downhole)  Green’s  functions. 

Misfit  between  the  calculated  and  observed  venical  seismograms  is  seen  in  the  first 
downw.'jd  swing  of  the  P  wave  cycle  followed  by  a  phase  slightly  smaller  in  amplitude.  On  the 
hodogram  plots  (shown  at  right  in  Figure  6),  it  is  clear  that  this  phase  is  oriented  at  approximately 
90®  t  j  the  primary  P  wave  pulse  and  probably  represents  the  initial  SV  pul.se  arriving  from  the 
source.  The  phase  can  be  seen  on  the  particle  motion  plot  of  the  observed  data  to  be  slightly  open, 
p<;ssibly  due  to  the  arrival  of  the  free  face  tensile  reflection  (at  approximately  2  ms)  which  would 
have  the  effect  of  deepening  the  first  downward  swing  of  the  P  wave  phase  on  the  n  iial 
component.  In  a  true  cylindrical  source  model,  both  P  and  S V  energy  are  generated  at  the  source 
(Heclan,  1953);  however,  we  are  only  approximating  a  cylindrical  source  with  the  superposed, 
half  space  Green’s  functions  and  so  are  unable  to  adequately  reproduce  either  the  relative 
amplitudes  between  the  P  and  SV  phases  or  the  quarry  face  reflection. 
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Choice  of  our  final  model  is  also  constrained  by  the  spectral  response  of  the  measured  data. 
However,  the  low  frequency  spectral  response  (below  the  comer  frequency)  is  primarily 
controlled  by  the  cavity  radius  (Grant,  1988),  and  our  final  parameter  adjustments  were  made  to 
match  both  the  low  frequency  response  and  the  comer  frequency  (controlled  by  elastic  radius)  of 
the  observations.  High  frequency  response  is  automatically  controlled  by  our  choice  of  source 
function  and  is  the  same  for  all  model  calculations.  Figure  7  compares  amplitude  spectra  of  the 
observed  and  synthetic  seismograms.  The  spectral  shape  and  ampilmdes  are  in  good  agreement 
from  the  comer  frequencies  at  40  Hz  (vertical)  and  45  Hz  (radial)  to  the  high  frequency  limit  at 
400  Hz.  The  synthetic  seismograms  from  our  final  model  give  the  best  fit  to  the  observations  at 
and  below  the  comer  frequency  for  both  components.  The  radial  component  of  the  synthetic  data 
contains  less  low  frequency  energy  between  6  and  35  Hz  than  the  observed  data  (maximum  factor 
of  four  amplitude  difference  at  28  Hz).  Between  10  and  30  Hz,  the  vertical  observed  amplitudes 
are  almost  a  factor  of  two  times  higher  than  the  synthetics. 

MuJtjpler-Source  SynMietic  Seismouram  Calculations 

For  the  explosion  source  model  utilized  in  the  present  study,  we  calculate  only  two 
components  of  motion,  radial  and  vertical.  The  observed  data  were  recorded  as  three-component 
seismograms  and  contain  a  considerable  amount  of  transverse  energy.  In  addition,  the  production 
shot  seismograms,  which  are  oriented  as  shown  in  Figure  2,  cannot  be  rotated  into  a  true  "radial” 
or  "transverse"  direction  due  to  the  spatial  extent  of  the  source  array.  Therefore,  it  should  be  noted 
for  the  supeipositioning  results  that  follow,  a  basic  difference  exists  in  that  for  each  Green’s 
function+source  convolution,  the  orientation  is  radial  with  respect  to  the  source  whereas  in  the 
observations,  only  the  borehole  closest  to  the  station  is  recorded  in  the  radial  sense.  We  refer  to 
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the  y  component  of  the  observed  production  shot  seismograms  as  the  ’’equivalent”  of  the  radiai 


component. 

For  shot  SVI6,  we  calculated  eight  Green’s  functions  at  the  different  depths  (the  same  as  for 
the  single  shot  model)  for  each  source-receiver  distance  in  the  array  (a  total  of  208  Green’s 
functions).  For  each  component  at  each  distance,  the  Green’s  functions  are  first  convolved  with 
our  source  model  (parameters  given  in  Table  2)  and  then  linearly  superposed  with  the  assumed 
downhole  detonation  time.  Then  the  synthetic  seismograms  for  each  source-receiver  location 
are  superposed  with  the  actual  delay  times  recorded  for  SVI6.  The  results  for  the  radiai  and 
vertical  components  of  station  D  are  shown  in  Figure  10.  Seismograms  are  all  plotted  relative  to 
the  same  scale.  The  first  0. 1 6  s  of  the  observed  radial  seismogram  contains  1 8  ms  delay  intervals 
complicated  by  the  fact  that  the  firing  order  for  the  first  four  shots  was  not  linear  sequential 
(Figure  2).  This  first  0.16  s  wave  packet  is  distinct  in  both  the  observed  and  synthetic 
seismograms  from  the  next  0.32  s  where  the  delay  intervals  alternate  between  18  and  27  ms.  The 
last  five  shots  (on  the  seismograms  from  0.50-0.64  s)  also  alternate  between  1 8  and  27  ms  delays, 
but  the  holes  were  located  5  m  perpendicular  distance  (y  direction)  fanher  from  the  station  in 
addition  to  the  increasing  horizontal  separation.  The  additional  distance  manifests  in  the 
observed  and  synthetic  seismograms  as  a  reduced-amplitude  wave  packet. 

For  comparison,  at  the  top  of  Figure  8  we  show  the  superposed,  measured  single  shot 
(SMSS)  signal  for  SVI6  using  the  seismogram  from  SVI4  measured  at  station  E.  Note  that 
although  the  first  0. 16  s  wave  packet  can  be  observed  in  the  SMSS,  the  focusing  effect  seen  in  the 
observed  and  synthetic  seismogram.^  is  r.ot  reproduced  here.  For  the  vertical  component,  the 
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attenuation  of  amplitudes  in  the  observed  seismogram  is  more  pronounced  than  in  the  synthetic 
although  the  interference  effects  of  all  three  wave  packets  are  qualitatively  reproduced.  Again, 
focusing  effects  are  not  reproduced  by  the  SMSS.  The  ratio  of  peak  radial  to  peak  vertical 
amplitudes  for  the  measured  data  is  1.1,  for  the  synthetics  1.4,  and  1.4  for  the  SMSS  signal. 

Superpositioning  results  are  shown  for  the  vertical  component  of  stations  B  and  C  for  shot 
SVI5  in  Figure  9.  Since  we  did  not  have  a  single  shot  signal  located  on  the  same  quarry  wall  as 
SVI5,  we  used  the  S VI4  single  shot  source  model  for  the  synthetics  and  the  SMSS.  We  obtained 
the  best  results  for  the  synthetics  by  convolving  the  single-source  model  with  the  Green’s 
functions  at  1 1  m  depth  for  each  source-receiver  distance,  rather  than  using  the  Green’s  functions 
at  all  depths.  This  result  may  be  attributable  to  propagation  path  diffeie.icos  at  the  different 
benches;  it  can  be  seen  from  Figure  1  that  the  two  shots  are  located  on  benches  perpendicular  to 
each  other.  Slight  differences  in  the  source  model  due  to  inhomogeneities  ir  the  source  region 
differences  in  source  coupling  also  cannot  be  ruled  out.  Amplitudes  are  plotted  to  the  same  scale 
for  all  seismograms  except  the  station  C  measured  seismogram.  Due  to  higher  noise  levels,  this 
signal  was  low-pass  filtered  to  200  Hz.  Also,  we  utilize  only  desired  firing  times  with  constant 
delay  intervals  of  20  ms  for  the  superpositioning,  since  actual  firing  times  were  not  recorded. 
Firing  direction  is  away  from  station  B  and  towards  station  C;  the  focusing  effect  of  firing 
direction  is  seen  in  both  the  measured  and  synthetic  seismograms.  Peak  amplitudes  for  station  B 
agree  well  for  the  vertical  component  although,  as  for  station  D  of  S  VI6,  relative  amplitudes  are 
not  accurately  reproduced.  At  station  C,  the  seismogram  shape  of  the  synthetic  signal  compares 
well  to  the  measured  data  even  though  absolute  amplitudes  do  not  agree. 
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Better  results  could  be  obtained  with  amore  detailed  propagation  path  model;  however,  for 
a  first-order  approximation,  the  half-space  assumption  works  surprisingly  well.  Reasonable 
comparisons  are  obtained  for  SVIS  and  SVI6  for  two  main  reasons:  ( 1 )  good  similarity  of  the 
seismic  source  function  between  the  measured  single  shot  and  the  individual  explosions  of  the 
multiple  array,  and  (2)  accuracy  of  the  delay  times  which,  in  the  case  of  SVIS,  allows  us  to  use  the 
desired  firing  times  in  our  superpositioning  model.  We  illustrate  this  last  point  by  examining  a 
spectral  domain  equivalent  representation  of  our  multiple-source  seismograms  given  by  (after 
Blair,  1988): 


n  °  1/2 

A(f)  =  S(f)  { [  I  aj  cos(2t;  tj  f )]  ^  +  [  aj  sin(2jt  tj  f )]  ^  )  (g) 

1=1  i=»i 


where,  A(0  = 

S(f)  = 

Oi  = 

fi  = 
n  = 
f  = 


amplitude  spectral  values 
source  funaion  (constant  for  each  shot) 
amplitude  weighting  term  (different  for  each  shot) 
delay  time 

the  total  number  of  explosions 
frequency. 


We  introduce  this  representation  to  illustrate  the  temporal  variations  as  observed  in  the 
velocity  spectra  due  to  rime  delay  blasting.  Peaks  in  the  spectra  correspond  to  the  harmonic 
frequencies  associated  with  each  delay  and  occur  in  multiples  up  to  the  Nyquist  frequency. 
Spectral  smearing  occurs  due  to  slight  changes  in  the  source-receiver  distances  (propagation 
path  effect)  which  also  enhances  high-frequency  damping  of  the  amplitudes  (Smith,  1989). 
Measured  and  superposed  spectra  of  the  venical  component  of  station  B  (SVIS)  and  the  radial 
component  of  station  D  (SVI6)  (Figure  10)  exhibit  the  scalloping  pattern  associated  with  the 
delay  time  firing  panem.  For  SVIS,  the  first  peak  at  SO  Hz  corresponding  to  the  constant  20  ms 
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delay  interval  is  quite  clear  for  the  measured,  synthetic  and  SMSS  spectra.  Peaks  at  100, 150,200, 
250  and  300  Hz  are  all  observable  in  the  measured  data  although  the  peak  at  150  Hz  is  split.  The 
harmonics  due  to  rippIe-flring  dominate  the  production  shot  spectra  so  that  source  spectral 
characteristics  of  the  individual  explosions  cannot  be  easily  discerned. 

For  SVI6,  the  firing  pattern  is  complicated  because  of  alternating  13  and  27  ms  delays.  In 
addition,  the  first  four  shots  were  not  fired  in  the  row  order  causing  slight  shifts  in  the 
source-receiver  offsets  in  addition  to  the  18  ms  delays.  What  is  seen  in  both  the  measured  and 
superposed  spectra  is  a  small  first  peak  at  22  Hz  corresponding  to  a  45  ms  delay  time,  which  is  just 
the  sum  of  the  two  alternating  delays.  A  second  smaller  peak  at  34  Hz  corresponds  to  the  first 
27  ms  harmonic.  The  second  harmonic  of  the  combined  delay  times  is  at  44  Hz  followed  by  the 
first  harmonic  of  the  1 8  ms  delay  at  56  Hz.  The  higher  order  harmonics  of  the  combined  delay  can 
be  picked  to  about  200  Hz  in  the  synthetic  and  SMSS  although  the  peaks  are  smeared  in  the 
measured  data.  Higher  order  harmonics  of  the  1 8  and  27  ms  delays,  although  present,  cannot  be 
discerned  in  the  measured  spectra. 


DISCUSSION 

We  have  successfully  modeled  the  temporal  and  spatial  variations  due  to  ripple-fired 
explosions  as  measured  in  seismic  data  acquired  at  near-source  ranges.  The  single-source  model 
was  calculated  using  the  analytic  Mueller-Murphy  (1971)  explosion  model  with  parameters 
constrained  by  results  from  chemical  explosion  studies.  We  tested  models  for  a  single  cylindrical 
charge  of  68  kg  ANFO  in  limestone  with  borehole  length  of  17.5  m  and  diameter  of  90  mm 
measured  at  80  m  distance  by  varying  cavity  radius  between  0.05  and  0.3  m  and  elastic  radius 
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between  10  and  20  m.  TTie  best-fit  model  is  obtained  with  a  final  cavity  radius  of  0.25  m,  ar. 
elastic  radius  of  1 8  m,  and  an  overburden  depth  for  the  peak  pressure  function  of  5  m.  In  order  to 
match  relative  amplitudes  between  the  P  and  Rayleigh  wave  phases  of  the  observed  data,  it  is 
necessary  to  linearly  superpose  the  seismic  response  of  Green's  functions  calculated  at  different 
depths  to  simulate  the  spatial  extent  (14.5  m)  of  the  cylindrical  source. 

As  shown  here,  temporal  and  spatial  finiteness  effects  of  ripple-fired  blasting  at  80  m 
distance  are  well  modeled  with  linear  supeipositioning  using  a  calculated  source  and  Green’s 
funaions  when  constraints  can  be  imposed  on  both  the  propagation  medium  and  the  seismic 
source  models.  However,  even  in  the  near-source  region  (80  m),  the  method  requires  accurate 
firing  times  ( 1  %  or  less  error).  It  has  yet  to  be  conclusively  shown  that  ripple-fue  effects  of  blasts 
with  small  delay  times  (<  40  ms)  can  provide  a  consistent,  singular  discrimination  critcrium  for 
quarry  blasts  at  regional  distances,  especially  when  one  considers  that  the  effects  of  blasting  cap 
firing  time  inaccuracies  (see  Appendix)  are  convolved  with  attenuation  and  local  site  effects. 
More  work  is  needed  with  explosions  employing  both  large  and  small  delay  times  in  controlled 
experimental  settings  to  gain  a  better  understanding  of  the  interaction  of  propagation  effects  and 
delay  time  variations  at  near-regional  and  regional  distances. 
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APPENDIX 


CAP  SCATTER  AND  SEISMIC  SPECTRAL  RESPONSE 

The  combination  of  attenuation  mechanisms  and  firing-time  scatter  can  obscure  the 
spectral  harmonics  associated  with  ripple  firing.  Smith  (1989)  observed  no  spectral  signature  at 
regional  distances  for  small-delay  (17  ms)  overburden  blasts.  Spectral  harmonics  at  regional 
distances  can  also  originate  due  to  propagation  effects,  either  at  the  source  or  the  receiver 
(Suteau-Henson  and  Bache,  1988;  Hedlin  et  al,  1989).  We  want  to  isolate  one  element  affecting 
the  spectral  modulation,  cap  scatter  effects.  Cap  scatter  is  here  defined  as  a  percentage  of  the 
desired  firing  time  and  represents  the  deviation  between  desired  and  actual  firing  time.  Firing 
time  deviation  has  the  same  effect  as  spatial  finiteness  of  the  explosive  anay,  namely,  smearing  of 
the  spectral  modulation  pattern.  Modification  of  equation  (8)  to  include  a  random  variation  in 
delay  times  results  in: 

**  **1/2 

A(f)  =  S(f)  { [Z  aj  cos(2jt(ti  +  n)  f  )]  2  +  [I  aj  sin(2K(ti  +  rO  f )]  2 )  (Al) 

i=l  i=l 

where,  q  =  noise  value  for  each  delay  time. 

Scatter  in  the  delay  times  is  not  a  simple  noise  term  added  to  the  seismic  data  but  is  inherent  to  the 
time  series  and  corresponding  amplitude  spectra.  Pattern  recognition  methods  such  as 
homomorphic  deconvolution  cannot,  therefore,  ’’see”  the  regular  scalloping  pattern  produced  by 
ripple-fired  blasting  if  the  delay  time  variations  become  too  large. 

To  quantify  the  influence  of  cap  scatter  on  the  seismic  spectra,  we  have  calculated  a  unit 
amplitude  impulse  series  representing  the  desired  firing  sequence  observed  at  station  B  for  S  VIS 
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(32  ms  delay).  Geological  structure  in  this  case  is  assumed  to  be  homogeneous  and  isotropic  with 
travel  times  between  individual  charges  and  the  seismic  station  added  to  the  delay  times  with  an 
assumed  compressional  wave  speed  of  4.5  km/s.  The  impulse  response  scries  and  corresponding 
amplitude  spectra  for  S  VI3  is  shown  in  the  bottom  trace  of  Figure  A1 .  The  spectra  of  the  impulse 
series  for  the  exact  times  show  the  first  peak  at  3 1.25  Hz,  corresponding  to  the  inverse  of  the  delay 
times  (32  ms)  and  higher  harmonics  occurring  at  multiples  of  the  fundamental  frequency. 

Cap  scatter  is  added  to  the  desired  firing  times  using  normally  distributed  random  values 
with  maximum  variance  levels  of  1%,  2%,  4%  and  6%  of  the  desired  firing  times  (Figure  Al). 
The  increasing  amount  of  cap  scatter  disturbs  the  scalloping  pattern  of  the  spectra.  With  1  %  delay 
time  variations,  the  maximum  possible  deviation  between  desired  and  actual  firing  time  is  +0.32 
ms  for  time  step  1  and  +5.44  ms  for  time  step  17.  Even  at  this  low  noise  level,  the  regular 
scalloping  structure  of  the  spectra  is  smeared  for  frequencies  higher  than  100  Hz.  With  2% 
maximum  cap  scatter,  only  the  first  peak  at  3  ..25  can  be  clearly  correlated  with  the  spectrum  for 
exact  firing  times.  For  noise  levels  of  ''<  /o  and  6%,  the  ripple  strjcture  is  completely  destroyed. 
Troughs  appear  at  frequencies  where  peaks  are  observed  in  the  spectra  for  exact  firing  times.  For 
example,  alarge  spectral  peak  appears  at  the  1%,  2%,  and  4%  noise  level  at  about  20  Hz.  As  these 
results  indicate,  even  a  small  value  of  1%  noise  added  to  the  firing  times  disturbs  the  spectral 
modulation  from  the  single-row  production  shots.  Usually  cap  manufacturers  do  not  give 
detailed  information  about  the  scatter  of  cap  firing  times.  Some  studies  have  shown  that  average 
cap  scatter  values  of  4%  are  fairly  representative  (Blair,  1988)  and  variances  as  high  as  20%  of  the 
desired  firing  times  can  be  realized  (Reamer  et  al,  1989). 
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Table  1.  Comparison  of  Desired  and  Actual  Firing  Times  for  SVI6 


Time 

Steo 

Desired  Firing 
Time  (ms) 

Actual  Firing 
Time  (ms) 

Percent 

.  Difference 

0 

0 

0.0 

0.0 

1 

18 

17.73 

1.5 

2 

36 

35.84 

0.4 

3 

54 

54.02 

0.0 

4 

72 

71.87 

0.2 

5 

90 

89.98 

0.0 

6 

117 

116.93 

0.0 

7 

135 

134.85 

0.1 

8 

162 

162.30 

0.2 

9 

180 

180.45 

0.2 

10 

207 

207.55 

0.3 

U 

225 

225.73 

0.3 

12 

252 

252.86 

0.3 

13 

270 

270.91 

0.3 

14 

297 

297.79 

0.3 

15 

315 

315.58 

0.2 

Id 

342 

342.78 

0.2 

17 

360 

360.62 

0.2 

18 

387 

387.90 

0.2 

19 

405 

406.40 

0.3 

20 

432 

432.90 

0.2 

21 

450 

451.20 

0.3 

22 

477 

476.91 

0.0 

23 

495 

495.08 

0.0 

24 

522 

517.50 

0.9 

25 

540 

536.05 

0.7 
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Table  2.  Model  parameters  and  results: 

8  m  overburden  depth  and  8  m  Green’s  function  depth 
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Table  3.  Model  parameters  and  results: 

5  m  overburden  depth**  and  superposed  Green’s  functions 
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FIGURE  CAPTIONS 


Figure  1.  Three  production  shots,  SVI3  and  SVI5  and  SVI6,  were  fired  at  neighboring  parts  of 
the  15  m  highwall  of  the  quarry.  Single  shot  SVI4  (open  circle)  was  located  dose  to  the  21st 
borehole  of  the  S  VI6  array.  There  were  18, 20  and  26  individual  e.xplosions  for  shots  SVI3,  SVI5 
and  SVI6,  respectively.  Seismic  stations  A-E  as  shown  here  were  located  on  the  quarry  bench. 
Station  A  corresponds  to  SVI3;  stations  B  and  C  correspond  to  SV15;  station  D  corresponds  to 
SVI6;  and  station  E  corresponds  to  SVI4. 

Figure  2.  The  plan  views  for  SVI5  and  SVI6  are  shown  in  the  upperpart  of  the  figure.  A  constant 
delay  time  of  20  ms  was  chosen  for  SVI5.  Delay  times  for  SVI6  alternated  between  18  and  27 
ms.  Spacing  and  burden  were  4  and  5  m,  respectively  for  both  SVI5  and  SVI6  for  total  array 
lengths  of  76  and  94  m.  Burden  and  spacing  for  the  last  five  boreholes  of  the  SVI6  shots  arc 
irregular  as  the  array  turns  a  "comer”  of  the  highwall.  The  crosscut  view  at  the  bottom  of  the 
figure  shows  the  typical  hole  loading,  borehole  inclination  of  30  degrees  to  the  vertical,  borehole 
length  of  17.5  m,  and  diameter  of  90  mm  for  the  produaion  shots  and  the  single  shot,  SVI4. 

Figure  3.  Measured  three-component  velocity  seismograms  are  plotted  at  the  top  of  the  figure 
for  firing  direction  away  from  the  station  (Station  B  of  SVI5  and  station  D  of  SVI6). 
Seismograms  from  stations  C  (SVI5)  and  station  A  (SVD)  are  plotted  in  the  lower  part  of  the 
figure  and  represent  the  signal  recorded  from  firing  in  a  direction  toward  the  seismic  station.  The 
two  perpendicular  horizontal  components  are  x  and  y,  where  x  is  parallel  to  the  explosive  array; 
z  is  the  vertical  component.  Seismograms  are  normalized  to  the  maximum,  given  in  mm/s  at  the 
end  of  each  trace.  Curves  plotted  below  the  seismograms  arc  the  normalized  cumulative  seismic 
trace  energy.  The  straight  dashed  lines  connect  the  first  arrivals  with  the  98%  energy  level. 

Figure  4.  Vertical  component  velocity  seismograms  of  shot  SVI3  (right)  and  SVI5  (left)  at 
station  A  and  B,  respectively.  The  upper  traces  (a)  consist  of  the  linearly  superposed  single-shot 
signal  from  shot  SVI4  with  a  time-delayed  series  of  unit  amplitudes.  The  middle  traces  (b)  arc 
the  recorded  seismograms  from  the  production  shots.  The  lower  traces  (c)  arc  the  weighted 
amplitude  results.  Weighting  factors  are  estimated  directly  from  measured  production  shot 
amplitudes  at  the  cumulative  delay  times.  The  measured  single-shot  signal  is  shown  as  inset  (d) 
at  the  same  time  scale  as  the  production  shots. 

Figure  5.  Half  space  Green’s  funaions  from  eight  depths  convolved  with  the  single  source  time 
function  and  linearly  superposed  with  a  total  of  3  ms  downhole  detonation  time  were  calculated 
for  travel  paths  from  SV16  to  station  D.  The  distance  increases  from  80  to  126.7  m  from  the  first 
to  the  26th  hole.  Every  second  time  series  is  plotted  here.  The  sections  on  the  left  and  right  are 
the  vertical  and  radial  components,  respectively  with  amplitudes  plotted  respective  to  the 
maximum  in  mm/s  (scale  shown). 
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P  ’’ire  6.  Measured  (a)  and  calculated  (b)  seismograms  of  the  single  shot,  SVI4  at  80.6  m 
distance  (station  E).  The  radial  and  venical  seismograms  on  the  left  and  right  side,  respectively, 
ar^  plotted  to  the  same  scale,  given  at  left,  and  the  maximum  amplitudes  are  plotted  above  each 
tra.e.  Hodograms  in  the  radial/vertical  plane  are  shown  at  the  right  of  the  figure. 

FigUi  e  7.  Velocity  spectra  of  the  measured  (at  station  E)  (thick  lines)  and  calculated  (thin  lines) 
singb  shot  signal,  SVI4.  Spectral  amplitudes  in  mm/s/Hz  are  shown.  Both  measured  and 
cala  lated  signals  are  filtered  with  a  low-pass,  2  pole  Butterworth  filter  at  400  Hz. 

Figure  8.  Superpositioning  results  are  shown  for  the  vertical  and  radial  component  seismograms 
from  SVI6  at  station  D.  The  upper  traces  (a)  are  the  superposed  measured  single  shot  (SMSS) 
signals  (SVI4).  The  middle  traces  (b)  are  the  observed  production  shot  seismograms.  The  lower 
traces  (c)  are  obtained  by  convolving  the  source  time  function  with  the  Green’s  functions 
calculated  for  each  shot-receiver  distance  and  linearly  superposing  the  resulting  seismograms 
with  alternating  18  and  27  ms  delay  times. 

Figure  9.  Superpositioning  results  are  shown  for  the  vertical  component  seismograms  from 
SVI5  at  stations  B  and  C.  The  upper  traces  (a)  are  the  superposed  measured  single  shot  (SMSS) 
signals  (S  VI4).  The  middle  traces  (b)  are  the  observed  production  shot  seismograms.  The  lower 
traces  (c)  are  obtained  by  convolving  the  source  time  function  with  the  Green’s  functions 
calculated  for  each  shot-receiver  distance  and  linearly  superposing  the  resulting  seismograms 
with  alternating  20  ms  delay  times. 

Figure  10.  The  spectral  amplimdes  of  the  vertical  component  at  station  B  for  SVI5  (left)  and 
radial  component  at  station  D  for  SVI6  (right)  are  plotted  in  the  middle  of  the  figure  for  the 
measured  production  shot  seismograms.  Spectra  of  the  corresponding  calculated,  superposed 
seismograms  and  the  SMSS  signals  are  given  in  the  upper  and  lower  traces,  respectively,  and  are 
shifted  upwards  by  two  decades  for  easier  viewbg. 

Figure  A 1.  The  unit  impulses  on  the  lower  right  of  the  figure  are  superposed  and  delayed  in  time 
by  32  ms  including  the  delays  for  travel  times  for  compressional  waves  (4.5  km/s)  in  a  whole 
space.  In  the  upper  four  traces,  normally-distributed  random  numbers  are  added  to  the  firing 
times  with  maximum  variance  levels  from  1%  to  6%  of  the  desired  firing  times.  The 
corresponding  amplitude  spectra  are  plotted  in  the  left  diagram  on  a  log-log  scale  from  5  to 
400  Hz.  All  amplitudes  are  normalized  to  the  maximum. 
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SV15:  20  ms  delay  (20  holes) 
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April  10, 1992 


Dear  Redpienc 

Enclosed  is  a  report  entitled  "Physical  models  of  spall  zone  ground  motions  and  the 
detennination  of  spatial  decay  rates"  by  B.W.  Smmp  and  T.A.  Weaver  prepared  for  the 
Source  Region  Program  at  Los  Alamos  National  Laboratory.  The  report  discusses  results 
from  a  preliminary  study  designed  to  document  relarionships  between  physical  processes  in 
the  spall  zone  around  explosions  and  the  radiated  seismic  waves.  This  is  imponant  for 
developing  realistic  spall  models  that  can  be  used  to  predict  effects  on  far-field  verification 
signals.  A  simple  physical  model  for  spall  is  proposed  the  includes  free-field  attenuation 
(decay)  of  the  wavefield  coupled  with  free-surfacc  interactions.  The  model  suggests  that 
spall  data  should  be  ploned  against  fixe-surface  range  and  interpreted  in  terms  of  two  decay 
rates.  The  decay  of  data  at  close  ranges  (out  to  about  100  m/lct*/^)  is  controlled  by  free- 
surface  effects  while  data  beyond  is  dominated  by  the  free-field  interaction.  It  is  illustrated 
that  decay  rates  (used  in  quantificarion  of  spall  momentum)  can  be  biased  if  an  improper 
physical  model  is  used.  The  work  contained  in  the  report  was  funded  by  the  Source 
Region  Program  at  the  Los  Alamos  National  Laboratory  for  the  DOE  Office  of  Arms 
Control  and  Nonproliferation  and  DARPA,  F19628-89-K-0025,  as  monitored  by  the 
Phillips  Laboratory.  This  report  is  being  distributed  to  researchers  at  the  DOE  laboratories 
and  other  government  agencies  and  univenities. 


ABSTRACT 


Spall,  the  tensile  failure  of  near-surface  layers,  which  is  observed  above 
contained  explosions,  has  been  Identified  as  a  possible  secondary  seismic 
source  contributing  to  teleseismic  and  regional'  signals.  The  relative 
importance  of  this  secondary  source  can  be  constrained  if  the  motion  field  in 
the  spall  zone  is  characterized.  Spall  zone  motions  from  nuclear  explosions 
detonated  above  the  water  table  at  Pahute  Mesa  are  analyzed  to  develop  these 
models.  Acceleration,  velocity,  displacement,  and  dwell  time  measurements  are 
made  from  gauges  placed  directly  above  the  explosion,  most  often  at  the  free 
surface.  Decay  of  peak  motions  are  strongly  affected  by  the  free  surface  with 
little  change  in  amplitude  out  to  a  free  surface  range  of  100  m/kf/3  followed 
by  rapid  decay  beyond.  Free  surface  interactions  are  assessed  with  first-order 
elastic  spherical  wave  calculations  that  match  observed  peak  velocity  decays. 
These  results  indicate  that  the  spall  zone  motions  may  be  strongly  affected  by 
the  scaled  depth  of  burial  of  the  explosion.  Spall  zone  velocities,  displacements 
and  dwell  times  are  compared  for  consistency  with  a  gravitational  model.  The 
data  is  in  agreement  with  the  functional  form  of  theoretical  models  although 
observed  displacements  may  be  as  much  as  a  factor  of  two  to  four  greater  than 
the  model  predicts  for  observed  velocities  and  dwell  times.  These  differences 
may  reflect  the  continuous  nature  of  the  spall  process  and/or  the  role  of 
material  strength  in  these  phenomena. 
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